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Resumen 



Los procesos de decoherencia y disipacion en sistemas mesoscopicos han sido 
objeto cn las ultimas decadas de gran interes dado su papel clave en el estudio y 
utilizacion de fenomenos cuanticos a escalas cada vez mayores. Los espectaculares 
logros alcanzados en las tecnicas de fabricacion y medicion han permitido el diseiio 
de nanoestructuras moviles que habilmente combinadas con otros protagonistas de la 
historia fisica reciente, tales como puntos cuanticos o gases electronicos de baja di- 
mensionalidad, otorgan la fascinante oportunidad de estudiar la interaccion electron- 
fonon cuanto a cuanto, con un potencial tecnologico prometedor. En la mayoria de las 
aplicaciones vislumbradas, asi como para la observacion de comportamiento cuantico 
en variables mecanicas macroscopicas tales como el centro de masas de resonadores 
micrometricos, es imprescindible una minimizacion de los procesos disipativos que per- 
turban la dinamica vibracional de estos sistemas. En consecucncia sc csta llcvando a 
cabo un gran esfuerzo actualmente en el analisis y modelizacion de los mismos, tanto 
desde un punto de vista experimental como teorico. 

Esta tesis forma parte de dicho esfuerzo. En concreto, dos de las tres partes 
en que se divide estan consagradas a la modelizacion de algunos de los principales 

mecanismos de friccion presentcs cn dos clascs de sistemas nanoclcctromecanicos 
(NEMS): nanoresonadores construidos a partir de heteroestructuras de semiconduc- 
tor, y nanoresonadores en los que la parte movil es un compuesto de carbono de baja 
dimensionalidad, bien sea una lamina de grafeno o un nanotubo. 

En el primer caso nos hemos centrado en los procesos superficiales de absorcion 
de energi'a mecanica asociados al acabado imperfecto, la presencia de impurezas y la 
estructura superficial parcialmentc dcsordcnada, invocando las similitudes existentes 
con Ids procesos dc atcnuamiento dc ondas aciisticas cn materiales amorfos. El re- 
sultado cs un modelo de espines acoplados al conjunto de modos vibracionales del 
resonador, analizado por medio de tecnicas desarroUadas para el estudio del Spin- 
Boson Model y de procesos de relajacion. Asimismo se modelizan los posibles efectos 
negatives que conlleva la deposicion de electrodos metalicos sobre las estructuras de 
semiconductor. 

Los osciladores de compuestos de carbono presentan una serie de peculiaridadcs 
que los diferencian de sus homologos semiconductores, destacando su alto grado de 
cristalinidad y caractcr mctalico o semimetalico. EUo implica la prevalencia de mecan- 
ismos de friccion distintos, varios de los cuales se modelizan, destacando la excitacion 
de pares electron-hueco en el resonador debido a cargas distribuidas por la estructura, 
estudiada con metodos perturbativos autoconsistentes. 
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La tercera y ultima parte se dedica a otro tipo de sistema mesoscopico, las 
particulas nanometalicas, cuya respuesta optica viene dominada por la excitacion 
colectiva denominada plasmon de superficie. La presencia en el espectro de excita- 
ciones electronicas de pares electron-hueco acoplados al plasmon provoca su progre- 
sivo atenuamiento tras ser excitado inicialmente por un campo electrico externo, por 
ejemplo un laser. Este proceso se puede consider ar como ejemplo de entorno disipa- 
tivo con un numero finito de grados de libertad en interaccion con el subsistema de 
interes, el plasmon. Los efectos que esta finitud conlleva en la modelizacion teorica 
de la dinamica electronica son analizados en detalle, justificando la validez de ciertas 
aproximaciones asumidas hasta ahora sin demostracion previa. 

Conclusiones. 

La friccion a bajas temperaturas asociada a los distintos procesos superficiales 
en resonadores de semiconductor ha sido modelizada adaptando el Standard Tun- 
neling Model utilizado para explicar las propiedades acusticas de solidos amorfos a 
~ 0.01 < T < 1 K. El orden de magnitud asi como la debil dependencia con la 
temperatura observada en los experimentos es cualitativamente reproducida. Sin em- 
bargo, llegar a un acuerdo cuantitativo se antoja dificil hasta que no se posea un 
conocimiento mas detallado de los procesos y defectos estructurales presentes en el 
resonador. 

Diversos mecanismos disipativos presentes en nanoresonadores de grafeno y nan- 
otubos han sido analizados, obteniendo, al igual que en el caso previo, dependencias 
parametricas de los mismos en funcion de las variables fisicas relevantes del sistema, 
como la temperatura y dimensiones caracteristicas. La importancia relativa de los 
mismos en funcion del regimen (amplitud de vibracion, temperatura) ha sido estable- 
cida, predominando a bajas temperaturas la disipacion ohmica debido a excitaciones 
electronicas en el seno del resonador y a altas temperaturas la atenuacion por proce- 
sos termoelasticos. 

Se ha presentado un modelo muy utilizado en el estudio de la dinamica de plas- 

mones superficiales y de la respuesta optica de clusters metalicos, analizando la co- 
herencia interna de dicho esquema teorico, justificando varios puntos clave del mismo 
que hasta ahora se habian tomado como validos sin justificacion rigurosa. De dicho 
analisis se obtiene asimismo informacion sobre los tiempos caracteristicos de evolucion 
del plasmon y del resto de excitaciones electronicas acopladas al mismo, lo cual per- 
mite fundamentar una aproximacion Markoviana a la hora de estudiar la dinamica 
del sistema acoplado plasmon - pares electron-hueco. 
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The physics of decoherence and dissipation in mesoscopic systems has attracted 
great interest, given its key role in the study and use of quantum phenomena at ever 
larger scales. The spectacular achievements in fabrication and measurement tech- 
niques have opened the doors to the design of mobile nanostructures which, combined 
in smart ways with other recent cornerstones like quantum dots or low-dimensional 
electron gases, grant the fascinating opportunity to study electron-phonon interac- 
tions at an individual level, with a promising technological potential. In most of the 
envisioned applications, as well as to observe quantum behavior in macroscopic me- 
chanical variables like the center-of-mass of micrometric resonators, a minimization 
of the dissipative processes perturbing the vibrational dynamics of these systems is 
compulsory. Consequently strong efforts are being made to analyze and model them, 
both from the experimental as well as from the theoretical sides. 

This thesis exemplifies these efforts. Specifically, two out of the three parts con- 
stituting it are consecrated to the modeling of some of the main friction mechanisms 
found in two kinds of nanoelectromechanical systems (NEMS): nanoresonators built 
from semiconductor heterostructures, and nanoresonators whose mobile part is a low- 
dimensional carbon compound, be it graphene or a nanotube. 

In the first case we have focussed in the surface processes absorbing mechanical 
energy, related to the imperfect finish, presence of impurities and partially disordered 
surface structure. Similarities with acoustic wave damping in amorphous solids have 
been invoked to build an adaptation of the Standard Tunneling Model which aims 
to provide a low temperature description of such processes. The result is a model of 
spins coupled to the ensemble of vibrational eigenmodes of the resonator, which is an- 
alyzed by means of techniques developed for the study of the Spin-Boson Model and 
relaxation processes. Possible negative effects due to metallic electrodes deposited on 
top of the semiconductor structures are also modeled. 

Oscillators made from carbon compounds, demonstrated for the first time this 
year, display a series of peculiarities distinguishing them from their semiconductor 
counterparts, specially their high degree of crystallinity and semimetal or metal char- 
acter. The prevailing friction mechanisms thus differ, several out of which are mod- 
eled, with the excitation of electron-hole pairs in the resonator due to charges dis- 
tributed throughout the device dominating at low temperatures and small vibrational 
amplitudes. 

The third and last part is devoted to another kind of mesoscopic system, namely 
nanometallic clusters, whose optical response is dominated by the so-called surface 
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plasmon collective excitation. The presence in the electronic excitation spectrum of 
one-body electron-hole pairs coupled to the plasmon causes its progressive damping 
after its initial excitation by an external electric field, e.g. a laser pulse. This process 
may be considered as an example of finite dissipative environment interacting with the 
subsystem subject of our interest, the plasmon. The consequences of this finiteness on 
the theoretical modeling of the electronic dynamics are analyzed in detail, justifying 
the validity of several approximations assumed until now to hold without proof. This 
study will also shed some light on the characteristic times of the different electronic 
degrees of freedom, providing a basis for the use of a Markovian approximation in the 
analysis of the plasmon dynamics of nanometric clusters. 
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Chapter 1 
NEMS 



1.1 Introduction: current research topics 

Nanoclcctromcchanical systems (NEMS) constitute a recent fascinating subject 
of current research, both from fundamental and apphcation-oriented points of view 
[CraOO, Cle02, Ble04, EROS, SR05]. Inheritors of the bigger microelectromechanical 
systems (MEMS), they are mesoscopic sohd-state devices with a special ingredient: 
one resonating component, with some of its characteristic dimensions lying in the 
submicron regime, whose mechanical motion is integrated into an electrical circuit. 
Typically this moving element can be a bridge, a pillar, or a cantilever, one of whose 
mechanical eigenmodes is excited and measured. Two features of these eigenmodes 
are essential to explain the potential of NEMS: First, due to its small dimensions, 
their modes vibrate at very high frequencies cjq) currently from tens of MHz to GHz 
(microwave frequencies). This is desirable if one wishes to study the quantum regime 
where thermal fluctuations are smaller than the energy of a phonon, kT < hwo (for 
GHz frequencies, milliKelvin suffice), and also for several high-speed signal processing 
components [Bh05, HFZ+05, BQKM07] . 

Second, and perhaps more important, they exhibit extremely low damping, mea- 
sured in terms of the so-called quality factor 

where Auq is the measured linewidth, reaching often Q ~ 10^ or higher. This trans- 
lates, for example in the case of NEMS-based single-electron transistors, into power 
dissipation rates several orders of magnitude below the ones of common CMOS- 
based transistors found in standard microprocessors operating at similar frequen- 
cies [BQKM07]. Thus, NEMS-based computing technology constitutes a very attrac- 
tive alternative for applications demanding GHz or lower clock-speeds, robustness 
and/or operation at high temperatures {T > 200°C). 
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Chapter 1: NEMS 



Apart from nanoscale versions of Babbage's 18th-century mechanical computers, 
these two features open the scope for many other technological applications, among 
whom ultrasensitive sensors and actuators clearly outstand. Their working princi- 
ple is easy: The interaction of the resonator with other elements of the circuitry, 
like charges if the vibrating element is charged, or with elements of the environment 
outside the device, cause measurable shifts in the frequency of the eigenmode. The 
sensitivity depends crucially on the small dimensions of the device and the linewidth 
AciJo of the mode, currently achieving sub-attonewton force detection [MROl] and 
mass sensing of individual molecules [KLER04]. Other successes include single-spin 
detection [RBMC04], study of Casimir forces at the nanometer scale [DLC^OS], high- 
precision thermometry [HKC^07] or in-vitro single-molecule biomolecular recogni- 
tion [DKEM06]. 

But leaving aside applications, they are also proving as very powerful tools in 
fundamental research: 

• They allow for a study of electron-phonon interactions at a single quantum 
level, with setups prepared for the study of single electrons interacting with 
single phonons, or the analysis of single electrons shuttled via mechanical motion 
[EWZBOl]. 

• At low enough temperatures several manifestations of a quantum behavior of the 
mechanical oscillator have been observed, or are expected to appear, like [Ble04] : 

— Zero-point motion detection and back-action forces due to the measure- 
ment process [LBCS04, NBL+06]. A displacement sensitivity of 10~^Ais 
needed, attainable with current cantilevers. 

— Avoiding back-action forces, quantum non-demolition measurements of en- 
ergy eigenstates [SDC04, JLB07,MZ07, Jac07]. 

— Quantization of thermal conductance [Ble99, SHWROO]. 

— Quantum tunneling between macroscopically distinct mechanical states, 
with several proposals for applications in quantum computing [CLWOl, 



GC05,SHN06,XWS+06], see fig. (1.1) 



Preparation of minimum-uncertainty quantum squeezed states with a po- 
sition uncertainty below that of zero-point fluctuations [RSK05a,RSK05b]. 



Figure (1.2) shows some beautiful examples of setups for some of the applications 
mentioned above. 

The study of quantum superpositions of mechanical degrees of freedom in NEMS 
is also a straight way to assess two important issues: i) The validity of quantum me- 
chanics for ever larger systems [Leg02] , and ii) The theoretical framework explaining 
environmentally-induced decoherence [LCD+87, Wei99], a part of which we will use 



1.2 Semiconductor NEMS fabrication 
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Figure 1.1: Preparation of quantum superpositions of macroscopically distinct res- 
onator states, a) From [ER05]: the cantilever interacts with a two-level system (TLS) 
in such a way that the flexure depends on the TLS state, | |) or | |). If the TLS is 
prepared in a superposition state, the mechanical resonator state will be a linear com- 
bination of the two flexure states, with a probability distribution P{x, t) as depicted 
to the right, b) When a doubly clamped bar is sufficiently compressed, it buckles, with 
two possible buckling positions. At low enough temperatures, tunneling through the 
energy barrier separating both states gives rise to eigenstates corresponding to super- 
positions of left and right buckled states, enabling its use as a qubit. From [SHN06] 



in this thesis to analyze the mechanisms limiting the quality factors Q of the eigen- 
modes, and with them the observation of such quantum behavior. 

As for most of these technological applications and fundamental science works 
low-dissipation, high-Q resonators are desirable, several works have been devoted 
to the analysis of the different sources of friction present in MEMS and NEMS 
[CR99,YSC+00,EOS+00,CR02,YOE02,MHE+02,AM03,HHP+03,ZGS+05,FZMR06, 
SCK"'"07], trying to determine the dominant damping mechanisms and ways to min- 
imize them. In the first two parts of this thesis we will analyze them, both for the 
most common semiconductor-based resonators, with vibrating structures made of Si, 
GaAs, SiN, SiC... (part I), and also for carbon-materials based resonators, like nan- 
otube [SYU^04] or graphene [BvdZV+07] devices (part II). But before we embark 
ourselves in such a task, it is convenient to learn a few basic facts about NEMS. 



1.2 Semiconductor NEMS fabrication 

The performance of a given nanomechanical oscillator depends to a great extent 
on its constituent materials and the manufacturing steps they have gone through, not 
to speak about reproducibility and scalability of results. A good review on MEMS 
and NEMS manufacturing techniques can be found in [Cle02]. We will focus on the 
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Figure 1.2: Several examples of nanoelectromechanical devices (NEMS). a) Mechan- 
ical single electron shuttle [EWZBOl]. b) Ultrasensitive magnetic force detector used 
to detect a single electron spin [RBMC04] . c) Setup using a torsional oscillator used to 
study Casimir forces to check the validity of Newtonian gravitation at the nanometer 
scale [DLC"'"05]. d) Parametric radio-frequency amplifier providing a thousandfold 
boost of signal displacement at 17 MHz [SR05]. e) A 116-MHz nanomechanical res- 
onator coupled to a single-electron transistor [KC03]. f) Tunable carbon nanotube 
resonator operating at 3-300 MHz [SYU^04]. 



most widely used technique for building NEMS, namely Epitaxial Heterostructure 
Nanomachining, applied to a (Si-Si02-Si) heterostructure, chosen frequently by ex- 
perimentalists due to the high degree of crystallinity and perfection of the resulting 
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resonators, and more concretely in a study of dissipation with which we will compare 
our predictions [ZGS+05, SGCOTa, SGC07b, SCK+07]. 

The starting heterostructure is fabricated by oxygen ion implantation in a crys- 
talline Si wafer, followed by a high-temperature annealing of several hours, and is 
commonly known as SIMOX (Separation by IMplantation of OXigen) wafer. During 
the second step the previously formed mixed buried layer Si-Si02 becomes purified 
through oxygen migration, constituting a separated Si02 smooth buried layer, below 
a single-crystal recrystallized Si top layer. Typically, this single- crystal top layer, out 
of which our resonator will be carved, contains a relatively large number of defects, 
mostly threading dislocations running from the Si02 layer to this one, with typical 
defect densities of about 10^ — 10^/ cm^. The top layer presents variations in thickness 
~ 10 nm, as does the silicon oxide underneath. 

To create the suspended structure, the wafer is subject to the nanomachining pro- 



cess depicted in fig. (1.3): first an etch mask is defined using lithography: a PMMA 
photoresist is coated over the wafer and then patterned using a nm-resolution electron 
beam, constituting a metal liftoff mask. Afterwards a layer resistant to Reactive Ion 
Etching (RIE) is deposited on top of it and lifted off, leaving behind a RIE-resistant 



mask, fig. (1.3 d). Then the RIE is applied on the sample enough time to remove 
the top Si layer except where the protective mask lies, and remove also part of the 
underlying Si02 layer, fig.( 1.3[ :;). The suspended structure appears when the device 



is submerged in a selective wet etch which removes the Si02, like hydrofluoric acid. 



fig. ( 1.3 i). Depending on the method to be used for the driving of the resonator, a 
top metallic coating can be deposited at will for example using lithography, before 
the RIE mask is created. 

Even though spectacular, the final result is far from perfect: the SIMOX wafer 
defects, together with surface roughness and imperfections due to the resolution lim- 
itations of the different steps result in strong dissipation, as we will later see. 



1.3 Motion transduction of mechanical eigenmodes 
in semiconductor NEMS 

Well stablished transduction setups for the detection and driving of mechani- 
cal modes in MEMS, like optical methods as fiber-optic interferometry and reflec- 
tion, or piezoresistive displacement transduction, cease to work properly for the 
submicrometer-sized NEMS. In the former case the light wavelengths exceed the res- 
onator transverse dimensions and diffraction effects spoil the scheme, while in the 
latter dissipation heats the sample, precluding its use for low-temperature applica- 
tions and studies. 

For semiconductor NEMS, two other schemes have proved quite successful: 
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a) SIMOX substrate d) Oxide selective wet etch 
Si . /<5 



b) Etch mask definition c) Anisotropic Si etch 




Figure 1.3: From [Cle02]. Building a nanoresonator using Epitaxial Heterostructure 
Nanomacliining. (a) The process starts with a previously fabricated heterostructure, 
such as the SIMOX (Separation by IMplantation of OXigen) Si-Si02-Si shown, (b) 
Using electron beam lithography a mask is etched, (c) Anisotropic dry etch such as 
Reactive Ion Etching, (d) The oxide layer is removed using a selective Si02 wet etch, 
obtaining in this way the suspended structure. Metallization of the top layer can be 
performed lithographically before the process. 



1.3.1 Magnetomotive scheme 



As shown in fig.( |1.4^ ), in this setup [MHE+02] the nanobridge is covered by a thin 

metallic layer, connected in turn to an external circuit whose AC current is controlled 
by the experimentalist, fixing its frequency uj at will. The whole circuit is embedded 
in a strong static magnetic field of several Tesla, perpendicular to the main axis of 
the bridge, as indicated in the figure. When the electrons of the AC current flow back 
and forth at frequency a;, they experience a Lorentz force due to the presence of the 
magnetic field perpendicular both to the field and the main axis, forcing the bridge 
to oscillate up and down. The area enclosed by the circuit is thus a function of time, 
and so the magnetic flux through it, therefore originating an induced electromotive 
force opposing the motion of the bridge, which is detected by a network analyzer. 
When the frequency is tuned to the one of a given mechanical eigenmode, the motion 
is strongly enhanced. With this scheme many modes can be excited, just by changing 
the circuit's frequency to the one of the mode of interest. The presence of the top 
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metallic layer adds extra dissipative sources which will be analyzed later. 

1.3.2 Coupling to SETs and other sensitive electronic probes 

The magnetomotive driving can be used in combination with more elements close 
to the resonator, to study electron-phonon interactions in a previously unexplored 
regime, with a single-electron single-phonon control. This can be achieved by cou- 
pling our magnetomotively driven charged resonator to a mesoscopic device whose 
conductance can be controlled and measured with great accuracy, like quantum dots, 
quantum point contacts or single-electron transistors (SETs) [KC03], as schematically 
depicted in fig.( 1.4[ 3). Using similar ideas but different geometries, fascinating devices 



like single-electron mechanical shuttles have been demonstrated [EWZBOl, SB04]. 

1.4 Nanoscale elasticity theory: Vibrational eigen- 
modes of long and thin rods 

Regarding the elastic properties of nanoresonators, it is clear that as size shrinks a 
point is reached when classical continuum mechanics ceases to be a good approxima- 
tion. Nevertheless, for the sizes of current semiconductor nanorods and wavelengths 
of the modes excited (~ 50nm - lO/xm), much bigger than typical interatomic spac- 
ing, continuum approximations continue to hold [SB98,LB99,SB02,SBN102,LBT"'"03, 
Ble04], and classical elasticity theory, using bulk concepts like Young's modulus E 
or Poisson's coefficient predicts with reasonable accuracy the eigenmodes observed 



in experiments. In Appendix A.l a brief review of continuum elasticity concepts and 
the derivation of the equations for the different modes present in quasi-lD resonators 
is given. Here we state the main results. For a long and thin rod, considering small 
deviations from equilibrium and wavelengths bigger than the transversal dimensions, 
three kinds of vibrational modes are found, namely longitudinal, fiexural and tor- 
sional. 

Longitudinal (or compressional) modes those found in bulk 3D solids 

but occur only along the direction of the main axis, chosen to be the Z axis, with 
the Z-component Uz of the local displacement field Ui = x[ — Xi^eq obeying the wave 
equation 

d'^Uz p d'^Uz 



dz"^ E dt^ 



1.2 



p being the mass density. Torsional vibrations, like the one of fig.( 1.5[ 3), are parametrized 



in terms of the twisting angle (j){z), which obeys the wave equation 
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Figure 1.4: (a) Scheme of the magnetomotive driving and detection setup with a 
nanobridge. The electrons flowing back and forth through the bridge experience a 
Lorentz force which drives the motion of the bridge. Due to the oscillation the area 
of the loop formed by the circuit varies, and with it the magnetic flux enclosed by it, 
inducing an electromotive force according to Faraday's law, measured as a function of 
the AC electron current frequency by a vector network analyzer [MHE"'"02]. (b) The 
vibrating bridge interacts with the electrons flowing through a mesoscopic detector 
such as a quantum point contact or quantum dot placed nearby, altering the current 
Ids-, which is monitored with extremely high precision. Figures from [ER05,SR05]. 



where C is the so-called torsional rigidity. Bending (flexural) modes differ from the 
previous two in that they do not obey a standard wave equation. They correspond 
to transversal displacements X{z,t),Y{z,t) following the equations of motion 



EI.^ - -PS^ . (1.4) 

Here Ij are the inertia moments of the X and Y axes, and S is the cross-section area 
of the rod. Even though these are not usual wave equations they admit plane wave- 
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Figure 1.5: Examples of cantilever modes not present in a bulk solid, with the black 
lines indicating the beam in absence of movement serving as guide to the eye. a) 
Second flexural (bending) mode, b) Fundamental (first) torsional (twisting) mode. 



kind of solutions X{z,t),Y{z,t) ~ e**^'^^ "^'^ but with a quadratic dispersion relation, 
^ji^) — a/-E/j/ {pS)k'^ . An example of bending mode for a cantilever can be seen in 

fig.(|tgi). 



1.5 Dissipative mechanisms damping the vibrational 
eigenmodes of semiconductor nanoresonators 

Numerous dissipative sources limit the lifetime of an externally excited mechanical 
eigenmode, absorbing its energy irreversibly. Some of them, like surface impurities, 
clamping imperfections or gas friction can be reduced devising solutions like working 
in high vacuum conditions [MHE"''02,MH07], annealing the samples at the right stage 
of the fabrication process [YSC"'"00, YOE02, WOE04,LVS"'"05] or engineering carefully 
the clamping points to minimize the coupling to bulk vibrational modes outside the 
resonator [HFZ"^05]. They are thus extrinsic mechanisms, as opposed to intrinsic 
mechanisms, present even in a theoretically perfect sample, such as anharmonic cou- 
pling with other modes, electron-phonon or electron-electron interactions. A brief 
description of the most important mechanisms follows. 



1.5.1 Extrinsic mechanisms 

Gas or liquid friction 

Collisions with surrounding gas molecules can impose severe limitations on the 
quality factor, specially in those setups designed to work in ambient conditions or 
in a liquid environment [BW04, DKEM06, MH07]. Depending on the fluid density 
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several regimes appear which have to be analyzed with different frameworks, ranging 
from free molecular to continuum limits. Q increases with decreasing pressure P, 
and for P < 1 mTorr it can be safely disregarded as the main frictional force for 
nanoresonators with Q > 10^ [MHE+02, ER05]. 

Heating linked to the actuation or detection setup 

In the magnetomotive scheme the flow of electrons through the metallic top layer 
is not baUistic, and heating of the layer is generated through their inelastic collisions, 
which in turn heats the semiconductor substrate underneath. Fortunately for the 
usual current densities this effect seems to be irrelevant compared to others [MHE"'"02]. 
In schemes using actuation and detection lasers illuminating the resonator, heating 
may play a more significant role, but we will not study it further in this thesis. 

Noise introduced by the external electrical circuit 

The external circuitry can be a source of extra damping and frequency shift, due 
to its finite impedance. This effect can be used for applications requiring an external 
control of linewidth, hke signal processing [CR99, Sch02]. 

Clamping (attachment) losses 

This mechanism corresponds to the transfer of energy from the resonator mode 

to acoustic modes at the contacts and beyond to the substrate through the anchor- 
age areas [J168,PJ04]. The motion of the resonator forces as well a motion of the 
atoms linking it to the rest of the device, irradiating in this way elastic energy to its 
surroundings. This mechanism can be of special relevance for short and thick beams, 
while for long and thin ones it plays a secondary role [HFZ+05, FZMR06]. 

Surface roughness and imperfections 

As mentioned previously, the structures obtained through current methods present 
a certain degree of roughness of ~ 10 nm. On top of this, adsorbed impurities, dis- 
locations and surface reconstruction processes are all coupled to the local atomic 
displacements linked to the resonator's motion, absorbing energy from the vibra- 
tional eigenmode and transmitting it irreversibly to the rest of degrees of freedom of 
the system. These imperfections and impurities also exist in the core of the resonator 
structure, but to a lesser degree as in the case of its surface, exposed to the envi- 
ronment and with atoms less bound than in the bulk, more susceptible to thermal 
rearrangements. Strong evidence has been accumulated of the increasing role played 
by surfaces as the resonator sizes shrink and surface-to- volume ratios grow, being the 
dominating friction source in thin nanobeams at low temperatures: i) The quality 
factor Q roughly scales hnearly with the former ratio [YSC+00, EROS] , as shown in 
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Figure 1.6: From ref. [ER05]. Evolution of reported quality factors in mono crystalline 
mechanical resonators with size, showing a decrease with linear dimension, i.e., with 
increasing surface-to- volume ratio, indicating a dominant role of surface-related losses. 



fig. (1.6), ii) Q has been observed to increase up to an order of magnitude when the 
system went through a thermal treatment [YOEOO, WOE04, LVS+05] or the surface 
was deoxidized [YOE02]. In the next chapter we will study in depth these surface- 
related losses, dominant for thin beams at low temperatures, precisely the interesting 
regime where quantum effects show up if a very high Q is reached. Surface waves 
can be also excited by the vibrational eigenmode due to the roughness, but at low 
temperatures this process becomes strongly suppressed [MHE+02], and thus will be 
ignored in our study. 



Generation of e-h pairs in the metallic electrode due to electron-phonon 
coupling 

In the magnetomotive scheme, where the semiconductor bridge is covered with 
a thin metallic top layer, the flowing electrons feel the Coulomb potential generated 
by static charges in the device and substrate. When the oscillating part is set into 
motion, this potential is time-dependent, and gives the chance to the electrons to 
absorb part of the mechanical energy stored in the eigenmode, creating e-h pairs. We 
will study this mechanism in the next chapter. 
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Temperature (K) 



Figure 1.7: From ref. [LROO]. Minimum dissipation imposed by the thermoe- 
lastic damping mechanism, as a function of temperature, for any resonator made from 
GaAs or Si. 



1.5.2 Intrinsic mechanisms 



Thermoelastic relaxation 



The strain field generated by the excited eigenmode modifies the local thermal 
equilibrium, with the coupling to the local temperature field characterized by the 
linear expansion coefficient a = {l/L)dL/dT. This translates into the local creation 
of thermal phonons by the long-wavelength eigenmode, which then diffuse to points 
with different temperatures leading to dissipation [Zen38, LROO, Cle02, U06]. This 
mechanism is strongly temperature dependent, being suppressed at low temperatures 
and dominating at high ones. In [LROO] a detailed study of this mechanism is carried 
out and the limits to the quality factor were established for any resonator made with 
Si or GaAs, two of the commonest materials, see fig.(1.7). 
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Anharmonic mode coupling (phonon-phonon interactions) 

When the nanoresonator is strongly driven the harmonic approximation for the 
lattice potential breaks down, and higher order terms which couple the former vibra- 
tional eigenmodes have to be included, giving a finite lifetime to a given externally 
excited mode. At a temperature T, the externally excited mode will interact with 
thermal phonons, loosing its energy [Cle02, U06]. We will assume that the experi- 
ments are done at low excitations amplitudes, within the linear regime, where this 
effect can be neglected. Nevertheless, it is worth mentioning the exciting physics 
and potential applications of strongly driven nanoresonators, currently a subject of 
intensive research [KRSL07, SDC04, URA04, BM05, SHN06, SIM07]. 

Surface reconstruction and amorphization 

Even if the fabrication processes and materials used were perfect, purely ther- 
modynamic considerations imply that at a temperature T > the sharp edges of 
a finite crystal begin to round off due to the increasing importance of entropy, in a 
process called creation of vicinal surfaces [Wei99]. The associated rearrangement of 
atoms and coupling of vibrational modes limit the maximum quality factor attainable. 

After this brief overview of semiconductor nanoresonators and friction mecha- 
nisms affecting them, we will focus our attention in Chapter |2] on the theoretical 
modeling of surface-related dissipation at low temperatures, borrowing ideas from 
sound attenuation in amorphous glasses. 
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Chapter 2 



Surface dissipation in 
semiconductor NEMS at low 
temperatures 



2.1 Introduction 

In this chapter we will try to provide a unified theoretical framework to describe 
the processes taking place at the surface of nanoresonators at low temperatures, 
which have been observed to dominate dissipation of their vibrational eigenmodes in 



this regime, see fig.(2.1). This is certainly a challenging task, since as we have seen 
many different dynamical processes and actors come into play (excitation of adsorbed 
molecules, movement of lattice defects or configurational rearrangements), some of 
whom are not yet well characterized, so simplifications need to be made. We will 
provide such a scheme, based on the following considerations: 

• Experimental observations indicate that surfaces of otherwise mono crystalline 
resonators acquire a certain degree of roughness, impurities and disorder, re- 
sembling an amorphous structure [LTBP99, Wei99]. 

• In amorphous solids the damping of acoustic waves at low temperatures is suc- 
cessfully explained by the Standard Tunneling Model [AHV72, Phi72, Phi87, 
Esq98], which couples the acoustic phonons to a set of Two-Level Systems 
(TLSs) representing the low-energy spectrum of all the degrees of freedom able 
to exchange energy with the strain field associated to the vibration. 

We will describe the attenuation of vibrations in nanoresonators due to their 
amorphous-like surfaces in terms of an adequate adaptation of the Standard Tunnel- 
ing Model. We start reviewing such model, and some concepts and techniques used to 
analyze the spin-boson hamiltonian appearing in it. Then we proceed to apply in the 
proper way the Standard Tunneling Model to the non-bulk nanoresonator geometry. 
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Figure 2.1: Sketch of the systems whose surface dissipation will be studied in this 
chapter. The inset shows a doubly clamped beam, while the main figure shows a 
cantilever characterized by its dimensions, width (ly), thickness (t), and length, (L), 
where w ^ t <^ L. The height above the substrate is h. A schematic view of 
the surface is given, highlighting imperfections like roughness and adsorbates, which 
dominate dissipation at low temperatures. 



Expressions for the frequency shift and dissipation in terms of the inverse quality 
factor for flexural and torsional modes will be given, showing the scaling with 
dimensions, temperature and other relevant parameters of these systems. We will 
see that qualitative agreement with experimental observations is obtained, discussing 
limitations and possible modifications of the scheme to reach quantitative fitting to 
experiments. For NEMS covered with metallic electrodes the friction due to electro- 
static interaction between the flowing electrons and static charges in the device and 
substrate will be also studied. 



2.2 Damping of acoustic waves in amorphous solids: 
Standard Tunneling Model 

The description and understanding of the properties of amorphous solids, where 
most of the simplifications linked to an infinite periodic structure are absent, has been 
full of unexpected results [Zal98]. In the low temperature regime, pioneering experi- 
ments by Zeller and Pohl [ZP71] revealed at T < 1 K striking differences between the 
thermal properties of dielectric amorphous solids as compared to those of crystalline 
solids. In a Debye dielectric crystalline solid it is well known that the heat capacity 
and the thermal conductivity due to phonons behave as C ~ T'^ and k ~ T^, but 
in the case of amorphous dielectrics none of them is true, instead a C ~ T^'^ (and 
much bigger in value as compared to a crystalline form of the material) and n ~ T^'^ 
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behavior is observed. Apart from this, a heat release as a function of time is ob- 
served on long time scales: after a sample is heated and then rapidly cooled down 
to some temperature Tq, in adiabatic conditions the sample begins gradually to heat 
itself afterwards. Moreover, measuring their acoustic properties at low temperatures, 
which a priori were expected to be similar to those found in crystalline solids due 
to the long-wavelength vibrations involved, the following dependence on the acoustic 
intensity was obtained [HAS"'"72]: at low intensities, as the temperature is lowered 
the mean free path of the phonon decreases, while at high intensities a monotonous 
increase is observed. 



Several theoretical scenarios were developed to explain these low-temperature fea- 
tures, the most successful of whom was the so-called Standard Tunneling Model, 
proposed simultaneously by Phillips and Anderson et al. [AHV72, Phi72]. 

The tunneling model claims the existence of the following states, intrinsic to 
glasses: an impurity, an atom or cluster of atoms within the disordered structure 
which present in their configurational space two energy minima separated by an en- 
ergy barrier (similarly to the dextro/levo configurations of the ammonia molecule), 
as depicted schematically in fig.( |2.2^ ). These states can be modeled as a degree of 
freedom tunneling between two potential wells. At low temperatures only the two 
lowest eigenstates have to be considered, characterized by the bias Ag between the 



wells and the tunneling rate Ag through the barrier, see fig. (2. 2d). This corresponds 
to the Two-Level System (TLS) hamiltonian: 

i/o = Aga, + A^a, , (2.1) 

where cTj are Pauli matrices. Due to the disorder of the structure, a broad range of 
values for the TLSs parameters Aq and Ag is expected to appear for the ensemble of 
TLSs present in the amorphous solid. In order to be operative, the Standard Tun- 
neling Model has to specify the properties of this ensemble in terms of a probability 
distribution -P(Aq, Aq). The form of P(Aq, Aq) can be inferred from general consider- 



ations [AHV72,Phi72], as detailed in Appendices B.1.1 and B.1.2, and is furthermore 
supported by experiments [Esq98]. The result is 

P(AS,Ag) = ^, (2.2) 

with Pq ~ lO^'^J^^m^^, Aq > Amm (Amm fixed by the time needed to obtain 
a spectrum around the resonance frequency of the excited eigenmode), and e = 
a/ (Ag)2 + (Ag)2 < Emax, estimated to be of the order of 5 K. 

A last feature necessary to calculate the damping of acoustic waves caused by 
these low energy excitations is the coupling TLSs - phonons. As justified in Appendix 



B.1.3| the lattice distortions caused by the propagation of the wave affect mainly the 
TLSs asymmetry Ag, with the following final form for the hamiltonian describing the 
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coupled system TLSs + phonons: 



H = Y^ hujk^jbl^bkj + +[^o + Yl ^k,Mj + bk,j)] ff.} (2.3) 

fc.i Ag,Ag k,j 

Here b^j are the creation operators of the strain field waves u(r, t) whose components 
Ui{r,t) obey the wave equations 



= p- 



dt^ 



(2.4) 



The interaction hamiltonian is more specifically [Esq98] 



int 



ki,s 



ha- 



fe.-.s 



7 



2pL^h 



(bL + bk 



(2.5) 



where s is the polarization and i the component of the wavevector. L is the lateral 
size of the sample, while 7 is the coupling constant, with an approximate value 7 ~ 1 
eV. 



With the low energy TLS ensemble appearing in hamiltonian of eq.(2.3) character 



ized by a statistical distribution given by eq.(2.2), the strange thermal and acoustical 



properties of glasses at low temperatures mentioned in the beginning of this section 
can be satisfactorily explained [Phi87]: C ~ and k ~ are obtained, close to 
the experimental ones, whose discrepancies can be corrected by modifying slightly 
the simplified probability distribution. The heat release experiments are as well ex- 
plained, because there is a broad distribution of relaxation times for the initially 
thermally excited TLSs, and most of them will relax very slowly, transferring gradu- 
ally their energy to the phonons of the solid, increasing its temperature. Concerning 
the dependence of the phonon mean free path behavior on acoustic intensity, it is 
explained by the saturation of the population of excited TLSs, that eliminates the 
scattering mechanism of resonant absorption and decay (to be explained later in this 
chapter) at high intensities. The scattering at higher temperatures is dominated by 
non-resonant relaxational processes, also to be described later, associated to the equi- 
libration of the levels after a change in population, which occurs at a rate resulting 
in a mean free path for the phonon Im-es T~^. 



The range of applicability of the Standard Tunneling Model is determined by 
the approximations involved. The first one concerns the assumption that the TLSs 
interact only indirectly via their interaction with the phonons. This is false below a 
certain temperature T ~ 50 mK, when interactions of dipolar character among TLSs 
cannot be disregarded [Esq98, YL88, AMOS, ERK04]. The high temperature limit 
of applicability of the model is fixed by the moment when the TLS approximation, 
restricting the Hilbert space to the one spanned by the groundstates of the wells, see 



2. 3 Theoretical analysis of the damping of phonons due to TLSs 
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Figure 2.2: a) From [JPAH76]. Schematic view of an amorphous quartz structure, 
showing three possible types of two-state defects (A, B and C). b) Those two-states 
can be represented by a degree of freedom in a double well potential with an energy 
barrier of height V , width d, and asymmetry between wells Aq. At low temperatures 
only the ground states of both wells play a role in the dynamics, becoming mixed 
by the tunneling through the barrier (green arrow), of amplitude Ag. Thus the 
description is reduced to that of a Two-Level System. 



fig.( 2.2[ 3), is not valid anymore, at about T ~ 1 K. The model also assumes weak 
intensities of the acoustic waves excited, or otherwise anharmonic coupling among 
vibrational modes and higher order terms coupling TLSs and phonons would become 
non-negligible. 



2.3 Theoretical analysis of the damping of phonons 
due to TLSs 

Given an externally pumped vibrational eigenmode of frequency uoq, we will divide 
its dissipation due to the presence of TLSs into three parts: i) Resonant dissipation, 
ii) Friction caused by symmetric non-resonant TLSs, and iii) Relaxation processes 
due to the bias Aq. 

For usual excitation amplitudes the first one can be ignored. To obtain the energy 
dissipation rate corresponding to the other two mechanisms, we will proceed for each 
of them in two steps: first we will study what is the effect of a single TLS, expressing 
the results in terms of its parameters Aq, Aq, and in a second step we will sum the 
individual contributions to dissipation over all the TLS ensemble, characterized by 



the distribution (2.2) 
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2.3.1 Resonant dissipation 

Those TLSs with their unperturbed excitation energies close to loq will resonate 
with the mode, exchanging energy quanta, with a rate proportional to the mode's 
phonon population n^^;, to first order. For usual excitation amplitudes ~ 0.1 — lA the 
vibrational mode is so populated (as compared with the thermal population) that the 
resonant TLSs become saturated, and their contribution to the transverse (flexural) 
wave attenuation becomes negligible, proportional to riujy^ [AHSD74]. 

2.3.2 Dissipation of symmetric non-resonant TLSs 

Dissipation caused by a single TLS 

We will analyze here the dynamical behavior of one TLS coupled to the phonon 
ensemble: 

H = Aga, + A^a, + J] A^,,- . + bk,j) + J] ^kAjhj (2.6) 

k,j k,j 

This hamiltonian is nothing but the well-known spin-boson model, which applies to 
a great variety of physical phenomena related to open dissipative systems, and has 
been subject of intensive theoretical work [LCD+87, Wei99]. We will use some of 
this knowledge to describe the dynamics of the single TLS in terms of the spectral 
function 

A{u;) = ^|(0|a,|n)|'(5(a;-a;„ + a;o) (2.7) 

n 

where \n) is an excited state of the total system TLS plus vibrations. We will see 
how this function describes the absorption properties of a symmetric TLS coupled to 
a vibrational ensemble through the operator cr^. 

The analysis of the spin-boson model for a TLS with asymmetry Aq 7^ is much 
more cumbersome than the symmetric TLS case, if one tries to use the kind of tech- 
niques described in [LCD+87] that work for the symmetric case. Thus, to learn about 
the effect of the asymmetry Aq in the absorption spectrum of the TLS we will make 
a different kind of analysis, studying the delayed response of a TLS due to its bias 
Aq, and finally see its signatures in the absorption spectrum. A natural division of 
dissipative mechanisms will arrive as a consequence of this analysis: one contribution, 
the relaxation mechanism, will correspond to the presence of the bias, while the other 
will describe the dissipation occurring independently of the presence or absence of 
bias. 

The techniques we will use for the analysis of the symmetric case are all de- 
scribed in [LCD+87], mainly i) Simple renormalization-group arguments to obtain 
effective or "dressed" values A^jj, A^g of the TLS parameters from Ag, Aq due to 
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the interaction with high energy vibrational modes, ii) Characterization of the in- 
teraction with the different kinds of vibrational modes based on the criteria derived 



in the non-interacting blip approximation (NIBA, see Appendix B.2 for a brief re- 



view on the basic concepts and results), iii) Second-order perturbation theory in the 



interaction term of (2.6). These techniques are adequate for the regime where the 
Standard Tunneling Model holds, described just before the beginning of this section, 
and work better for very low vibration amplitudes of the externally excited vibra- 
tional mode. Even within the linear regime of the mode, above a certain amplitude 
other approaches designed to study strongly driven TLSs become more suited to 
describe the dynamics of at least part of the ensemble of TLSs, namely numerical 
methods as Quantum Monte Carlo [EW92, LWW98] or Numerical Renormalization 
Group [BTV03,BLTV05, ABV07, Wv07], or theoretical studies of the driven spin- 
boson model [GSSW93, GH98]. We will assume that the TLSs in resonance with the 
externally excited mode are saturated and thus do not contribute significantly to its 
dissipation, while for the rest of TLSs the vibration is assumed to be sufficiently weak 
for the techniques i), ii) and iii) to hold. 

Spin-boson model for an almost symmetric TLS interacting with the nanores- 
onator modes 

Now we take a TLS with Aq -C Aq. In principle it is coupled to three kinds 



of vibrational modes, those described in section L4 The analysis of the resulting 



spectral function (2.7) would be greatly simplified if we could somehow compare the 
effect of the different kinds of modes and conclude that two of them have a negligible 
influence as compared with the third. Fortunately the path integral NIBA approach 



developed by Leggett et al. [LCD+87] (see Appendix B.2) reveals that the degree of 
influence of a given bosonic bath huj^b^bk on the dynamics of a TLS coupled to it 

through ifint = cr^ J2k '^k {^k + temperatures and small enough coupling 

Afc, is completely determined by the spectral density of the bath 

/oo 
dte"''{H,^tit)H^ntm = Y,\>^k\^Siuj-uJk) (2.8) 
k 

This function corresponds to the spectral decomposition of the operator i/mt, telling 
us the rate at which transitions between the two eigenstates of the isolated TLS 
accompanied by excitation of a phonon happen, as a function of u. The three kinds 
of modes coupled to our TLS are characterized by a J{uj) ~ u'^. As explained in 



the end of Appendix B.2 this result, obtained for T = 0, changes at T > because 
the thermal energy increases the number of possible initial and flnal states, and J{uj) 
changes to J{lj,T ^ cj) ~ T*. Then, if temperature is lowered the modes with a 
J{uj) with the lowest power s (weakest T dependence) will affect the TLS dynamics 
more strongly than the others, which to a flrst approximation can be ignored. 
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Focussing on the case of a TLS coupled to the vibrations present in a rod, the 
starting hamiltonian is 

H = A^a^ + Ag(T, + (J,F{diUj) + i/eiastk (2.9) 

where diUj is a component of the deformation gradient matrix, F is an arbitrary func- 
tion, and inelastic represents the elastic energy stored due to the atomic displacements, 
which will correspond in second quantization to the last term in eq.(|2.6p. Changing 



basis to the energy eigenstates of the TLS, eq.(2.9) becomes H= eaz + [(Ag/£:)o"s + 
{/\ye)az]F{diUj) + inelastic- We are considering for the moment only slightly biased 
TLSs for whom Aq ^ Aq, so the last term can be ignored. A further expansion of F 
to lowest order in the displacement, together with a 7r/4 rotation of the eigenbasis, 
leads to: 

H = ea^ + 7(Ag/e)(T,9,Mj + //elastic , (2.10) 

where 7 is the coupling constant, assumed to be as in amorphous 3D solids 7 ~ 1 eV. 
The atomic displacements can be decomposed into the normal vibrational modes. 



which for the cuasi-lD case of nanorods are the ones obeying eqs.(1.2 1.4). The 
vibrational modes of a beam with fixed ends have a discrete spectrum, but we will 
approximate them by a continuous distribution. This approximation will hold as long 
as many vibrational modes become thermally populated, kT ^ hwfund, where Ufund 
is the frequency of the lowest mode. The condition is fulfilled in current experimental 



setups. The derivation of the second quantized version of the sort of eq.(2.6) for 



each of the modes and the corresponding spectral functions J{uj) are described in 



Appendix B.2.1 Here we state the main results: 

The compression and twisting modes lead to an ohmic spectral function for u <^ 
27Tc/R {R being a typical transversal dimension of the rod and c the sound velocity), 
when the rod is effectively ID. In terms of the Young modulus of the material, E, 
and the mass density, p, we get: Jcomp(t^) = ttc|c<j|, where, 

a, = ^\27r^ptw)-\E/p)--'/\ (2.11) 

A factor (Aq/£:)^ has been thrown away, as we consider almost symmetric TLS with 
Ag ~ e. We will proceed similarly with the other modes. The twisting modes 
are defined by the torsional rigidity, C = pt^w/3 (/i is a Lande coefficient), and 
I = J dSx^ = t^w/ 12 (where S is the cross-section). The corresponding spectral 
function is given by: Jtorsion(^^) = Cit\^\, where 

at = -i^C{M^ptwpI)-\pI/Cf'^ . (2.12) 

The analysis of the flexural (bending) modes differs substantially from the other 
ones, because they correspond to two fields that do not follow usual wave equations. 
The normal modes have a quadratic dispersion, ujj{k) = ^jEIj/{ptw) k^, leading to 
a sub-ohmic spectral function [LCD+87], 



J&eA^) = aby^y/uj , (2.13) 
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Figure 2.3: Sketch of the contributions to the spectral function which determines the 
dynamics of the TLS. 



with, 



7' (1 + z/)(1-2z/) ^pV/4 
a.v^ = 0.3^ ^(3_5,) bj ' (2-14) 



where z/ is Poisson's ratio and ojco — ^Elyl{ptw){2Txltf is the high energy cut- 
off of the bending modes (frequency for whom the corresponding wavelength is of 
the order of t, thus indicating the onset of 3D behavior). This example of sub- 
ohmic spin-boson model is the first result of this thesis to be stressed, because the 
problem of a TLS interacting with a sub-ohmic environment is interesting in its own 
right [LCD+87, KM96, Wei99, Sta03, VTB05, GW88, SG06, CT06, Khv04, IN92], and 
the systems studied here provide a physical realization. 

The bending modes prevail over the other, ohmic-like, modes as a dissipative 
channel at low energies. One may ask at what frequency do the torsional and com- 
pression modes begin to play a significant role, and a rough way to estimate it is to 
see at what frequency do the corresponding spectral functions have the same value, 
-''flex(t<^*) = -^comp,tors(^*)- The results are u;* ^ 'i^{\\vf{\-2vf{E j pfl^ l\t{Z-^vf\ 
for the case of compression modes and cj* ~ 300(1 — 2z/)^(£'/p)^/^/[t(3 — 5z/)^(l + z/)] 
for the torsional. Comparing these frequencies to the one of the onset of 3D behavior, 
Wco, they are similar, justifying a simplified model where only flexural modes are con- 
sidered. Collecting the previous results, we find the spectral function J{oS) plotted in 



Fig. 2.3 . 



TLS's dynamics 

The interaction between the bending modes and the TLSs affects both of them. 
When a single mode is externally excited, as is done in experiments, the coupling 
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Figure 2.4: a) Schematic representation of the irreversible flow of energy from the 
externally excited mode to the TLS ensemble, and from the ensemble to the rest of 
vibrational modes. This process can be viewed as a flow of energy from the excited 
mode to an ensemble of "dressed" TLSs, with their dynamics modified by the pres- 
ence of the vibrational modes, b) Left: Spectral function A{uo) of a single dressed 
TLS, weakly damped (F < Ag). A peak around u; = arises if the system is bi- 
ased, corresponding to the relaxational dissipation mechanism. Right: Total spectral 
function Atotioj) of the ensemble of dressed TLSs. 



to the TLSs will cause an irreversible energy flow, from this mode to the rest of the 
modes through the TLSs, as depicted in fig.( |2.4^ ). The dynamics of the TLSs in 
presence of the vibrational bath determines the efficiency of the energy flow and thus 
the quality factor of the excited mode. Taking a given TLS plus the phonons, its 
dynamics is characterized by the Fourier transform of the correlator {az{t)az{Q)), the 



spectral function A{uj), eq.(2.7). We will see that, using second-order perturbation 
theory, the dissipation caused by an ensemble of TLSs on a given externally excited 
vibrational mode of frequency c^o is proportional to Ylij ^jl'^o), with the index j run- 



ning over all the TLSs in the resonator, cf. eq.(2.18). We therefore study now A{uj 



High-energy oscillators with ^ adjust instantaneously to the value of cr^, 
"dressing" the wells' eigenfunctions and decreasing the overlap between them, re- 
ducing accordingly the tunneling amplitude Ag [LCD+87]. This adiabatic renormal- 
ization of Ag can either reduce Ag to a flnite value or suppress it completely. For 
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sub-ohmic coupling, eq.(2.14), the self-consistent equation for the renormalized value 

^TP.rt. iS- 



Aq exp{— afoY^u^ / dijj J [yj] I up'} 

^ren 



(2.15) 



This equation has no solutions other than A^en = if Aq <^ aluJco, so that the tunnel- 
ing amplitude of the low energy TLSs is strongly suppressed [SD85,LCD+87, Wei99, 
KM96]. The remaining TLSs experience a shift and a broadening of the spectral 
function function A[ij). For a typical Si nanoresonator with L ~ 1/xm, t ^ w ^ 100 
nm, ttf, ^ 1; and the shift is small, A^gn ~ ^q, and will be neglected in the following. 
In addition, A{uj) acquires low and high energy tails, whose derivation can be 



found in Appendix B.3 At zero temperature the low energy tail is: 



A{uj) oc ah 



(A; 



< Ag 



(2.16) 



for uj 3> Ag. 



while the high energy part takes the form A{uj) oc a6y^a^(Ag)^ci;~''/^, 
The main features of A{uj) are shown in fig. ( |2. 4) 3). Finally, we obtain the width of 
the resonant peak, r(AQ), using Fermi's golden rule, r(AQ) = 16ahy/u!co \/Ag. This 
description is valid for wavelengths such that 1/ L <^ k <^ l/max(t(;, t). For a biased 
TLS, the value of Aq is not renormalized by the phonons. 

As mentioned before, the coupled system TLSs -|- vibrations can be viewed, taking 
the coupling as a perturbation, from the point of view of the excited mode c^o as a 
set of TLSs with a modified absorption spectrum. The TLSs, dressed perturbatively 
by the modes, are entities capable of absorbing and emitting over a broad range 
of frequencies, transferring energy from the excited mode to other modes. The 
contribution to the value of the inverse of the quality factor, Q~^{(jJq), of all these 



non-resonant TLSs will be proportional to A"°^ 



^0 



E 



-r(a;o) 
=0 



A(Ag,Ag,^o) + 



EL™+r(cjo) ^(Aq) Aq, 1^0) ~ '^Poib\/^co/^Q , a quantity measuring the density of states 
which can be excited through i^int at frequency tuo, see Appendix B.3| for details. 
For an excited mode uq populated with phonons, Q~^{ujq) is given by 



AE 



(2.17) 



where Eq is the energy stored in the mode per unit volume, Eq ~ n^^^^hujQ/twL, and 
AE is the energy fluctuations per cycle and unit volume. AE can be obtained from 
Fermi's Golden Rule: 



AE 



tot 

off— res 
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X hw. 



2n 

X —n, 



A 



^off-res (^0) 



(2.18) 



and the inverse quality factor of the vibration follows. For finite temperatures the 

uJo,T) is done in Appendix B.3 It has to be noted that in 



calculation of ^off_res 
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experiments the observed linewidth is due to the total amount of fluctuations, so at a 
finite temperature it corresponds to the addition of emission and absorption processes. 
Thus in this context dissipation means fiuctuations, and not net loss of energy. The 
combined contribution of these processes is 



{Q 



^)off-rcs(^0 



, T) ~ lOP^t^/' 



w 



1/4 2 



coth 



hujQ 



(2.19) 



Before we continue with the third dissipative mechanism due to TLSs, a couple of 
words about the importance of multi-phonon processes is in order: Until now preva- 
lence of one-phonon processes in the interaction among TLSs and vibrational modes 
has been assumed, but at temperatures much higher than the frequencies of the 
relevant phonons, multi-phonon processes should be taken into account. A useful 
indicator to estimate if one-phonon processes suffice to describe the interaction be- 
tween the TLS and the bath is the Fermi's Golden Rule result for the linewidth: if 
the linewidth r(Ag,T) calculated in this way is much smaller than Ag, the approx- 
imation can be taken as good, while for r(Ag,T) > Ag multiphonon processes have 
to be considered. The expression for the linewidth at a finite temperature is 



T{e,T) = 16ab 



/i coth 



2kBTi 



(2.20) 



One can estimate, using the probability distribution (2.2), the total number of over- 



damped TLSs, T{e, T) > Aq, in the volume fraction of the resonator presenting amor- 
phous features, Vamorph- With Vamorph ~ Vfot/10, and using T{e,T) ^ 2rr(e,T = 
0)/e, the number of overdamped TLSs, e < T(e,T) ^ e < [SOaby/uJ^T]'^/^, is 
PQtwL [30aby/uJcoT] , which for typical resonator sizes L ~ l/im, t,w O.l/xm 
is less than one for T < 1 K. Therefore, unless the resonator is bigger and/or Pq too, 
the TLSs can be assumed to be underdamped. 



2.3.3 Contribution of biased TLSs to the linewidth: relax- 
ation absorption 

This very general friction mechanism arises due to the phase delay between stress 
and imposed strain rate. A detailed derivation for the case of TLSs in amorphous 
solids can be found in Appendix |B. 4 In our context, for a given TLS the populations 



of its levels take a finite time to readjust when a perturbation changes the energy 
difference between its eigenstates [J72,Esq98]. This time r corresponds to the inverse 



of the linewidth (2.20). For a coupling as the one described by Hint, a component k 
of the deformation gradient matrix, duk ~ {k,nk\{k/ y/LJk){al. + ak)\k,nk) associated 
to a vibrational mode \k,nk) modifies the asymmetry Aq between the energy levels 
£1^2 = =F^ a/ (^o)^ + (^0 + ikdukY {ik the corresponding coupling constant with the 
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Figure 2.5: Schematic representation of the levels' population evolution of an ensem- 
ble of 5 identical TLSs with a response time r ~ l/cj, where uj is the frequency of the 
bending mode excited, whose evolution is also depicted in the lower part of the figure. 
The delay r plus the bias give rise to the relaxational energy loss mechanism of 
the mode, see text. 




proper dimensions). So to first order the sensitivity of these energies to an applied 
strain is proportional to the bias 



d{dUk) 



and the response of the TLS will be oc (A*°*/^)^- derived in Appendix B.4, the 
imaginary part of the response, corresponding to Q^^ is also oc r/(l which 
in terms of A{ijj) is the lorentzian peak at = of fig. ( 2.4) 3). 



The mechanism is most effective when r ~ where uj is the frequency of 

the vibrational mode; then, along a cycle of vibration, the following happens (see 



fig. (2.5)): When the TLS is under no stress the populations, due to the delay r 
in their response, are still being adjusted as if the levels corresponded to a situation 
with maximum strain (and therefore of maximum energy separation between them, cf. 
£i,2(Mfc)), so that the lower level becomes overpopulated. As the strain is increased to 
its maximum value the populations are still adjusting as if the levels corresponded to a 
situation with minimum strain, thus overpopulating the upper energy level. Therefore 
in each cycle there is a net absorption of energy from the mechanical energy pumped 
into the vibrational mode. 

If on the other hand r ^ l/tu the TLSs levels' populations are frozen with respect 
to that fast perturbation, while in the opposite limit r <C l/cu the levels' populations 
follow adiabatically the variations of £12 and there is neither a net absorption of 
energy. 

For an ensemble of TLSs, the contribution to is [Esq98] 



Q-Uiu,T) = ^ I de / du^^^-^ — (2.22) 
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The derivation of eq.(B.43) only relies in the assumptions of an existence of well 



defined levels who need a finite time r to reach thermal equilibrium when a perturba- 
tion is applied, and the existence of bias | Ag| > 0. This implies that such a scheme is 
applicable also to our ID vibrations, but is valid only if the perturbation induced by 
the bath on the TLSs is weak, so that the energy levels are still well defined. Therefore 



we will limit the ensemble to underdamped TLSs for whom hT{e, T) < e. In eq.( B.43 
the factor cosh~^ {e/2kBT) imposes an effective cutoff e < ksT, so that in eq.(2.20 
one can approximate coth[e/2kBT] ~ 2kBT/e, resulting in T{e) ~ and thus 

the underdamped TLSs will satisfy e > [30ah-\/uJ^T]'^/^ . For T ^ [32ah^y^J^]'^ , which 



is fulfilled for typical sizes and temperatures (see Appendix B.5 for details): 



20Po7^ (1 + z/) (1 - 2i/) / p \ 1/4 VT 



Here V„ 



amorph 



t^/^W 

Vtot/W was assumed. 



E2(3-5i/) \EJ 



(2.23) 



2.3.4 Comparison between contributions to Q 
prevalence. 



Relaxation 



It is useful to compare the importance of the contributions to coming from 
the last two mechanisms. For that sake, we particularize the comparison to the case 
of the fundamental flexural mode, which is the one usually excited and studied, of a 
doubly clamped beam, with frequency ujq ~ 6.5{E / pY^'^t/L'^ (for a cantilever these 
considerations hold, with only a slight modification of the numerical prefactors; the 
conclusions are the same). The result is: 



Qre\ 



(T) 



off-res (^) 



fund 



300^1/2 (3_5;y) E 1 

~L (l + z/)(l-2z/)7^ 



(2.24) 



For a temperature T = IK, the result is as big as 10^ even for a favorable case 
t = Inm, L = l/xm, E = 50GPa, = 0.2, p = 3g/cm^. Then it can be concluded 
that for any reasonable temperature and dimensions the dissipation is dominated by 
the relaxation mechanism, so that the prediction of the limit that surfaces set on the 



quality factor of nanoresonators is, within this model for the surface, eq.(2.23): 



a"urfacc(^o,T)^Q-[(u;o,T) 



2^1/2 



(2.25) 



For typical values L ~ 1pm, t,w ~ 0.1pm, 7 ~ 5 eV, PoVamorph/Vtot ~ lO'^^J/m^, and 
T in the range lmK-0.5K the estimate for Q^rfacc ~ gives the observed order of 
magnitude in experiments like [ZGS+05, SCK+07], and also predicts correctly a sub- 
linear dependence, but with a slightly higher exponent, 1/2 versus the experimental 



fit 0.36 in [ZGS+05], or 0.32 in [SCK+07], see fig. (2.6). 
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T(K) 



Figure 2.6: Experimental results for the dissipation of the fundamental flexural (bend- 
ing) mode of three different doubly-clamped resonators. The dissipation is expressed 
in terms of Q~^, the inverse of the quality factor. (a) Results obtained in ref. [ZGS"'"05] 
for two silicon doubly clamped beams of dimensions 0.2 x 0.3 x 6/im and 0.2 x 0.3 x 7/im. 
Their fundamental bending modes, 14.586 and 12.028 MHz respectively, show a 
Q-i(T) ~ optimum fit for 0.1 < T < 1 K. (b) Resuhs of ref. [SCK+07] for 

a similar type of beam made of GaAs, of dimensions 0.5 x 0.5 x 15/im, exhibiting an 
power law dependence Q~^{T) ~ f^j. g gg < T < 3 K. These behaviors have to 
be compared with our model, who predicts Q^^(T) ~ T^/^, see text for details. 



2.4 Extensions to other devices 

2.4.1 Cantilevers, nanopillars and torsional oscillators 

The extrapolation from doubly clamped beams to cantilevers [MSL"^05] and nanopil- 
lars [SB04] is immediate, the only difference between them being the allowed (k,ti;(k)) 
values due to the different boundary conditions at the free end (and even this differ- 
ence disappears as one considers high frequency modes, where in both cases one has 
kn ~ (2n + l)7r/2L). All previous results apply, and one has just to take care in the 
expressions corresponding to the of the fundamental mode, for whom there is 
more difference between the frequencies of both cases, the cantilever one being cUq'^"* ~ 
{E / pY^H/L"^ as compared to the doubly clamped case, if^°-"^'P^'^ ~ Q.h^E / pY^'^t/ L"^ . 

2.4.2 Effect of the flexural modes on the dissipation of tor- 
sional modes 

The contribution from the TLSs + subohmic bending mode environment to the 
dissipation of a torsional mode of a given oscillator can be also estimated. We will 
study the easiest (and experimentally relevant [SRCR04]) case of a cantilever. For 
paddle and double paddle oscillators the geometry is more involved, modifying the 
moment of inertia and other quantities. When these changes are included, the analysis 
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follows the same steps we will show. 

Relaxation absorption. We assume, based on the previous considerations on 
the predominant influence of the flexural modes on the TLSs dynamics, as compared 
with the influence of the other modes, that the lifetime r = of the TLSs is 
The change in the derivation of the expression for comes in 

is different, which translates simply, in 



given by eq.(|2.20|) 
eq. 



eq. 



jtors 



2.21[) , where the coupling constant ^\ 

7^^Qj,g, and the corresponding prediction for 



B.43), into substituting 7 



20Po7?ors7' (1 + z/) (1 - 2Z/) / p X 1/4 



t^/^w 



UJQ 



(2.26) 



where now ujq is the frequency of the corresponding flexural mode. The range of 
temperatures and sizes for which this result applies is the same as in the case of an 
excited bending mode. 

Dissipation of symmetric non-resonant TLSs The modified excitation spec- 
trum of the TLS's ensemble, A^^_^^^{u!o,T), remains the same, and the change hap- 
pens in the matrix element of the transition probability of a mode |fco,no) appearing 
in eq.(2.18), {Xk^/ y/uJ^y. The operator yielding the coupling of the bath to the 



torsional mode which causes its attenuation is the interaction term of the hamil- 



tonian, which for twisting modes is given by eq.(B.30). In Appendix B.6 a de 



tailed derivation of eq.(2.27) is shown 

^1 



The main results are stated here. Again, 
Q~^{ujj) = AE/27iEq, where the energy Eq stored in a torsional mode (j)j{z,t) = 
y4sin[(2j — l)7rz/{2L)] sm{ujjt) per unit volume is Eq = A'^Ujp{t'^ + m;^)/48 {z is the 
coordinate along the main axis of the rod). Expressing the amplitude A in terms of 
phonon number Uj, the energy stored in mode \kj,nj) is 



Eo{kj,n) 



2 {t^w + wH)L 



+ w^){2nj + I) 



(2.27) 



the energy fiuctuations in a cycle of such a mode is 



27r 



27r '-f'^hujj 



(2.28) 



and the inverse quality factor 



T 



(2.29) 



For sizes and temperatures as the ones used for previous estimates the relaxation 
contribution dominates dissipation. 



2.5 Frequency shift 
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2.5 Frequency shift 

Once the quality factor is known the relative frequency shift can be obtained via 
a Kramers-Kronig relation (valid in the linear regime), because both are related to 
the imaginary and real part, respectively, of the acoustic susceptibility. First we will 
demonstrate this, and afterwards expressions for beam and cantilever will be derived 
and compared to experiments. 



2.5.1 Relation to the acoustic susceptibility 

In absence of dissipative mechanisms, the equation for the bending modes is given 
by —{12p/t^)d'^X/dt^ = Ed^X/ dz'^ . The generalization in presence of friction is 

l2pd'^X ,^ a^X 
-lf^-(E^x)^ (2.30) 

Where x is a complex-valued susceptibility. Inserting a solution of the form X{z, t) = 
j^^i{kx-ujt) ^ where k is now a complex number, one gets the dispersion relation oj = 
^t'^{E + x)/il2p)k^. Now, assuming that the relative shift and dissipation are small, 
implying Re(x) << E , Im(k) << Re(k), the following expressions for the frequency 
shift and inverse quality factor are obtained in terms of x'- 

g-i =Auj/uj = -lm{x)/E 

6uj/uj =Re(x)/2E ^ ' 

Therefore a Kramers-Kronig relation for the susceptibility can be used to obtain the 
relative frequency shift: 

^^(.,r,.-ipr..'«:;(^, ,2.32) 

UJ 2n uj' -UJ 

where P means here the principal value of the integral. 



2.5.2 Expressions for the frequency shift 



Relaxation processes of biased, underdamped TLSs dominate the perturbations of 
the ideal response of the resonator, as we have already shown for the inverse quality 
factor. For most of the frequency range, uj > [30aby/uJ^T]'^/^ , Q~^(u,T) ?a Ay/T/u , 
with A defined by eq.(2.23). The associated predicted contribution to the frequency 
shift, using eq.(2.32), is 




(2.33) 
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For low temperatures, u ^ [SOab^/uJ^T]'^/^ , the negative shift grows towards zero as 
5uj/lu(uj,T) ~ — v^T log[a;/T^/^], reaching at some point a maximum value, and de- 
creasing for high temperatures, cu < [SOab^yuJ^T]'^^^ , as 5uj/uj{uj,T) ~ Even 
though the prediction of a peak in 5u}/uj(T) qualitatively matches the few experimen- 
tal results currently available [ZGS+05, SCK+07], it does not fit them quantitatively. 

2.6 Applicability and further extensions of the model. 
Discussion 

As mentioned, the predictions obtained within this theoretical framework do 
match qualitatively experimental results in terms of observed orders of magnitude 
for Q~^{T), weak subhnear temperature dependence, and presence of a peak in the 
frequency shift temperature dependence. But quantitative fitting is still to be reached, 
while on the experimental side more experiments need to be done at low temperatures 
to confirm the, until now, few results [ZGS"''05, SCK"^07]. 

Applicability. The several simplifications involved in the model put certain con- 
straints, some of which are susceptible of improvement. We enumerate them first 
and discuss some of them afterwards: i) The probabihty distribution P(Aq, Aq), bor- 
rowed from amorphous bulk systems, may be different for the case of the resonator's 
surface, ii) The assumption of non-interacting TLSs, only coupled among them in an 
indirect way through their coupling to the vibrations, breaks down at low enough 
temperatures, where also the discreteness of the vibrational spectrum affects our pre- 
dictions iii) When temperatures rise above a certain value, high energy phonons with 
3D character dominate dissipation, the two-state description of the degrees of free- 
dom coupled to the vibrations is not a good approximation, and thermoelastic losses 
begin to play an important role, iv) For strong driving, anharmonic coupling among 
modes has to be considered, and some steps in the derivation of the different mech- 
anisms, which assumed small perturbations, must be modified. This will be the case 
of resonators driven to the nonlinear regime, where bistability and other phenomena 
take place. 

The solution to issue i) is intimately related to a better knowledge of the sur- 
face and the different physical processes taking place there. Recent studies try to 
shed some fight on this question [CRB07], and from their results a more realistic 
P(Ao, Aq) could be derived, which remains for future work. Before that point, it is 
easier to wonder about the consequences of a dominant kind of dissipative process 
which corresponded to a set of TLSs with a well defined value of Ag and a narrow 
distribution of Aq's of width Ai, as was suggested for single-crystal silicon [Phi88]. 
Following [Phi88], a Q~^{T) ~ x/T behavior is obtained for low temperatures T < Ai 
ff r(Ag, T) < uj, and a Q-^{T) ~ l/^/T ii r(Ag, T) > uj, while at high temperatures 
T > Ai a constant Q~^{T) ~ Qq is predicted for both cases. These predictions do not 



2.1 Dissipation in a metallic conductor 



35 



match better with experiments than the results obtained with P(Ag, Aq) ~ Pq/Aq, 
so issue i) remains open. 

We try now to have a first estimate of the temperature for which interactions 
between TLSs cannot be ignored. Following the ideas presented in [Esq98], we will 
estimate the temperature T* at which the dephasing time Tint due to interactions is 
equal to the lifetime t(T) = T~^{T) defined in eq.(2.20), for the TLSs that contribute 
most to dissipation, which are those with e ~ Ag ~ T. For them r~^(T) = T{T) ^ 
4006^/0^ a/T. The interactions between the TLSs are dipolar, described by ifmt = 
T.ijUi^'icrla^, with Ui^2 = &12M2, &12 verifying (612) ^ 0, (I612I) = Uq ^ 'y^/E 
[YL88, Esq98]. From the point of view of a given TLS the interaction affects its 
bias, Aj = (Aq)^ + UijCj, causing fluctuations of its phase Sej(t), which have an 
associated Tint defined by Sej{Ti^t)'Tmt ~ 1- These fluctuations are caused by those 
TLSs which, within the time Tint, have undergone a transition between their two 
eigenstates, affecting through the interaction ifjut the value of the bias of our TLS. 
At a temperature T, the most fluctuating TLSs are those such that e ~ Ag ~ T, and 
their density can be estimated, using P(Ag, Aq), as ut ~ P^kT. They will fluctuate 
with a characteristic time t(T) [A{)ai,^/uj^\/T]~^ , so for a time t < t{T) the amount 
of these TLSs that have made a transition is roughly n{t) ~ PQkTt/T{T). For a 
dipolar interaction like the one described above, the average energy shift is related 
to n{t) by [BH77] Se(t) ~ U^nit). Substituting it in the equation defining Tint, and 
imposing rint(T>K) = r(T*) gives the transition temperature T* ^ [Qah^/uJ^/ (f/oPo)]^- 
For example, for a resonator like the silicon ones studied in [ZGS"'"05], L = 6/im, 
t = 0.2/im, w = 0.3/im, the estimated onset of interactions is at T* lOmK. 

An upper limit Thigh of applicability of the model due to high energy 3D vibrational 
modes playing a significant role can be easily derived by imposing Thigh = ^^'^n- The 
frequency uj'^^ corresponds to phonons with wavelength comparable to the thickness 
of the sample, oj^^ = It:^ E j pjt. The condition is very weak, as the value for 
example for silicon resonators reads Thigh ~ 400/t, with t given in nm. At much 
lower temperatures the two-state description of the degrees of freedom coupled to the 
vibrations ceases to be realistic, with a high temperature cutoff in the case of the 
model applied to amorphous bulk systems of T ~ 5K. 



2.7 Dissipation in a metallic conductor 

Many of the current realizations of nanomechanical devices monitor the system 
by means of currents applied through metallic conductors attached to the oscillators. 
The vibrations of the device couple to the electrons in the metallic part. This coupling 
is useful in order to drive and measure the oscillations, but it can also be a source 
of dissipation. We will apply here the techniques described in [GJS04, Gui05] (see 
also [WP95]) in order to analyze the energy loss processes due to the excitations in 
the conductor. 
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Figure 2.7: Sketch of the distribution of charges in the device. When the system oscil- 
lates, these charges induce time dependent potentials which create electron-hole pairs 
in the metallic layer deposited on top of the beam, absorbing part of the mechanical 
energy of the flexural mode. See text for details. 



We assume that the leading perturbation acting on the electrons in the metal are 
offset charges randomly distributed throughout the device. A charge Q at position 
R creates a potential: 

(2.34) 

eo|R(t) - r(t)| 

at a position r inside the metal, see Fig. 12.7] . As the bulk of the device is an insulator. 



this potential is only screened by a finite dielectric constant, eo- The oscillations of 
the system at frequency ojq modulate the relative distance |R(t) — ?(i)|, leading to a 
time dependent potential acting on the electrons of the metal. 

The probability per unit time of absorbing a quantum of energy uoq can be written, 
using second order perturbation theory, as [GJS04, Gui05]: 



j dHrdtdt'V{i, t)V{r, t')lmx[i - f , t - t']e^'^«(*-*') 



(2.35) 



where Imx[r — r', t — t'] is the imaginary part of the response function of the metal. 

Charges in the oscillating part of the device. We write the relative distance 
as R(t) — r(t) = Rq — ?o + S'Rit) — 6rit) and expand the potential, whose time- 
dependent part is approximately 



V{r,t) 



Ro - fo 



6R{t) - 5f(t) 



For a flexural mode, we have that 



Rn 



(2.36) 



■ [5R(t) — 5r{t)] has turned into [Xq — 
Xo] ■ [6X{t) — 6x{t)] t- A - sin(co't), where t = ti^s + ^motai is the thickness of the beam 
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and A is the amplitude of vibration of the mode. Thus the average estimate for this 
case for the correction of the potential is 



SV{f,t) 



Q^t- A- sm{ujt) 



L3 



(2.37) 



L being the resonator's length. 

The integral over the region occupied by the metal in eq.(2.35) can be written as 
an integral over tq and Fq. The RPA dielectric constant of a dirty metal is: 



Imx(r — r', t — t') 



giq(f-f')gia;(t-t') 



\UJ\ 



e^Du\q\ 



(2.38) 



where e is the electronic charge, D = hv-pl is the diffusion constant, v-p is the Fermi 
velocity, / is the mean free path, and u ^ (fcFimetai)^/(^'i^F) is the one dimensional 
density of states. 

Combining eq.(2.35) and eq.(2.38), and assuming that the position of the charge, 
Ro is in a generic point inside the beam, and that the length sca les ar e such that 
^F^'^metai,^ "C L, we cau obtain the leading dependence of T in eq.(2.35) on L: 



r 



DuL \L 



t 



(2.39) 



where we also assume that \Q\ = e. The energy absorbed per cycle of oscillation and 
unit volume will be AE = {27i / ujQ)huJoTpi,/t'^ L = 27ihTpii/t'^L, and the inverse quality 
factor Qph(t^o) will correspond to 



QphK) 



1 AE hTpt 1 



2n Eo twL \puolA^ 



(2.40) 



where Eq is the elastic energy stored in the vibration and A the amplitude of vibration. 
Substituting the result for V one obtains 



Q 



-1 



2h 



IkpLH'^pUjQ \tmctal 



t 



(2.41) 



In a narrow metallic wire of width tmetai, w e exp ect that / ~ tmetai- 

Typical values for the parameters in eq.(2.41) are kp^ ^ lA, A ^ lA, I 

^ ^metal ~ 

lOnm ~ 10^ A, t ^ lOOnm ^ lO^Aand L ^ 1pm ^ lO'^A. Hence, each charge in the 
device gives a contribution to of order 10~^°. The effect of all charges is obtained 
by summing over all charges in the beam. If their density is uq, we obtain: 
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For reasonable values of the density of charges, hq — Iq^, Iq > lOnm, this contribution 
is negligible, Q^^ < 10^^^. 

Charges in the substrate surrounding the device. Many resonators, however, 
are suspended, at distances much smaller than L, over an insulating substrate, which 
can also contain unscreened charges. As the Coulomb potential induced by these 
charges is long range, the analysis described above can be applied to all charges within 
a distance of order L from the beam. Moreover, the motion of these charges is not 

correlated with the vibrations of the beam, so that now the value of 5R(i) — Srit) 

has to be replaced by: 

5R(t) - 5T{t) fti ^e^^°* , (2.43) 

and the value of |Ro ~ fo| ~ L. Assuming, as before, a density of charges uq — Iq, 
the effect of all charges in the substrate leads to: 

(0-1 ~ ~ n 



where the second result corresponds to the fundamental mode, ujq ^ Q.b{t/ L'^)\J E / p. 
For values L ^ l/xm, A ^ lA, kp^ ~ lA, I ~ tmctai ~ lOnm and Iq ~ lOnm, we 
obtain Q'^ ~ 10^^. Thus, given the values of Q^^ reached experimentally until now 
this mechanism can be disregarded, although it sets a limit to Q^^ at the lowest 
temperatures. It also has to be noted that this estimate neglects cancelation effects 
between charges of opposite signs. 



2.8 Conclusions 

Disorder and configurational rearrangements of atoms and adsorbed impurities at 
surfaces of nanoresonators dominate dissipation of their vibrational eigenmodes at 
low temperatures. We have given a theoretical framework to describe in a unified 
way these processes. Based on the good description of low temperature properties of 
disordered bulk insulators provided by the Standard Tunnehng Model [AHV72,Phi72, 
Esq98], and in particular of acoustic phonon attenuation in such systems, we adapt 
it to describe the damping of ID fiexural and torsional modes of NEMS associated to 
the amorphous-like nature of their surfaces. 

We have explained the damping of the modes by the presence of an ensemble of 
independent Two-Level Systems (TLSs) coupled to the local deformation gradient 
field diUj created by vibrations. Fiexural modes, with a high density of states at low 
energies, lead to sub-ohmic damping, which can modify significantly the distribution 
of TLSs. The problem of a TLS interacting with a sub-ohmic environment is interest- 
ing in its own right [LCD+87,KM96,Wei99,Sta03, VTB05,GW88,CT06,Khv04,IN92], 
and the systems studied here provide a physical realization. 
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The different dissipation channels to which this ensemble gives rise have been 
described, focussing the attention on the two most important: relaxation dynamics 
of biased TLSs and dissipation due to symmetric non-resonant TLSs. The first one 
is caused by the finite time it takes for the TLSs to readjust their equilibrium pop- 
ulations when their bias Aq is modified by local strains, with biased TLSs playing 
the main role, as this effect is oc [Ag/^(Ag)2 + (Ag)2]2. In terms of the excitation 
spectrum of the TLSs, it corresponds to a lorentzian peak around u; = 0. The second 
effect is due to the modified absorption spectrum of the TLSs caused by their cou- 
pling to all the vibrations, specially the flexural modes. A broad incoherent spectral 
strength is generated, enabhng the "dressed-by-the-modes" TLSs to absorb energy of 
an excited mode and deliver it to the rest of the modes when they decay. 

We have given analytical expressions for the contributions of these mechanisms 
to the linewidth of flexural or torsional modes in terms of the inverse quality fac- 
tor Q~^{u}q) = Aujq/ujq, showing the dependencies on the dimensions, temperature 
and other relevant parameters characterizing the device. We have compared the 
two mechanisms, concluding that relaxation dominates dissipation, with a predicted 
Q~^{uJo, T) ~ T^/^/cJo • Expressions have been provided for damping of flexural modes 
in cantilevers and doubly-clamped beams, as well as for damping of their torsional 
modes. 

Analytical predictions for associated frequency shifts have been also calculated. 
Some important successes have been achieved, like the qualitative agreement with a 
sublinear temperature dependence of Q~^{T), the presence of a peak in the frequency 
shift temperature dependence 6uj/uj{T), or the observed order of magnitude of Q~^{T) 
in the existing experiments studying flexural phonon attenuation at low temperatures 
[ZGS+05, SCK+07]. Nevertheless, the lack of full quantitative agreement has led to 
a discussion on the assumptions of the model, its links with the physical processes 
occurring at the surfaces of NEMS, its range of applicability and improvements to 
reach the desired quantitative fit. 

Finally, we have also considered the contributions to the dissipation due to the 
presence of metallic electrodes deposited on top of the resonators, which can couple 
to the electrostatic potential induced by random charges. We have shown that the 
coupling to charges within the vibrating parts does not contribute appreciably to 
the dissipation. Coupling to charges in the substrate, although more significant, still 
leads to small dissipation effects, imposing a limit at low temperatures Q^^ ~ 10^*^, 
very small compared to the values reached in current experiments [ZGS"*"05,SCK+07]. 



40 Chapter 2: Surface dissipation in semiconductor NEMS at low temperatures 



Part II 

Friction mechanisms in graphene 
and carbon nanotube-based 
resonators 
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Chapter 3 



Graphene-based structures: 
promising materials 



3.1 From 3D to OD, ID and 2D 

The versatility of carbon atoms to form all kinds of compounds with completely 
different properties can be hardly overemphasized [Pie93]. When carbon atoms com- 
bine to constitute a solid, they can form two types of bonds [MGM06, GMM06], see 



fig.(3.1): the very strong sp'^ bonds present in diamonds, where each atom has four 
first-neighbors located in the vertices of a tetrahedron whose center is occupied by 
that atom; and the more versatile sp^ bonds, with three identical coplanar sp bonds 
forming 120° angles between them, plus a fourth p-type orbital perpendicular to that 
plane. As a result an hexagonal honeycomb lattice with closest-neighbor distances 
of 1.4A tends to form, with in-plane a orbitals giving rigidity to the lattice and 
out-of-plane tt orbitals providing delocalized charge carriers. 

The stacking of these planes creates graphite, a highly anisotropic material thanks 
to the differences between the intra and inter plane bonds. The former are strong 
covalent bonds, whereas the latter are weak van der Waals forces, explaining the 
easiness with which the planes can slide one over the other, basis of the lubricant 
properties of graphite. Indeed, a band-structure calculation renders hopping ampli- 
tudes tab = 2.4 eV between nearest neighbors within a plane and tc = 0.3 eV between 
atoms of consecutive planes, separated by 3.5A. This anisotropy led Wallace already 
in 1947 [Wal47] to propose a simplified model where the interaction between planes 
is neglected, and a single graphite plane is studied to obtain approximately the elec- 
tronic properties of graphite, viewing it as a 3D solid formed by the stacking of 2D 
independent hexagonal crystals. 

Fortunately for the mesoscopics community this is not the end of the story of 
carbon compounds. The angles formed by the in-plane a orbitals can be modified 
to a certain extent, so that not only hexagons but also pentagons or heptagons of 
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C atoms can be formed. This automatically leads to non-planar structures, some of 
whom may be energetically stable. This is the case of some closed structures, carbon 
fullerenes [DDE96], whose main representative is the Cgo, a "football" ~1 nm in di- 
ameter composed of carbon hexagons and pentagons, called buckminsterfullerene, see 



fig. ( 3.2 i). The experimental discovery of this new class of 0-dimensional carbon al- 
lotropes in 1985 by a team leaded by Kroto, Curl and Smalley [KH0^85] was a break- 
through awarded with the Noble Prize in Chemistry in 1996, opening several research 
lines in mesoscopics still active. Studies have focussed on the physical and chemical 
properties of such nano-objects, like heat conductance and high-temperature super- 
conductivity [Had91], and applications have been proposed, for example in biomedics 
as vehicle for other molecules. They have been even used in double-slit experiments 
to test the validity of the wave-particle duality for such big molecules [ANVA+99]. 

Some years later a second breakthrough took place: Ijima's paper on the syn- 
thesis of carbon nanotubes [Iji91] caught the attention of the scientific community 
and boosted research on these carbon structures, an example of whom is depicted in 



fig.(3.2i). This figure shows a single- walled nanotube, the result of the folding of a 
single graphite plane, but there are as well multiple-walled nanotubes, with several 
concentric folded graphite planes. Soon their amazing properties, with much more 
technological potential than OD fullerenes, became apparent: depending on the way 
the carbon plane is cut and pasted when folding it to make the tube, it can be semicon- 
ductor or metal, and by doping them properly one can extend the scope to insulating 
and superconductor states; they are also good candidates for nanoelectronics, with ex- 
tremely low scattering rates due to their high degree of crystallinity, allowing electrons 
to travel micrometers ballistically at room temperature [DDAOl, FSOl]. Coming to 
their mechanical properties, they are already used in the fabrication of high strength 




Figure 3.1: The two main types of hybrid orbitals formed by C atoms to cre- 
ate a solid. Left: sp^ bonds. Right: sp^, in plane, bonds. Image from 
http : //en. wikipedia.org/wiki/OrbitaLhybridisation . 
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Figure 3.2: Eight examples of allotropes of carbon: a) Diamond, b) Graphite, c) 
Lonsdaleite, d) Ceo (BuckminsterfuUcrcne or buckyball), c) C540, f) C70, g) Amor- 
phous carbon, and h) single-walled carbon nanotube or buckytube. Image from 
http : //en. wikipedia.org/wiki/Allotropes_of_carbon . 



composites thanks to their ~1 TPa (!) Young's Modulus and capacity to relax to their 
initial form after severe deformations. Their sharp-ending geometry makes them ideal 
field-emitters, with application in flat panel displays, and good nanotools, serving as 
tips in several kinds of microscopes (Atomic Force, Scanning Tunneling, Magnetic 
Resonance Force and Scanning Nearfield Optical, Chemical/Biological Force Micro- 
scope tips). Moreover, from the fundamental research point of view they are also 
very attractive: the propagation of electrons only along the axis makes them effective 
1-dimensional systems where a Luttinger liquid picture substitutes the conventional 
Fermi liquid [EBF+01] apphcable to bulk metals [NP99], due to the lack of screening 
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of the Coulomb interaction, as compared to the 3D case. Quantum-mechanical in- 
terference effects, like the Aharonov-Bohm effect in presence of magnetic fields, have 
been also observed and studied [BSS+99]. 

The discovery of OD fuUerenes and ID carbon nanotubes was believed to put an 
end to the quest for carbon allotropes, as there were convincing arguments and strong 
experimental evidence against the thermodynamic stability of 2D crystals [Pei35, 
LanSO, Mer68, ETB06]. The surprise came in 2004, when a single plane of graphite, 
the so-called graphene, was identified over a Si02 substrate [NGM+04]. Its existence 
over such an amorphous substrate or even suspended in a different setup [MGK+07] 
can be reconciled with the previous arguments by noting that the carbon monolayer 
shows a certain amount of corrugation, with "mountains" ("valleys") typically ~ 1 
nm high (deep) and with a basis of ~ 30 nm diameter, fig. ( 3. 3) 3) . The associated 



total free energy is minimized, thanks to the strong reduction of thermal vibrations 
linked to the loss of planarity, which compensates the increase in elastic energy. 



As depicted in fig. ( 3.3 1), graphene can be considered the "mother" of 3D graphite, 
OD fuUerenes and ID carbon nanotubes, increasing the attractiveness of its study. We 
will focus in this second part of the thesis mainly in graphene-based devices, so now 
a brief overview on its properties and virtues follows. 



3.2 Graphene 

Graphene is a remarkable material from many points of view. The first one refers 
to its exceptional electronic quality: in spite of the primitive methods employed to 
produce it, graphene crystals are extremely pure, able to carry huge current densities 
~ 10^ A/cm^ (two orders of magnitude that of copper) [Wil06], with electron mobil- 
ities fi exceeding 15.000 cm^V~^s~^ under ambient conditions [NGM"'"04, NGM"'"05] 
and high carrier concentrations n > lO^^cm^ in both electrically and chemically doped 
devices. In terms of mean free paths ballistic transport is observed on the submicrom- 
eter scale (/mfp ~ 0.3;um at T = 300K). For a recent comprehensive review on the 
electronic properties of graphene see [CGP+07]. 

A second remarkable property: it is stable down to nanometer sizes, where the 
confinement of electrons leads to the opening of semiconductor gaps, changing com- 
pletely its electrical properties. Proceeding in this direction, Walt de Heer has pro- 
posed the use of electron-beam lithography to pattern a graphene sheet the proper 
way, carving in it nm ribbons serving as waveguides, and quantum dots to create, 
from a single starting sheet, a transistor. A semiconducting gap can be as well in- 
duced simply by growing epitaxially graphene on a SiC substrate, according to recent 
studies [ZGF+07]. Circuits using carbon nanotubes, however, need at some points to 
have them connected to the circuit through highly-resistive metal contacts, a major 
drawback also limiting the downscaling of silicon microchips, absent in graphene- 
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Figure 3.3: a) Image from [GN07]. The detailed understanding of graphene is im- 
portant for the comprehension of the other carbon materials deriving from it: OD 
buckyballs, ID nanotubes and 3D graphite, b) Pictorial representation of the ob- 
served crumpling of graphene monolayers giving them thermodynamical stability. 
Extracted from [MGK+07]. 



taylored circuits. These features raise hopes on graphene as a candidate for future 
electronics. 

But even if those hopes were not fulfilled there is another aspect which justifies 
by itself research on graphene: Using the band-model of Wallace [Wal47] one can see 
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that the electrons in a honeycomb lattice can be described at low energies in terms of 
quasiparticles obeying the (2+l)-dimensional Dirac equation. This means that at low 
energies the charge carriers in graphene behave as though they were massless Dirac 
fermions, but with a "speed of light" reduced with respect to vacuum photons by a 
factor 300. This is clearly reflected in fig. ( 3. 4) 3). 

The main steps for its derivation go as follows: In the hexagonal lattice there are 
two inequivalent atoms per unit cell, so the lattice can be seen as composed of two 



interpenetrated triangular sublattices A and B, see fig. (3. 4 1). An atom of a given 



sublattice is surrounded by three nearest neighbors belonging to the other sublat- 
tice. Once the filled a orbitals determine the honeycomb geometry, the properties of 
graphene can be calculated considering only the vr orbitals with energies close to ep- 
Starting with a tight-binding hamiltonian with nearest-neighbor interactions 



H = —t a\aj 



(3.1) 



<ij> 



where a\ is the creation operator for site i, so that a generic wave function can be 
written in this localized basis as 



(3.2) 



one can find the eigenfunctions and eigenenergies of H by solving the matrix equation 
for the coefficients ca and cb 



-t ^ . e*'^-"' 



= E{k) ( 
Cb J \ Cb 



(3.3) 



The vectors Uj and Vj are defined in fig. ( 3.4^ ). Solving the equation one obtains the 
dispersion relation represented in fig. ( 3.4[ 3~namely 



^(k) = ±tWl + 4cos2 



+ 4 cos 



akx 



cos 



2 '^^y 



(3.4) 



with a lattice constant a = 1.42A. With one electron per vr orbital (half-filling) the 
Fermi level is located just where the Fermi surface collapses into 6 points, two of 



whom inequivalent. Expanding eq.(3.4) for low energies around these so-called Dirac 



points the dispersion relation is linear in k 

E{k) = hvF\k\ 



(3.5) 



This is in stark contrast with the usual quadratic dispersion relation found in bulk 
metals or semiconductors. The fermi velocity vp is given hj hvp = (3/2)ta. Also 
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Figure 3.4: a) Graphene lattice can be seen as composed of two interpenetrated 
triangular sublattices A and B, with sublattice vectors Uj and Vj. b) Plot of graphene's 
bandstructure as derived with a simple band-model for the honeycomb lattice [Wal47] , 



eq.(3.4). When the graphene layer is neutral, the Fermi level is situated just at the 



6 so-called Fermi points, one of whom is zoomed at the right. At low energies, the 
dispersion relation around those Fermi points is linear, relativistic-like, E = hvF\\<-\ ■ 
Image from [Wil06]. 



within this low energy subspace around the Dirac points, one can Taylor-expand and 
derive a low-energy effective model: 

2 V - / f^VF{K<yx + kyay) , (3.6) 



50 



Chapter 3: Graphene-based structures: promising materials 




Figure 3.5: a) Schematic side view of the graphene-based resonator setup of ref. 
[BvdZV"'"07]. b) Scanning electron microscope image from [BvdZV"^07], showing a 
few-layer (A^ ~2) graphene resonator. Scale bar, 1 fim. 



arriving as anticipated at the Dirac hamiltonian in (2+1) dimensions, cfj being the 
Pauli matrices. Experiments confirm this relativistic behavior of pseudoparticles close 
to the Fermi point, most clearly seen in recent ARPES results [BOS^07], which 
reproduce fig. ( 3.4 d). 

The vanishing density of states at the Fermi points means that in neutral graphene 
there are few quasiparticles to screen the Coulomb interaction. The resulting strong 
Coulomb interaction results, according to some works [GGV99,DHT07], into a marginal 
Fermi liquid behavior, with a quasiparticle lifetime T{\e — eF\) ~ je — ^fI, as compared 
to the quadratic dependence present in Fermi liquids. Nevertheless it has to be noted 
that this question is not yet fully settled. When the graphene layer is doped with 
carriers the Fermi level separates from the Dirac points and locates in a region where 
the density of states is no more vanishing, so that now the Coulomb interaction can 
be efficiently screened, and a standard Fermi liquid picture applies. 

When one does not have a monolayer but a stack of them, the interlayer interac- 
tion modifies the dispersion relation, displacing slightly the conduction and valence 
bands creating the so-called electron and hole pockets, keeping a very low density of 
states (though no more vanishing) at neutrality. Numerous studies address the evo- 
lution of the physical properties as the number of layers of the system N is increased, 
starting from a single 2D graphene monolayer [NGM+04,NMM+06,Tho06,PP06]. It 
turns out that the 3D graphite electronic structure is very soon recovered, with about 
just 10 layers. Before that, the evolution of the electronic spectrum can be roughly 
classified into single layer, bilayer and 3 < iV < 10 cases, becoming increasingly com- 
plicated as grows. We will not move further into this issue, and limit ourselves 
as well to mention that many other interesting phenomena occur in graphene as a 
consequence of the interplay of its peculiar bandstructure and finite size effects, like 
the observation of an anomalous quantum Hall effect [NGM"^05,ZTSK05,NMM"'"06], 
a finite value of the electric conductivity at zero field [NGM+05, GN07], or unclear 
results about its magnetic response [MP06]. 



The extraordinary mechanical and electrical properties of graphene mono and 
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multilayers have not passed unnoticed to the NEMS community. Recently Bunch 
et al [BvdZV+07] demonstrated the first example of electromechanical nanobridge 
where the vibrating bridge is a single or multilayer graphene sheet, see fig. (3. 5). The 
advantages of the use of this material as compared with semiconductors like silicon 
or GaAs are clear: the strength of the a bonds provides stability even to a single 
suspended carbon plane, as opposed to crystal Si or GaAs devices, allowing for much 
lighter resonators (the lightest one can think of) with thus higher mass or force sen- 
sitivity. Moreover, those bonds result in an extremely high Young's Modulus, which 
together with the large surface area of these devices enhances the detector's sensitivity 
as well. As to the friction mechanisms limiting the quality factor of the resonator's 
modes, the high degree of crystallinity of graphene minimizes the dissipative pro- 
cesses coming from impurities and imperfections, which as we saw dominate damping 
in semiconducting nanoresonators at low temperatures. 

Graphene resonators of course present some drawbacks too, in the form of new dis- 
sipative mechanisms absent in semiconductor-based devices, or increased importance 
of some of the decoherence processes already described. The detailed description of 
these processes is the subject of the following chapter. As a final remark, we mention 
another interesting proposal for graphene as the basis for ultrasensitive gas sensors, 
based on the detection of local changes in the carrier concentration when a molecule 
is adsorbed [SGM+07]. This is only possible thanks to the extremely low intrinsic 
electronic noise in graphene. 
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Chapter 4 

Dissipative processes in graphene 
resonators 



4.1 Introduction 

In this chapter we will analyze several dissipative mechanisms affecting graphene 
resonators, both mono- and multilayer, like the ones fabricated by Bunch et al 
[BvdZV"'"07], comparing them with the semiconductor-based NEMS we have stud- 
ied in the first part of the thesis. We start with a brief description of these devices. 



4.1.1 Experimental setup 

The heart of the setup consists of a mono- or multilayer graphene sheet exfoliated 
on top of a previously patterned trench etched in a Si02 surface, figs. (3. 5) and (4.1 1). 
This Si02 layer, ~300 nm thick, is thermally grown from a starting underlying crys- 
talline Si wafer. The movement of the sheet is detected with a fast photodiode by 
looking at variations in the intensity of the reflected light from a 632.8 nm He-Ne laser 
focussed on the resonator, which are due to the size variation of the interferometer 
formed by the suspended graphene sheet and the silicon back plane, see fig.( 4.1[ 3). 
Bunch et al use two different driving schemes to induce motion: 



Electrical drive 



As shown in Fig.(4.1i), a gate voltage Vg is applied between the graphene and 
the silicon layers, with both a static (DC) component and a small time-varying (AC) 
one, inducing a charge q = CgVg, where Cg is the capacitance between the two layers. 
The attraction between the charges of the two layers causes an electrostatic force 
downward on the graphene. The total electrostatic force on the sheet is 

^-e. = Ic'X - ic-;n""(K,°" + = + c;ok, (4.1) 
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Figure 4.1: (a) Scheme of the resonator and the electrical actuation setup, (b) 
Optical actuation and detection scheme. From [BvdZV"'"07]. 



where C'g = dCg/dz is the derivative of the capacitance with respect to the distance 
between the plates of the capacitor. The second term of the electrostatic force leads 
to the motion of the oscillator, with increased amplitude if the driving frequency uj is 
close to a resonance ujq. 



Optical drive 

A 432 nm diode laser is intensity modulated at a frequency defined by the net- 
work analyzer and focussed on the suspended graphene, causing a periodic contrac- 
tion/expansion of the layer that leads to motion (see Fig. ( 4.1) 3)). 



4.1.2 Quality factors measured 

The quality factors Q measured at room temperature for the fundamental mode 
of several resonators, with thicknesses ranging from the monolayer case to ~ 100 
layers, are shown in fig.( |4.2^ ). The most prominent feature noted by the authors is a 
surprising lack of dependence of Q on the thickness of the resonators, leading them to 
suggest that the dominant dissipation mechanism is different from that of standard 
semiconductor NEMS. Comparing with results obtained for carbon nanotubes, they 
mention that the similarity of the structure and observed magnitude of Q perhaps 
points to a common dissipation mechanism, unknown until now for CNT resonators. 
They suggest also that an extrinsic mechanism such as clamping loss or fluctuating 
charge noise may dominate dissipation in graphene resonators. In a preliminary study 
of the temperature dependence of Q, they observe a sharp increase of the quahty as 
T is lowered, as shown in fig.( 4.2[ 3). 
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Figure 4.2: (a) The quality factor Q of the fundamental mode for all resonators stud- 
ied in [BvdZV+07] does not show a clear increase with thickness, contrary to semicon- 
ductor NEMS. Thus surface processes like the ones described in previous chapters do 
not seem to dominate damping in these devices, in agreement with graphene's well- 
known high degree of crystallinity. (b) Temperature dependence of Q for a resonator 
20 nm thick, showing a strong increase of the quality as the temperature is lowered. 
Courtesy from A. van der Zande (unpublished). 



4.1.3 Damping mechanisms present in graphene resonators 

Prior to analyzing in detail several of the mechanisms limiting the quality of the 
bending modes, we will briefly enumerate and describe the (known) most prominent 
ones, comparing them with the case of semiconductor NEMS: 

Novel or enhanced mechanisms with respect to semiconductor resonators 

The ability of graphene to sustain large carrier densities when an external voltage 
is apphed is a feature distinguishing it from the NEMS studied in the first part of 
this thesis, and as we have seen is used in one driving scheme. When the oscillating 
graphene layer is charged, the presence of charge in the Si electrode gives rise to 
friction associated to the excitation of electron-hole pairs in the graphene 
layer. This is due to the time-dependent Coulomb potential that the charges in the 
graphene layer feel, created by the charges in the electrode. 

And viccversa, when the Si electrode is charged, the time-dependent Coulomb po- 
tential due to the oscillating charged graphene layer causes excitation of electron- 
hole pairs in the electrode, absorbing mechanical energy. These processes are 
present in semiconductor devices only when the setup includes a metallic layer de- 
posited on top of the resonator (magnetomotive technique), and play usually a sec- 
ondary role due to the low charge density in the non-mobile part of the device. 
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Graphene resonators, on the contrary, have a higher charge density both in the 
graphene bridge and in the Si electrode in the setup with electrical drive, increasing 
the importance of this dissipative channel. Moreover, this mechanism is also present 
in the optical driving setup, as graphene is known to be, even in absence of applied 
voltages, locally charged, with so-called electron and hole puddles corresponding to 
the sheet corrugation [MAU"'"07], though its importance will be smaller as compared 
to the electrical drive setup case. 

Apart from their presence in the graphene bridge and the Si electrode, trapped 
chcirge in the Si-Si02 interface, as well as charge distributed throughout the volume 
of the oxide is known to exist [Sze81,NB82], limiting for example the performance of 
MOS-based devices. In our case it will add an extra Coulomb potential to the one 
generated by the charge in the Si electrode, enhancing the creation of e-h pairs in 
the graphene layer with the corresponding damping of the vibration. 

A process also missing in current semiconducting resonators which can damp the 
oscillation is the breaking and healing of bonds gluing the graphene sheet to the Si02 
substrate, conveniently denominated Velcro effect [VZO]. We will do some estimates 
based on the known chemical properties of Si02 surfaces and conclude that it may 
become a strong source of energy loss if the resonator is strongly driven into the 
non-linear regime, being otherwise negligible. 

Other mechanisms a priori non-negUgible in both kinds of resonators 

Some of the trapped charged impurities throughout the device are not static enti- 
ties, but have their own dynamics, coming frequently associated to the same defects 
causing surface dissipation in semiconductor NEMS. A clear example corresponds to 
charges switching between metastable trapping sites, which can be modeled by a set 
of two-level systems (TLSs). They generate the previously mentioned fluctuating 
chcirge noise, a real headache for the quantum computing community in need for 
long coherence qubit times [DH81,Wei88,PFFF02,BRS+04,SSMM07]. We will model 
the ones of the amorphous Si02 substrate with the usual TLS distribution employed 
in studies with amorphous materials [Phi88], concluding that this mechanism can be 
neglected. 

Attachment losses due to the coupling of the bridge's oscillation to the phonons 
of the rest of the device through the anchorage areas will as well be estimated. An- 
other source of dissipation which can be of importance, specially at high temperatures, 
is the thermoelastic coupling of the strain field of the excited resonator's bending 
mode to the rest of its modes, causing local variations in the temperature field. Its 
contribution will prove to be as well of secondary importance. 
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System properties 


Dimensions 




Thickness t 


10 ■ 10-9 m 


Width w 


10-*^ m 


Length L 


10-6 m 


Height above substrate d 


300 ■ 10-9 m 


Frequency /o 


100 MHz 


Amphtude A 


0.5 nm 


Carrier density pc 


lO^^cm-^ 


Properties of graphite 


Mass density 


2200 kg/m^ 


Elastic constants 




E 


10^2 Pa 


V 


0.16 


Debye temperature 


~ 570 is: 


Specific heat Cp 


700 J / Kg. K 


Thermal conductivity k 


390 W / m . K 



Table 4.1: Parameters used in the estimates presented in the chapter, adapted to the 
systems studied in [BvdZV"'"07]. Bulk data taken from [Pie93]. 



For all of the aforementioned mechanisms we will give a numerical estimate of 
their contribution to the quality factor at 300 K, and determine its temperature 
dependence, using a prototypical resonator characterized by the parameters given in 
table 14. 1[ The results are summarized in table 14.21 



Heating by the lasers 

In the optical drive setup, the intense actuation laser induces a local heating 
of the resonator, reducing the quality factor. To a lesser extent this mechanism is 
also present in the electrical drive setup due to the weak He-Ne monitoring laser. We 
will not analyze in detail this effect, concentrating in the electric drive setup. The 
interesting study of this heating process is yet a task to be performed. Nevertheless, 



a first observation can be made, based on fig.(4.2D): in that experiment the intensity 
of the lasers used is not varied as T is lowered, so its effect on Q should be roughly 
independent of T, but as a drastic increase of Q is observed, it leads us to think that 
the heating by the lasers is not the dominant source of dissipation, in that temperature 
range at least. 
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Q-i(T = 300i^') 


Temperature 
dependence 


Charges in the Si02 


10-^ - 10-6 


T 


Charges in graphene sheet 


10-2 


T 


and metaUic gate 






Velcro effect 


Absent 




Two-level systems 


10-22 


A + BT 


Attachment losses 


10-6 - 10-5 




Thermoelastic losses 


10-^ 


T 



Table 4.2: Contribution of the mechanisms considered in section l4l2l to the inverse 
quality factor Q~^(T) of the systems studied in [BvdZV"'"07]. 

4.2 Estimates for the damping caused by the dif- 
ferent mechanisms 

4.2.1 Coupling to fixed charges in the Si02 substrate 

We first calculate the dissipation induced in the graphene layer by a single fixed 
charge located in the Si02 substrate, and afterwards sum over all fixed charges. 

The time dependent component of the unscreened potential induced by a charge 
separated by a distance d in the vertical direction from the graphene layer, acting on 
an electron at position r in the graphene layer is given, approximately, by: 

(4.2) 

where A is the amplitude of the fiexural mode, and c^o its frequency. Its Fourier 
transform is: 

i>2tt poo iqr cos 8 . 

= 2ne^Ae-'^'^6{uj-uJo) (4.3) 
This potential is screened by the polarizability of the graphene layer [WSSG06]. In an 



RPA self-consistent calculation (see Appendix C.1.2 and Section 6.4.2) the screening 



is expressed in terms of a renormalization of the charge, replacing e"^ as follows: 

(4.4) 





q| 




l + (e2/|q|)Re[xo(|q 


,uj)] Re[xo(|q 





where xo is the susceptibility of the graphene layer. Rigorously this result is obtained 
only when 2D translational invariance is assumed, something we will do, ignoring 
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border effects. At low energy and momenta tfie value of the susceptibility tends to 
the compressibility (density of states) of the electrons in the layer. For the monolayer 
and multilayer cases we have: 



lim Re[xo(|q|,w)] 

|q|— >0,aJ— *0 



N = 1 
A^ ^ 1 



(4.5) 



where A^ is the number of layers and 7 is the interlayer hopping element. The result 
for the monolayer can be immediately obtained assuming a linear dispersion relation 
e{k) = hvpk, while for a stack with A^ layers, the model with one interlayer hopping 
element [GCP06] has been used, which gives rise to 2N low energy bands, most of 
which show a quadratic dispersion. 



Using Fermi's golden rule applied to the perturbation potential of eq.(4.2) screened 



as described in eq.(4.4) one arrives at eq.(C.6) of Appendix C.1.1 for the width of the 
graphene mode. Fourier transforming to momentum space (v(q, cu) = f(— q, tu)): 



rph ~ / rf2q|v(q)|2lmxo(q,o;o) 



(4.6) 



The low energy and momentum limit of graphene layers lmxo{^, ^^o) is given by 
[WSSG06,GCP06] 

Imxo(q,a^o)^<; (4.7) 



where, for A^ 7^ 1, p is the total carrier density. This last expressions are valid for 
lengths bigger than the mean free path, / ^ Imfp- 

The energy absorbed per cycle of oscillation and unit volume will be AE = 
(2TT/tuo)hiuoTph/twL = 27C hTpi^/twL, and the inverse quality factor Qp^i^^o) will cor- 
respond to 



1 AE hT 



ph 



2hT 



ph 



2n En 



twL IpulA^ 



(4. 



where Eq is the elastic energy stored in the vibration and M is the total mass of 



the resonator. Calculating eq.(4.6) and substituting in eq.(4.8) for a single graphene 



layer, a single charge turns out to give a contribution to the inverse quality factor of: 



Q 



-1 



1 



2h 



kpd Muj^d^ 



(4.9) 



where fcp = Tx^f)c-, and pc is the density of carriers in the graphene sheet. We can 
estimate this density as a function of the gate voltage Vg treating our system as a 
capacitor [NGM+04], 

eoeV^ 
e d 



Pc 



(4.10) 
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where eo is the permittivity of free space, e the relative permittivity of the material 
between the charged plates of the capacitor (in this case the graphene layer and the 
Si), and d the distance between the plates. In the region where the layer is suspended, 
e = 1, and for typical voltages applied Vg ~ 100 V the resulting carrier density turns 
out to be pc ~ 2 ■ 10^^ cm"^. Thus k^d ~ 10^ — 10^ > 1. Eq. (|4.9|) can be generalized 



to a graphene sheet with layers using eq.(4.7): 



2h 



W^d Mcuod"^ 



(4.11) 



The suppression with the number of layers is due to the increased screening in this 
system. 

The total contribution to the inverse quality factor is obtained by multiplying 
eqs.(4.9) or (4.11) by the total number of charges A^^ch- An upper bound to the ef- 



fective surface density of local charges, deduced from some models for the electric 
conductivity of graphene [NM07, HAS06], is pch ~ lO^^cm"^. In Appendix C.3 



some 



facts about these charges are collected from the literature, and an independent es- 
timate for Pch is made, obtaining a smaller value, eq.(C.28). The reason for this 



discrepancy is still unclear to us. Using the first estimate for pch and the parameters 
in Table 4.1 , we find Nch ~ 10^ and ~ 10~^^ at low temperatures. 

This mechanism leads to ohmic dissipation, as the energy is dissipated into electron- 
hole pairs in the metallic graphene layer. Hence, as deduced in Appendix |C.2.1 the 
temperature dependence of this mechanism is given by Q~^{T) ~ Q~^{0) x (kT/huo) , 
and Q-^ ~ 10"^ at 300 K. 

The analysis presented here does not consider the additional screening due to the 
presence of a metallic gate. In that case, one needs to add to the potential from a 



static charge, eq.(4.2), a contribution from the image charge induced by the gate. 



This effect will reduce the coupling between the graphene layer and charges in the 
vicinity of the gate. 



4.2.2 Ohmic losses at the graphene sheet and the metaUic 
gate. 

The electrons in the vibrating graphene layer induce a time dependent potential 
on the metallic gate which is sometimes part of the experimental setup. The energy 
is transferred to electron-hole pairs created at the gate or at the graphene layer. 
These processes contribute to the energy loss and decoherence of electrons in metallic 
conductors near gates [GJS04, Gui05]. 

The coupling between charge fluctuations in the two metallic systems is due to 
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long range electrostatic interactions. The corresponding hamiltonian is 



H 



Vscriz,V,t)p^\z,V,t) + / Vscr{0,r',t)p^(j') 



C 



G 



+ 



1 , lEt^w[rd^^^^ 



2pMtw 



2 12 



(4.12) 



where the indices G and C stand for the gate and graphene layer, respectively, pm 
is the mass density of the graphene sheet, and t, w, E its thickness, width and Young 
modulus, whereas (f){r, t) represents the vibrating amplitude field of bending modes 
and n = dL/d(j) is its conjugate momentum (L is the Lagrangian). The self-consistent 
screened potentials Vscr{z,r,t), Vscr(0,r,t) are calculated as a function of the bare 



potentials vo{z,r,t), fo(0,r,t) in Appendix C.4 



As in the case of eq.(4.2), the time-dependent part of the bare potentials couples 
the electronic degrees of freedom and the mechanical ones through the charge p(r) 
and amplitude of the vibrational mode, Aq, and would give rise to a term in the 
quantized hamiltonian of the form 



Hm oc p(f)Aq oc {bl + 6q) J][4+ii'^k + h-c] 

k,k' 



(4.13) 



where Aq and p(r) have been expressed in terms of creation and annihilation operators 
of phonons (q, Uq) and electrons of a 2D Fermi gas, respectively. 

But a realistic model requires taking into account the screening of the potential 
associated to these charge fluctuations. In terms of the screened potentials, the in- 



duced broadening of the mode (q, tUq) of the graphene layer is, extrapolating eq.(C.6) 
of Appendix C.1.1 to this [GJS04] 



a=G,C 



d^r'i ReV3"„(f , uj^) x ReV^^^.^r , tUq 



xlmx"[f-f,cjq]} (4.14) 



The static screening properties, limq^o R-ex"(q, 0), of the graphene layer and the gate 
are determined by their electronic compressibilities (densities of states), u'^ and 
respectively. We will assume that the distance between the graphene and the gate 
is much larger than the electronic elastic mean free path in either material, so that 
their polarizability is well approximated by the one of a diffusive dirty metal, see 



Appendices C.1.3 and C.1.4 



x"(q,^) 





a 


q 


2 




q 


2 


+ 


iu 



(4.15) 



where Z)" = fp/" is the diffusion constant, and is the elastic mean free path. The 
two dimensional conductivity is g'^ 
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We have assumed the gate to be quasi two dimensional. This approximation is 
justified when the distance between the gate and the graphene layer is much larger 
than the width of the gate. Appendix |C.3.2 shows how for typical oxide thicknesses 
and voltages applied it is fulfilled. In this situation, the broadening of the mode, 



eq.(4.14), can be expressed, analogously to eq.(4.6), as 



d k<^ \vscr{d,\i,uj^)\ Imx + |vs 



,(0,k,cuq)|2lmx^ 



(4.16) 



The screened potentials for a graphene layer oscillating in an eigenmode (q, tUq) 
of amplitude Aq, have in a first approximation only one momentum component, 
fscr(k, cUq) = fscr (q, '^q)5(k — q), and these components are (see Appendix C.4) 



fscr(c?, q,t^q) 

t;scr(0,q,u;q) 



|q| 



(4.17) 



2uCuG l-e-29d 



where (? = |q| and pc is the charge density in the graphene layer. The results for 



r(co'q) and Q 



can be formulated in terms of the total charge in the graphene 



layer, Qc = f d'^rpc ^ L x w x pc- 

In the limit of short separation between the layers, d L, which is the situation 



present in current experimental setups, substitution of eq.(4.17) into eq.(4.16) leads 
to 



(4.18) 



The limit -D|qp ^ uj for the imaginary part of the susceptibility of a dirty metal, 
Imx(q, u;) ^ cjzz/Dlqp, has been used, and is justified in Appendix C.1.3 The first 



term in the summation describes losses at the graphene sheet, and the second at the 



gate. The associated inverse quality factor, according to eq.(4.8), is given by 



hQ 



+ 



G.2 1 



(4.19) 



To make numerical estimates, we use the parameters in Table 4.1 , with u (E) 



E/2'Kh'^Vp, vf ~ 10^ m/s for a single layer of graphene, and i7\E) = {N'y)/h?Vp 
for a stack of layers [GCP06]. Carriers in graphene stacks have large mobilities 
[NGM"''04], and we take D'~^v^ k, 10^. Typical charge densities for the graphene 



layer discussed in the previous subsection, pc ~ 10 cm , leading to a total 

charge Qc ~ 10^. For these parameters, the contribution of the graphene sheet 
is ~ 10"*^. The relative contribution from the gate depends on the distance 
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to the graphene sheet. For a Si layer with D'^u'^ ^ 10^ and at short distances, the 
contribution to the damping from the gate is of the same order as that of the graphene 
sheet. 

Damping is associated to the creation of e-h pairs in a metal, which implies that 
this mechanism is ohmic. The inverse quality factor should thus increase linearly with 

-1 -10-2 at 300 K. 



temperature as described in Appendix C.2.1 leading to Q 



4.2.3 Breaking and healing of surface bonds: Velcro effect. 

In the fabrication process of the device, the graphene flake is deposited on the 
silica (Si02) substrate. The surface properties of silica are largely determined by the 
nature of silanol (SiOH) groups present there [MM90], due to the hydroxyl functional 
group that has chemical reactivity, as compared to the inert siloxane (Si-O-Si) surface 



[SG95]. In Appendix C.3.3 more information relevant to us is given on silica and its 
surface. 

Graphene becomes linked to the silica surface through hydrogen bonds created 
by those silanol groups. When the flake is set into motion, some of this bonds may 
repeatedly break and heal (the Velcro effect [VZO]), causing dissipation of the energy 
stored in the vibration. Numerical estimates are difficult to make, but nevertheless 
two qualitative arguments showing that its role in the damping is probably negligible 
can be presented: 

i) This mechanism is expected to be temperature independent, in contrast with the 



strong decrease of friction observed as temperature is lowered [BvdZV"''07], fig.(4.2D). 

ii) The elastic energy stored in a typical graphene oscillator of lateral dimensions 
w ~ l/xm is Eq = {1 / 2) pvo'^tu'lA^ ~ lOeV, when the amplitude is A ~ Inm. This 
means about ~ lO^^eV per nm^. On the other hand, the energy per hydrogen bond 
is about 10-^eV, and typical radical densities at Si02 surfaces are [DPX98] ~ Inm-^. 
Hence the elastic energy available on average for each hydrogen bond is much less than 
the energy stored in the bond. Only rare fluctuations, where a significant amount of 
energy is concentrated in a small area will be able to break bonds, and in this way to 
induce energy dissipation. Note, however, that this argument ceases to be valid for 
very large amplitudes > 30nm. For higher amplitudes, this mechanism can induce 
significant losses. 



4.2.4 Fluctuating charge noise: Dissipation due to two-level 
systems. 

As we have seen in Part I of this thesis, in an amorphous material an atom or a few 
atoms can have two nearly degenerate configurations. A vibration modifies the energy 
difference between these situations, and this coupling leads to the damping of acous- 
tic phonons in amorphous Si02 [Esq98]. We have seen this mechanism is expected 
to dominate friction in many semiconductor NEMS [MHE"'"02]. But with graphene 



64 



Chapter 4-' Dissipative processes in graphene resonators 



resonators, we expect the graphene sheet to show a high degree of crystalhnity, so we 
will only consider two-level systems (TLSs) in the rest of the structure. 

The TLSs can only dissipate energy if they are coupled in some way to the vi- 
brating graphene sheet. A possible mechanism is the existence of charge impurities 
associated to these defects (fluctuating charges), which are electrostatically coupled 
to the conducting electrons in the graphene. 

We expect this mechanism to be less effective in the device considered here than 
in NEMs made of semiconducting materials, as now the TLSs reside in the Si02 
substrate, not in the vibrating structure. The coupling, arising from long range 
forces, will be in comparison accordingly suppressed, by a factor of order (a/d)", 
where a is a length comparable to the interatomic separation, and n describes the 
decay of the coupling {n = 1 for the Coulomb potential between charged systems). 

The temperature dependence of the contribution of TLSs to is determined by 
the density of states of the modes coupled to the TLSs and the distribution of TLSs 
in terms of their paramenters (tunneling amplitude Ag and bias Aq) [SGC07a].The 
hamiltonian describing the coupling of the effective TLS's and the oscillating graphene 
sheet is given by [SGC07a] 

H = ea, + 7^cr, J] Ak(6k + + I] ^^k&L&k (4.20) 

k k 

where e = v^(AgpT(A§)2, 7 is the coupling constant, which will be strongly sup- 
pressed in these devices as compared to attenuation of acoustic waves in amorphous 
materials, 7 ~ leV x (a/d)"^, h\. represent the phonon creation operators associated to 
the different vibrational modes of a sheet, and Ak(6k + &!k) represents the coupling 
to the strain tensor Uik- There are two types, compression modes (longitudinal waves) 
and bending modes. The damping is due to the initial transfer of energy from the 
vibrational mode studied by the experimentalists to the TLSs, which in a second step 
transfer this energy to the rest of the modes. The properties of the spin-boson model, 
eq.(4.20), are fully determined by the power-law s of the spectral function [LCD+87], 
J{uj) = \^\k^Q/ 5{ijj — Uk) ~ aul'^^uj''^ , where ujk is the frequency of mode k, 
a is an adimensional constant and is the upper cutoff of the phonon bath. For 
this system, compression modes gives rise to a superohmic, s = 2, bath, while the 
bending modes constitute an ohmic bath, s = 1, and thus will prevail as a source of 
dissipation at low temperatures [LCD+87]. We will therefore restrict our analysis to 
the dissipation caused by the ohmic component of the vibrational spectrum. 
A detailed calculation of the dissipation due to this mechanism is carried out in 
Appendix C.5 Here we state the main results. One arrives at J{uj) = au, with 

a^4,^YPl^l±g!!(l^ ,4.21) 



ht^E^/^{9 + J 



3u_ 
-2u' 



Here v is the Poisson ratio of graphene. Choosing fairly symmetrical TLSs, Ag/e ~ 1, 
for the parameters in table 4.1, a ~ 10~^ x [a/dy^, very small. The relaxation of the 
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TLSs dissipates the energy of the vibration, giving rise to the following expression for 
the inverse quality factor [Esq98] , 



where u = Ar/e, e-max ~ 5 K, and (mVI — u'^)^^ comes from the probability density 
of TLS's in an amorphous solid, like Si02. Q~^{lj,T) is a function of C(a;,T), the 
Fourier transform of the correlation function C{t,T) = {az{t)az{0))T- For biased 
TLSs and ck <S 1 an extensive analysis of C{u,T) is performed in [Wei99], where 
several expressions are provided in different limits. Using them, the estimate for 
g-i(cu,T) follows: 



In the range of temperatures of current experiments (5K<T<300K), the depen- 
dence of dissipation with T is weak, and ~ 10~^ x (a/d)^ ~ 10~^^. The main 
uncertainty of the calculation has been the use of the TLSs' distribution assumed 
for amorphous solids [Phi87], but due to the small value of a a weak dissipation is 
expected also with a modified distribution. Thus the conclusion is that the relative 
importance of TLSs damping is much smaller for graphene than for other NEMs 
devices [CR02, Ble04, ER05] . 

4.2.5 Other friction mechanisms 

Attachment losses. 

The main expressions needed are given in [PJ04]. When d ^ t, and d is much 
smaller than the wavelength of the radiated elastic waves in the Si02 substrate, the 
contribution to the inverse quality factor is given by 



where the superscript O applies to the sihcon oxide, and stands for Poisson's ratio. 
The range of values of the quality factor varies from ^ 5 ■ 10~^ for a graphene 
monolayer, to Q^^ ~ 5-10^^ for a stack with 30 layers and t = lOnm. These quantities 
probably overestimate the attachment losses, as they do not include the impedance 
at the Si02-graphene interface. 

This damping process due to energy irradiated away from the resonator should 
not depend on temperature. 




(4.22) 




(4.23) 




(4.24) 
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Thermoelastic effects. 

When the phonon mean free path of the acoustic phonons is shorter than the 
wavelength of the mode under study, the acoustic phonons can be considered a dis- 
sipative environment coupled to the mode by anharmonic terms in the ionic poten- 
tial [Zen38, LROO, U06]. These anharmonic effects are described by the expansion 
coefficient, a, and the thermal conductivity, k. We follow the analysis in [Zen48]. For 
a rectangular beam vibrating at a frequency u the inverse quality factor is 

where E is the Young Modulus, Cp is the specific heat at a constant pressure, and 
Tz is the thermal relaxation time associated with the mode, which in the case of 
a flexural vibration is given by tz = t^Cp/ {tt'^k). This estimate assumes that the 
graphene sheet is weakly deformed, and that the typical relaxation time is associated 
to the diffusion of phonons over distances comparable to the thickness of the sheet. 



Although better approximations are available in the literature [LROO], eq.(4.25) is 



enough for an estimate of the order of magnitude of Q ^ . Using the parameters from 



table 4.1 , for t = 10 nm and / ~ 100 MHz, we find that ujtz <^ 1, and 



g^i(r = 300K)^ 5-10-^ (4.26) 



4.3 Extension to nanotube oscillators. 

The analysis presented here can be extended, in a straightforward way, to systems 
where the oscillating part is a nanotube. 

We expect in these devices a larger impedance between the modes of the nanotube 
and those of the substrate, so that attachment losses will be suppressed with respect 
to the estimate presented here for graphene. 

The damping mechanisms which require long range forces between the moving 
charges in the nanotube and degrees of freedom of the substrate (fluctuating and 
static charges) will not be significantly changed. A nanotube of length L at distance 
d from the substrate will interact with a substrate area of order (L + rf) x rf. A similar 
estimate for a graphene sheet of length L and width w gives an area ~ {L+d) x {w+d). 
As L ^ w ^ d ^ 1/im, the two comparable. 

On the other hand, ohmic losses induced in the nanotube will be reduced with 
respect to the two dimensional graphene sheet, as the number of carriers is lower in 
the nanotube. 

Finally, we expect a longer phonon mean free path in the nanotube, which implies 
that thermoelastic effects will be reduced . 
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4.4 Conclusions. 

We have considered six possible dissipation mechanisms which may lead to damp- 
ing in a graphene mesoscopic oscillator. The main results are summarized in Ta- 



ble 4.2 . We expect that the calculations give the correct order of magnitude and 
dependence on external parameters. 

We find that at high temperatures the leading damping mechanism is the ohmic 
losses in the metallic gate and the graphene sheet. This effect depends quadratically 
with the total charge at the graphene sheet, which can be controlled by the gate 
voltage. 

When speaking about the room temperature prevailing damping source, however, 
a disturbing fact obtained in ref. ( [BvdZV"'"07]) has to be noted: the independence of 
the magnitude of Q at room temperatures for both the optical and electrical driving 
setups. Unless a mechanism like optical heating were of the order of the predicted 
dissipation due to ohmic losses, we are missing something. In principle optical heating 
does not seem to be an important friction source, as it should lead to an almost 
temperature independent Q{T), contradicting results like the one shown in fig.( 4.2[ 3). 



At low temperatures attachment losses limit the quality of the vibration. If the 
resonator is strongly driven, a new damping mechanism may come into play, the 
Velcro effect, which may limit substantially the quality factor as compared with the 
slightly driven case. The high crystallinity of the resonators eliminates the main 
source of dissipation in semiconducting resonators, namely surface-related effective 
TLSs coupled to the local strain field. 

These conclusions apply with only slight modifications to carbon nanotube-based 
resonators. 

A deeper knowledge on the charge distribution throughout the device as a func- 
tion of the parameters characterizing the fabrication process and the operation of 
the resonator for the different setups would prove of great interest, as it may lead 
to substantial improvements in the quality factor through the design of charge min- 
imization strategies, and shed some light on the prevailing damping mechanisms for 
these devices. More measurements of the quality factor as a function of temperature 
would be also desirable for these purposes. 
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Part III 

Dissipation of collective excitations 
in metallic nanoparticles 
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Chapter 5 

Metallic clusters: optical response 
and collective excitations 



5.1 From atoms to bulk 



^ Nanometer-sized systems are the natural playground where the crossover from 
atomic-like to bulk behavior is observed. Among the simplest and most studied ones 
are clusters of alkaline or noble metal atoms, with radii a between 0.5 and several nm 
(few atoms to hundreds of thousands, respectively, see Fig. 5.1). The simpler elec- 
tronic structure of alkaline atoms renders the theoretical modeling of alkaline clusters 
easier as compared to the case of noble metals (with both s and d electrons), but their 
high reactivity imposes stringent conditions on their production: clusters cannot be 
formed embedded in a solid matrix, so they are created, starting from an alkaline 
vapor at a high temperature, by expansion into vacuum within an inert carrier gas 
flux, a method called seeded beam expansion [KKS82]. Noble metals circumvent this 
problem, allowing for a more controlled fabrication and study by ion implantation 
in an inert matrix and subsequent high temperature annealing [DYN+QQ, OKNSOO]. 
Moreover, their chemical stability is required for their use as biological markers, one 
of their most promising applications [BTM"'"02, CTB+03]. We will comment further 
on this issue later. 



The crossover size to observe bulk behavior depends on the physical property un- 
der study, and is moreover never clearcut. Two ingredients are the main responsibles 
of several new and interesting phenomena occurring in this interregnum. On the one 
hand, the spectrum of the electronic states remains discrete, with the energy levels 
broadened by the interactions with the environment. On the other, interferences and 
coherent behavior of electrons happen due to their coherent motion inside the system: 

^For a very good and detailed review on the issues covered in this chapter see [Wei06]. It has 
been the main bibhographic source of this chapter and the beginning of the next. 
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inelastic or phase-breaking scattering lengths are larger than the size of the cluster. 

The strongest evidence of the quantization of the electronic spectrum in metallic 
nanoparticles is the so-called electronic shell structure observed for the first time in 
1984 by Knight et al. [KCdH+84]. In the nuclear shell-model [BB75, RS80] a given 
nucleon of a nucleus can be considered to a first approximation as confined by a mean- 
field potential created by the rest of nucleons. Analogously for clusters, the valence 
electrons originating metallic bands in the bulk can be thought of as dwelling in a 
mean-field potential created by the rest of valence electrons and all the ionic cores 
comprising the cluster. An immediate consequence of such a model is the existence of 
certain "magic numbers" of valence electrons (or correspondingly, atoms) for whom 
the system is energetically more stable, namely those filling angular-momentum shells, 
as in the nuclear shell-model. In [KCdH+84], the histogram of clusters as a function of 
their size detected in a molecular beam of sodium seeded in argon showed a clear pref- 
erence for A/" = 8, 20, 40, 58 and 92, which was explained with the use of such a shell 
model with a spherically symmetric potential V{r) of the Woods-Saxon form [WS54]. 
It turns out that a sphere is a reasonable approximation for electronically closed-shell 
nanoparticles [dH93, Bra93], while for open-shell structures distortions occur due to 
the Jahn- Teller effect [dH93, JT37] . We will focus our attention in this thesis on the 
first, simpler case. 

Bulk-like behavior is frequently observed for surprisingly small clusters. Further- 
more, the size dependence of quantitative features can already be smooth for very 
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small systems. For instance, the binding energy of Nag clusters is quite close to 
that of Nag, even if these two systems are very different from the molecular point of 
view [dH93]. Such a continuity points, together with the aforementioned evidence of 
electronic shells [KCdH+84, KCdHS85], towards the relatively minor importance of 
the ionic cores and supports the descriptions based on the jellium model, where the 
conduction electrons are subject to a uniform neutralizing background [Bra93]. 

Electronic properties of metallic nanoclusters owe their uniqueness to two facts: 
first the already mentioned size range, where an evolution from discrete to continuum 
spectrum takes place, and second, the existence of a surface, a finite boundary that, 
due to the small size of these systems, has profound implications, the most relevant of 
which is the appearance of the so-called surface plasmon, a collective excitation that 
dominates the optical response for all clusters except the smallest ones, for whom 
transitions between single particle levels dictate the optical response. Indeed, the 
question of how large the size of the nanoparticle has to be in order to observe this 
collective excitation in such a finite system has the surprising answer that a very 
small cluster, like for instance Nag, may already be enough [SKM+91, YB91]. 

As we will see in next subsection, the surface plasmon corresponds to the ex- 
citation of a periodic displacement of the electronic center of mass with respect to 
the positive ionic background. The jellium approximation for this background allows 
to decompose the electronic Hamiltonian into a part associated with the electronic 
center of mass, a part describing the relative coordinates (treated in the mean-field 
approximation), and finally a coupling between the two subsystems. The coupling 
between the center of mass and the relative coordinates causes decoherence and dis- 
sipation of this collective state once it is excited. 

The treatment of the surface plasmon as a quantum particle (the electronic center 
of mass degree of freedom) therefore provides a model system for the study of deco- 
herence and quantum dissipation in confined nanoscopic systems, where the role of 
the electronic correlations is preponderant and the excitation spectrum arises from 
a finite number of particles. The finite number of electrons makes the term "en- 
vironment" not completely justified in this situation, and leads us to the following 
question, which will be addressed in this thesis: How large in size do we need to go to 
be allowed to describe the relative coordinates of the electron gas as an environment 
damping the collective excitation? The existence of several experimental techniques 
to observe the surface plasmon dynamics, and the use of such excitation in biomedical 
applications add even more interest to the subject. 

Based on the mentioned decomposition and mean-field approach an extensive and 
powerful analytical study of the plasmon dynamics is feasible, incorporating the ef- 
fects of finite size, a dielectric material, or the finite temperature of the electron 
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gas. With the use of semiclassical expansions it allows to i) obtain the paramet- 
ric dependence (in size, temperature, Fermi energy, etc) of the addressed quanti- 
ties like the surface plasmon lifetime or its resonance frequency, ii) compute the 
optical properties of arbitrarily large clusters, and with some extensions also de- 
scribe smaller ones, iii) be extended to nonlinear processes such as those resulting 
from the interaction with strong laser fields iv) address coherence and dissipation of 
the surface plasmon in both weak and strong laser field driving of the nanoparti- 
cle [YB91, YB92, MW J02, MWJ03, WMWJ05, WIJW06, WIWJ07, WWIJ07] . 

Such a wealth of analytical results and gain of knowledge demands the use of 
certain approximations, some of which have been assumed until now to be valid 
without proof. In the context of this thesis we will first present the theoretical scheme 
and then justify its approximations, confirming the validity of previous results and 
putting the approach on a more firm basis. In particular, we will determine which are 
the energies of the electronic excitations that are active in the damping of the surface 
plasmon, and what is the characteristic response time of the large enough electronic 
environment. The latter is important in justifying the Markovian approximation that 
is assumed to describe the dynamics of the surface plasmon coupled to particle-hole 
excitations [WIWJ07]. Other complimentary techniques, like the time-dependent 
local density approximation (TDLDA), will not be analyzed. We will start by giving 
a first physical image of the surface plasmon, and summarizing available experimental 
techniques for its study and use in biomedical applications. 



5.2 Surface plasmon excitation 

The study of the response of metallic nanoparticles to externally applied electro- 
magnetic fields reveals the main features of the electronic spectrum just outlined, 
with the optical response dominated by a collective resonance, the surface plasmon. 
What happens to the system when such a collective state is excited? For optical 
fields, the wavelength is bigger than the diameter 2a of the system and the frequency 
is such that the light electronic cloud is able to react to the variations of the field, 
whereas the ionic background, three orders of magnitude more massive, cannot follow 
the perturbation and remains approximately static. The result is the uniform relative 
displacement of electron and ionic clouds, creating unbalanced charges at the surface 
of the nanoparticle (a dipole) which try to restore the equilibrium situation in absence 
of an external field. The dipole is forced to oscillate with the frequency of the applied 



field (depicted in fig. 5.2), and its magnitude peaks at the so-called surface plasmon 
frequency, the natural frequency of the dipole, fixed by the electronic density, the 
spherical geometry, and the charge and mass of the electrons. So the surface plasmon 
is nothing but the self-sustained state of motion for which the center of mass of the 
negative charges oscillates with respect to the center of mass of the more massive 
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2a « K 



E(t) = E cos((ot) 




Figure 5.2: Schematic representation of the surface plasmon resonance in small metal- 
lic clusters. The wavelength A of the exciting electrical field E(t) is much larger than 
the diameter 2a of the cluster. Thus, the electrical field throughout the nanoparticle 
is nearly uniform. The shaded regions symbolize the ionic background, while the 
dashed circles represent the electronic cloud. The "-" and "+" signs account for the 
uncompensated negative and positive charges. The electrical field exerts a force on 
the positively charged ions and an equal force on the electrons in the opposite direc- 
tion, separating their centers of mass mainly due to the motion of the light electrons. 
The ions, being much more massive, remain almost at rest. The resulting surface 
charges create a dipolar restoring field (thick arrow) opposing E(t). 



ionic cores. 



Applying classical Maxwell's equations to the case of a dielectric metal sphere 
with appropiate boundary conditions, Mie obtained in 1908 [KV95, Mie08] a first 
estimate for this frequency, the classical Mie frequency lum = ^p/ V^, where Up = 
[A'Knf.e'^ /rrieY/'^ is the bulk plasma frequency and e, me, and rie denote the elec- 
tronic charge, mass, and density, respectively. In a real cluster surface effects lead 
to a small reduction of the surface plasmon frequency with respect to the Mie value 
ujy[ [Bra93,GGI02,WIJW06]. Typical values for the surface plasmon energy are about 
several eV, belonging thus to the energy range corresponding to bound single-particle 
states [BB94]. 



76 



Chapter 5: Metallic clusters: optical response and collective excitations 



5.2.1 Probing and using the 

Several techniques show and use the 
electronic properties of metallic clusters. 



surface plasmon 

existence of the surface plasmon to study 
We will mention three of them: 



Photoabsorption experiments 

The cluster is illuminated by a light of a certain frequency, and the amount of 
absorbed photons is determined and given in terms of the photoabsorption cross 
section (t{uj), which is the probability that the system absorbs a photon times the 
area illuminated by the photon beam. The presence of the single-particle states 
derived in the electronic shell model should be signalled by several peaks in a{uj), but 
they are usually overwhelmed by the peak of the collective excitation [BCLS93], thus 
indicating a substantial transfer of spectral weight to build the surface plasmon. 



Noble metal nanoclusters as biological markers 

An ideal optical label for large molecules should generate an intense optical signal; 
it should also be small, durable, chemically inert, and apt to bind to the molecule of 
interest in a controlled manner. Noble metal clusters meet all of these requirements 
when illuminated with a frequency close to the one of their surface plasmons. They 
overcome the main shortages present in other candidates like fluorescent dyes or 
nanocrystals of II- VI semiconductors, as they do not photobleach like the former 
and do neither optically saturate at reasonable exciting intensities like the latter. 
Several techniques take advantage of the excitation of the surface plasmon as a means 
to track the nanoparticle adhered to the molecule of interest: nanoparticles bigger 
than 40 nm in diameter can be imaged, based on Rayleigh scattering, in an optical 
microscope by illuminating in dark-field at the plasmon frequency [SSMSOO], with 
differential interference contrast (DIG) and video enhancement [JG88], or with total 
internal reflection [SGH+00]. For smaller nanoparticles, whose Rayleigh scattering is 
too small to be detected, local heating associated to its strong light absorption causes 
a temperature change which can be monitored by a sensitive interference method 



similar to DIG [BTM+02], see fig. 5.3 



Pump-probe experiments: plasmon dynamics 

The use of femtosecond pulsed lasers in time-resolved pump-probe spectroscopy 
has rendered possible to experimentally address the surface plasmon dynamics [BHMDOO, 
FVF"'"00, LKLA99, LUGSOl]. An initial ultrashort laser pulse lasting around a hun- 
dred femtoseconds pumps energy into the nanoparticle by exciting the surface plas- 
mon mode. After a certain time delay a second, much weaker, pulse, the probe laser, 
tests the absorption spectrum of the hot electron gas in the cluster. By repeating 
the experiment for several values of the delay an image of the time evolution of the 
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Figure 5.3: Noble metal nanoparticles can be used as biological tracers, thanks to 
their strong light scattering and absorption in the optical range due to the existence 
of the surface plasmon collective excitation. From ref. [BTM"'"02]. (A) Differential 
interference contrast (DIG) based on Rayleigh scattering of light and (B and C) pho- 
tothermal images measuring local heating due to strong light absorption, of a sample 
containing 300-nm-diameter latex spheres, 80-nm-diameter gold spheres, and 10-nm- 
diameter gold spheres. The heating intensity was 30 kW/cm^ (B) or 1.5 MW/cm^ 
(C). The 80-nm gold spheres appear on all three images and allowed direct compar- 
isons of the same imaged sample area. In the DIG image, the 10-nm particles are 
not visible, whereas the 300-nm latex spheres give very strong signals and the 80-nm 
gold spheres give weaker ones. In (B), a photothermal image at low power shows 
only the large metal particles. In (C), at high heating power, the 10-nm particles are 
clearly visible, whereas the 80-nm gold spheres saturate the detection capabihty. The 
strongly scattering latex spheres are completely absent in the photothermal images, 
as they do not absorb light as efficiently as a metallic cluster. 



relaxation process is obtained. 

Usually the experiments are performed on an ensemble of noble metal clusters 
embedded in an inert matrix, taking care that typical separations between clusters 
guarantee independence of their response. Several timescales are involved in the re- 
laxation, linked to the different characteristic times of the degrees of freedom partic- 
ipating: surface plasmon, single-particle electron excitations, nanoparticle phonons. 



and embedding matrix phonons. As schematically depicted in fig. 5.4, the energy is 
initially stored in the surface plasmon collective mode, causing a non-thermal, highly- 
correlated, electronic distribution /(e). In a femtosecond scale the surface plasmon 
decays through electron-electron and electron-surface scattering, transferring its en- 
ergy to the rest of (single-particle) electronic degrees of freedom, which accordingly 
heat and reach a thermal equilibrium distribution with temperatures up to hundreds 
of degrees, depending on the intensity of the pump laser. Phonons of the nanoparti- 
cle, with typical frequencies ~ 10^^ Hz, come into play on a much longer timescale. 
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Figure 5.4: Prom ref. [Wei06]. Sketch of the relaxation processes in a metaUic 
nanoparticle after an excitation by a strong and ultrashort femtosecond laser pump 
pulse, /(e) is a sketch of the electronic distribution, while fi is the chemical potential. 
The arrow in the uppermost image symbolizes the effect of the pump laser on the 
electronic system. 



of the order of a picosecond, and absorb part of the energy stored in the electronic 
distribution, cooling it. The relaxation process ends with the transfer of energy to 
the phonons of the matrix where the nanoparticle is embedded, on the timescale of 
several picoseconds, reaching the system as a whole thermal equilibrium. 

In experiments, the normalized difference of transmission of the sample with and 
without the pump laser field, called differential transmission AT/T, is the quantity 
measured to study the relaxation process: AT/T = {%n — %ff)/%ff. An example is 
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Figure 5.5: Experimental differential transmission AT /T as a function of the delay 
time between the pump and probe pulses, for different probe photon energies. The 
nanoparticles in this experiment are made of copper and embedded in a glass matrix. 
They have an average diameter of 10 nm. The plasmon resonance is located close to 
hwsp ~ 2.2 eV. (Reproduced from Ref. [BMCD95]) 



given in fig. 5.5, showing the differential transmission as a function of the time delay 
between the pump and probe pulses, for different probe energies in the vicinity of the 
surface plasmon resonance. The decrease of AT /T as a function of time reflects the 
relaxation of the electronic energy to the lattice. 
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Chapter 6 



Electronic dynamics in metallic 
nanoparticles: theoretical model 

This chapter describes a very successful theoretical framework for the analytical 
description of the electron dynamics in metallic clusters, expressed in terms of two 
types of excited states, the single-particle e-h excitations, and the many-body surface 
plasmon collective state. In a basis of relative -|- center of mass electronic coordinates, 
e-h and surface plasmon excitations will correspond to the relative and center of mass 
coordinates, respectively. Some key assumptions needed for the conceptual soundness 
of the model will be underlined, whose justification will be provided in the final part 
of the chapter. 

6.1 Electronic hamiltonian: relative and collective 
coordinates 

The systems we will analyze are neutral, spherical, metallic nanoparticles of ra- 
dius a containing N valence electrons. They thus have, as mentioned in the previous 
chapter, closed angular momentum shells {N — 8,20,40,58..., the so-called "magic 
numbers" [Bra93]). We will work within the adiabatic Born-Oppenheimer approx- 
imation [AM76], treating the ions as static, compared with the lighter and faster 
valence electrons. Moreover, we will work within the jellium approximation, which 
neglects effects of the ionic structure and treats the ions as a continuous, homoge- 
neous positively charged background of total charge +Ne, producing a charge density 
in the case of a spherical nanoparticle Pi{r) = {Ne/V)Q{a — r), with V = Ana^/d the 
volume of the cluster and Q{x) the Heaviside step function. The subsequent model 
is simple enough to be applied for the description of big clusters containing up to 
several thousands of atoms, and the effects linked to the ionic structure turn out to 
be important only for the smallest cases, < 10 [Bra93]. 
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The electronic hamiltonian for an alkaline-metal nanoparticle in vacuum is in this 
context given by 



^ r ^2 

1=1 *J=1 



where e and mg are the electronic charge and mass, Fj is the position of the i*'' con- 
duction electron, = |rj|, Pj is the associated momentum, and U{r) represents the 
electrostatic interaction between the electrons and the ionic background. 

For a noble metal cluster embedded in a dielectric medium, the electron-electron 



interaction in (6.1) is modified by the presence of the d-electrons and the matrix in 
terms of two dielectric constants and e^, and the expression for the Coulomb po- 
tential is much more cumbersome due to the loss of translational invariance. For a 
thorough study of this case see [Wei06, WMWJ05]. 

The shape of the confining potential U{r) for a charged sphere of radius a and 
total charge +Ne, harmonic inside the nanoparticle with the Mie frequency um = 
N jmeO? and Coulomb-like outside, can be readily obtained using Gauss's theo- 
rem, 

U(r) = — -(r' - 3a^)Q(a - r) -Qir - a) (6.2) 

2a'^ r 



If the many-body eigenstates |/) and eigenenergies Ef of the hamiltonian (6.1) are 
known, the photoabsorption cross section <7{uj) can be calculated using Fermi's golden 
rule [Mer70] 

A 2 

a{u) = ^ 5^ K/|;.|0)p5(na; - + E,) (6.3) 
/ 

The ground state is |0), and the electric field is assumed to be along the z axis. From 
o{uS) the frequency oj^i and linewidth of the surface plasmon are immediately 
obtained. But the calculation of the many-body eigenstates becomes exceedingly 
difficult for clusters bigger than a few electrons, and further approximation schemes 
have to be developed. 

The most useful one starts by changing the electronic coordinate basis to the 
center of mass R = ^ • rj/A^ and relative coordinates r'- = rj — R. The corresponding 



conjugated momenta are P = XliPi Pi = Pi ~ P/^? respectively, and (6.1) 
becomes 

H = ^ + Hrei + X^[f/(|r: + R|) - U{r[)] (6.4) 

6 ^ 1 
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where the hamiltonian for the relative-coordinate system reads 

1 



N 



H, 



rel 



i=l 



r n' 2 



+ 



E 



2 ^ \r' 



(6.5) 



A first glance comparing eqs. (6.1 ) and (6.5) may make us wonder where is the differ 



ence between them, and what is our motivation to go from (6.1) to (6.4). Regarding 



the difference between hamiltonians, note that one of the relative coordinates is a 
function of the rest, as Yli^'i ~ 0- '^^^ change of variables, on the other hand, is 
justified as it is the right one if we are interested in the study of the surface plas- 
mon collective excitation, which was shown in the previous chapter to correspond to 
electronic center-of-mass oscillations. 

Assuming small displacements for the center of mass, |R| <C a, the dependence 



on the relative coordinates in the last term of (6.4) drops out, f/(|r' + R|) — f/(r') ~ 
R- V[/(r') + (1/2)(R ■ V)^f/(r'), where derivatives are taken at r' = r. Choosing the 
oscillation axis of the center of mass in the ^-direction, R 



ZSz, one obtains with 



(6.2) 



R • Vf/(r') 
(R- V)'[/(r') 



z'e{a - r') + ^e(r' - a) 



\e{a-r') + - 



3cos2 9' 



r'3 



e(r' 



(6.6) 
e(a-r') 



where in the last line a second-order coupling in Z between the center-of-mass and 
the relative-coordinate system for r' > a has been neglected as compared with the 
first-order one expressed in the second term of the right-hand side of the first line. 



Inserting (6.6) into (6.4) 



H 



INe 



N 



2Nm, 2 a 



(6.7) 



i=l 



where the coupling between the center-of-mass and the relative coordinates is, to 
first order in the displacement R of the center of mass 



N 



1=1 



R=0 



(6.8) 



R=0 



given in (6.6). Notice that due to the spread of the electrons' 



with R ■ V?7(r^) 

wavefunctions outside the geometrical boundary of the cluster, r = a, there is a finite 
probability for them to have > a, which due to the sum over 0(a — r^) in (6.7) 



results in a reduction of the center-of-mass oscillation frequency from the classical 
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Mie value, ujm = t^Af a/1 — Nout/N < ujm- This quantum correction is the so-called 
spill-out effect [dH93,Bra93,KV95]. 

The structure of (6.7) is typical for quantum dissipative systems [Wei99]: The sys- 
tem under study {Hem) is coupled via He to an environment or "heat bath" described 
by Hrei- The system-environment coupling results in dissipation and decoherence of 
the collective excitation. 



6.2 Mean-field approximation and second quanti- 
zation of the hamiltonian 

The presence in H^ei of electronic interactions makes the problem still a formidable 
one and demands further simplifications. A mean-field approximation for the relative- 
coordinates' degrees of freedom, together with a second quantization procedure, will 
render a tractable hamiltonian suited for the study of collective excitations in metallic 
clusters. 



The harmonic oscillator hamiltonian Hem for the center-of-mass is quantized as 
usual with bosonic creation and annihilation operators [BF04], 




NmeCjM y J i 

NrrieUJM y i 

2h ^/2NmJSr^i 



Pz 
Pz 



(6.9) 



where Pz is the conjugated momentum to Z. Eigenstates \n) of if, 
y/n\n — 1) and b^n) 



+ l\n 

subtracting the zero-point energy /icDm/S, as 



satisfy b\n) = 

1), and Hem is expressed in terms of b and b\ 



He 



huj]\jb^b 



(6.10) 



In the mean-field approximation for the hamiltonian Hrei, eq.(6.5), we forget the 



fact that one of the relative coordinates is not independent of the rest (hopefully 
causing 1/A^ errors which become negligible as the number of electrons N grows). 



and substitute Hrei by the mean-field approximation of eq.(6.1), mimicking the effect 



of the rest of electrons on a given one with an effective one-body potential ^(r) felt 
by the latter: 



N 



N 



Hrei-^J2HscKp[] = J2\^ + Vir'^ 



2mp 



(6.11) 



We have used V{r) = V{r), valid for a spherical nanocluster. Numerical calculations 
using time-dependent local density approximation (TDLDA) within the spherical jel- 
lium model at zero temperature provide information about the shape of the potential 
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r/ao 

Figure 6.1: From [WMWJ05]. Self-consistent mean-field potential \^ as a function of 
the radial coordinate (in units of the Bohr radius = 0.53A) from TDLDA numerical 
calculations (at T = 0) for a nanoparticle with = 832 valence electrons. The 
mean distance between electrons = (3/47me)^/^ = 3.02ao (silver), corresponding 
to a ~ 28.5ao; for Na = 3.93ao [AM76], yielding similar results for V^(r), as the 
e = 1 case is equivalent to the one of an alkaline-metal cluster in vacuum. For most 
purposes in analytical calculations V can be approximated by a step function with 
discontinuity at r = a. 



V{r) [Eka84a,Eka84b,GK85,GDP96,CRSU00]. In this approach the original hamil- 



tonian (6.1) is replaced by a self-consistent hamiltonian like (6.11) where V{r) is a 
local function containing exchange and correlation terms. Results show that for most 
purposes V{r) can be very well approximated in analytical calculations by a step 
function, with the jump located at r = a. An example of TDLDA result for V{r) can 



be found in fig. 6.1 



From the self-consistent single-particle hamiltonian (6.11) a second-quantized ver 



sion immediately follows in terms of fermionic creation and annihilation operators 
and c„, 



rel 



4, 



(6.12) 



where cj^ creates the one-body state \a) whose wavefunction in the coordinate repre- 
sentation satisfies Hgci^ai^) = eaipai''^). 



The coupling term He between the center-of-mass and the relative coordinates 
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becomes in terms of the bosonic and fermionic operators and cj^ 



H^ = A{b^ + b)J2daf3clc(3 

afi 



where A = ^JhmeUj\i/2N , and 



(6.13) 



za 



dap = {a\ zQ{a — r) H —Q{r — a) 



(6.14) 



is the matrix element between two eigenstates of the unperturbed mean-field problem. 
In the constant A the approximation ojm ~ ujm has been included, as it leads to 1/a^ 
negligible corrections in the quantities one wants to compute, like the lifetime of the 
surface plasmon, using the second quantized approximation for (6.1), 



H = mjMb'^b + A(6t + h)Y^ dapcicp + ^ e^c^c^ 

Q /3 a 



(6.15) 



6.3 Surface plasmon linewidth and finite size ef- 
fects 

6.3.1 Mechanisms causing the decay of the surface plasmon 

The finite lifetime of the surface plasmon excitation can be related to several 
processes. When Mie performed his study of the optical response of metal clusters 
using Maxwell's equations, he used a simple free-electron Drude model [AM76] for 
the dielectric constant describing the metal. Within this approximation, 

e{u;) = e,-—^— (6.16) 

where Up = \/ A.Txnf.e^ jmf. is the plasma frequency, 7^"^ the relaxation or collision time, 
while e, mg and n^, stand for the electron charge, mass and bulk density, respectively, 
is the dielectric constant of the metal without the free-carrier contribution. Insert- 



ing (6.16) in the expression he had derived for the photoabsorption cross section, he 



obtained 

""^"^^ 2c (e, + 2e^)2(^-c^A^)2 + (7./2)2' ^ 
where is the dielectric constant of the medium in which the particle is embedded. 



and V the volume of the cluster. Expression (6.17) predicts the existence of the surface 
plasmon at the Mie frequency cum, with a linewidth 7^. This value of 7^ accounts for 
the various processes leading to elastic or inelastic scattering in hulk metals, such 
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as interactions with phonons, electrons, impurities, etc. Experiments give, however, 
a much larger value for the linewidth, indicating that the plasmon lifetime is due 
mainly to quantum mechanical effects arising in such a confined system, absent in 
Mie's treatment. For example, for Sodium h'^i ~ 20 meV, while the observed linewidth 
for Na clusters with a = 2.5 nm is more than 10 times this value. 



We have already described one such finite-size effect, see eq.(6.15) and fig. (5.2): 



the motion of the electronic center of mass linked to the excitation of the surface 
plasmon results in a net surface charge creating a time-varying dipolar electric field 
inside the cluster, which couples to the mean-field one-particle states calculated in the 
absence of such a field. Through this coupling, represented by the second term of the 
right-hand side of (6.15), the one-particle states become excited, absorbing energy 
from the field, and thus damping the surface plasmon. This decay of the surface 
plasmon into one-particle one- hole (Ip-lh) uncorrelated single-particle excitations 
which have the same energy as the collective mode is called Landau damping, just 
as the decay in bulk metals of the bulk plasmon into Ip-lh pairs [BB94, NP99]. The 
associated linewidth was first calculated using linear response theory [Kub57] and 
the fluctuation-dissipation theorem [CLOO] in [KK66,BS89], and can also be obtained 
applying Fermi's Golden Rule with the coupling hamiltonian in ( 6.15[ ), [WMWJ05] 



7 



27rA- 



ph 



M 



+ eh) 



(6.18) 



where \p) and \h) represent, respectively, particle and hole states of the mean-field 
problem. The result renders a plasmon linewidth with a characteristic 7(a) ~ 1/a 



dependence, observed in experiments as shown in fig. (6.2): 

^(a) - — ^ ( 

Aa \fyjJM 



(6.19) 



with e-p = h'^kp/2me and fp the Fermi energy and velocity, respectively. The mono- 
tonically increasing function qq is defined as 



^|v/x(a; + l)(4a;(x + l) + 3) - 3(2a; + 1) In + Vx + { 
^x{x - l){Ax{x - 1) + 3) - 3(2x - 1) In {^/^ + ^/x - l) 6 (x — V) \ , 



0(x) being the Heaviside step function. It is very well approximated by 

^o(a;) = (1 + (225x/64)2)-i/4 



(6.20) 



(6.21) 



In the most interesting size range 0.5 nm< a < 10 nm, where the transition from 
atomic to bulk behavior of physical properties takes place, the Landau damping 
mechanism dominates, limiting the lifetime of the collective excitation. 
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Figure 6.2: From [Dor65]. Experimental results for the inverse surface plasmon 
linewidth of silver nanoparticles embedded in a glass matrix. The quantity repre- 
sented on the abscissa is proportional to the radius a of the nanoparticles (between 
~4 and 10 nm), while the vertical axis shows the inverse of the linewidth, r = 7"^. 
It clearly demonstrates the linear size-scaling of the surface plasmon lifetime. The 
circles and triangles correspond to two different samples. 



For large clusters, a > 10 nm, radiation damping due to the interaction of the 
surface plasmon with the external electromagnetic field dominates our mode's atten- 
uation [KV95] . The lifetime associated to this mechanism can be shown to scale with 
size as '^rad ~ CL^ [KK66], explaining why it is dominant for large clusters, becoming 
negligible for smaller ones in comparison with the Landau damping, who presents 
a weaker dependence of 7 on a. For radii smaller than 0.5 nm the interaction with 
the ionic background becomes dominant, and the jellium model no longer provides a 
useful description. We will focus our attention in this thesis in the size range where 
Landau damping dominates. 



6.3.2 Landau damping and the double-counting problem 



The final form of the hamiltonian (6.15) is an almost perfect starting point for 
the description of the surface plasmon dynamics, in particular to address the Landau 
damping mechanism, with the plasmon (center-of-mass) degree of freedom coupled 
through He to the single-particle excitations, which can be viewed as a dissipative 
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environment for the plasmon. If the bosonic operator and the fermionic operators 
cj^ corresponded to well-defined, independent degrees of freedom, one could proceed 
without further considerations to calculate for example perturbatively the effect of 
the relative coordinates on the plasmon lifetime, but as we have already noted, the 
final form of Hrei is obtained forgetting the fact that there are not N but — 1 
independent relative coordinates. Thus, it seems that from original degrees of 



freedom in (6.1) we have arrived at + 1 after the mean-field approximation and 
second quantization procedure, double-counting in some way the plasmon degree of 
freedom in Hrei- From the point of view of the study of the plasmon dissipative 



dynamics, this spoils the virtues of our approximate hamiltonian (6.15), as we have 
to guarantee that the environment to which we couple our system does not include 
our system. 

Fortunately, previous knowledge coming from a related problem of nuclear physics, 
the description of the giant dipolar resonance of nuclei [RS80, YDG82], gives a hint 
on how to circumvent this inconsistency. This nuclear resonance corresponds to a 
collective motion of the (charged) protons with respect to the (neutral) neutrons, when 
an external electromagnetic field is applied, analogously to what happens in metallic 
clusters. Using a slightly different approach to obtain the approximate eigenstates 
of the nuclei, it was shown that the collective state is to a very good approximation 
a linear superposition of low-energy Ip-lh excitations of the Hartree-Fock ground 
state. The restriction to this kind of excitation, neglecting multiple p-h states, is a 
good approximation to describe the low energy physics of the system. The Ip-lh 
excitations are not true eigenstates due to the non-Hartree-Fock part of the Coulomb 
interaction, named residual interaction, which mixes them giving rise to more realistic 
eigenstates of the system, most of which are very similar to the ones obtained at the 
Hartree-Fock level, except for one linear combination, whose high energy lies outside 
the low-energy sector where its Ip-lh constituents are located, and which corresponds 
to the collective resonance. We will give a detailed derivation of these results in the 
following section. 

The surface plasmon is thus basically originated from a combination of low-energy 
one-particle mean-field states, showing up when the Coulomb interaction is included 
in its full complexity. Coming back to the description of its dynamics, if we want 
to calculate its Landau damping, keeping a clear separation between the surface 
plasmon and the degrees of freedom constituting its dissipative environment, the 
previous results suggest the following two-step procedure [YDG82, YB92]: 

1. Divide the Hilbert space spanned by the Ip-lh Hartree-Fock states into a low- 
energy (called restricted subspace, Sr) and a high energy (called additional 
subspace, Sa) subspaces. Introduce now the residual interaction in Sr, and 
obtain the eigenstates, one of whom will be the surface plasmon, having an 
energy lying in Sa- 



2. To calculate the lifetime of the surface plasmon due to Landau damping, just 
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apply Fermi's Golden Rule using the residual interaction coupling the plasmon 
to the quasi-continuum of states in Sa degenerate with it. 



In the context of the calculation of Landau damping using (6.1), the residual inter- 
action is given by He, and to avoid the mixing of the surface plasmon and what we 
consider its environment, we just have to determine a safe energy cutoff Aec below 
the surface plasmon frequency Upi such that (i) the value of ujpi is insensitive to the 
increase of Aec, showing that the relevant low-energy Hartree-Fock Ip-lh excitations 
building up the surface plasmon lie below it, and (ii) only Ip-lh excitations above 
Aec, which can be thus regarded as belonging to a subspace independent of the sur- 
face plasmon, are included in Fermi's Golden Rule sum. The latter excitations are 
by definition Sa, so that the linewidth has an expression 

7CX ^ \{p\z\h)\^5ihLUpi-ep + eh) (6.22) 

pheSA 

Here ph is a particle-hole excitation composed of a particle with energy ep and a hole 
with energy eh- In this calculation of 7 the concrete estimate for the value of Aec is 
irrelevant as long as Aec < ^pi [YB92], but the approximate theoretical description 
of other physical quantities might depend on the cutoff. For example, if one uses 



(6.15) to study the decay and decoherence of the surface plasmon, the timescales that 
characterize the dynamics of the electronic environment in 5*^ [WIWJ07] depend on 
the cutoff energy. Moreover, it has been shown in [WIJW06] that the environment- 
induced redshift of the resonance frequency depends logarithmically on the cutoff. 

We will in the following section derive an estimate for this cutoff based on sound 
physical arguments. 



6.4 The surface plasmon as a superposition of low- 
energy particle-hole excitations 

To show this fact and find an estimate for Aec we will forget in this section about 
the model based in a separation of center-of-mass plus relative coordinates described 



in the beginning of previous section, whose main result was eq. (6.15), and use 
instead the approach borrowed from the description of giant dipolar resonances in 
nuclei which has already been briefly introduced. 



6.4.1 Separable interaction ansatz 

We now obtain the simplest model which starting from a mean-field approach in- 
cludes the essential features of the Coulomb interaction originating the surface plas- 
mon. The starting point is the simplest, mean- field approximation of the hamiltonian 
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(6.1), where the one-body electronic states obtained come from a self-consistent cal- 



culation like the TDLDA one (fig. 6.1 ). The problem with this naive approximation is 
that it fails to describe the ability of the Coulomb interaction to create self-sustained 
collective excitations, like the surface plasmon we want to study. We have to introduce 
an additional component 6V{r) to the self-consistent potential V{r) which captures 
this feature of the Coulomb interaction. To obtain SV{r), we recall the fact that 
the occurrence of collective vibrational modes in a system governed by independent- 
particle motion can be understood in terms of the variations in the average one-body 
potential produced by an oscillation in the nanoparticle electronic density. Such vari- 
ations in the one-particle potential give rise to excitations in the electronic motion, 
and a self-sustained collective motion is obtained if the induced density variations are 
equal to those needed to generate the oscillating potential [BB75]. 

In our case, from the several collective modes susceptible of appearing in our 
system we focus our attention on the simplest, dipolar in character, surface plasmon. 
Starting from the knowledge of the existence of such a mode, the procedure to obtain 
the associated SV{r) is then straightforward: calculate the variation in the average 
one-body potential V{r) produced by an electronic density variation 6n{r) linked to 
the excitation of the collective mode. This was done in [WMWJ05], as follows: 

Assuming that at equilibrium the electron density is uniform within a sphere of 
radius a, ne(r) = neQ{a — r), a rigid displacement with a magnitude Z along the 
^-direction changes the density at r from ne{r) to ne(r — R) = ne{r) + 6ne{r). To 
first order in the field R = Ze^, 5ne(r) = — R ■ Vne(r) = Znf.cos{6)6{r — a), where 
6 = arccos(r ■ e^/lrl). The related change in the mean-field potential seen by the 
electrons is thus 

6V{r,R) = j dh'6n,{r')Vc{r,r') (6.23) 
where Vc{r, r') denotes the Coulomb interaction between electrons. The term 6H 



which has to be added to the mean-field hamiltonian (6.11) is then 

N 



5i/ = ^5y(R,r:) (6.24) 

i=l 

A simple expression for 6V{r, R) can be obtained, given the spherical symmetry of 
the problem, using a multipolar decomposition for the Coulomb interaction [Jac75], 
namely 

Attt) r 

bV(v, R) = Z — zQia - r) + — e(r - a) , (6.25) 

3 L r-^ J 



so (5V(r,R) is actually the same as R ■ Vf/(r), eq.(6.4). A further simplification 
is obtained if we remember that, even though the electronic wavefunctions spread 
outside the geometrical boundaries of the nanocluster, giving rise to small quantum 
corrections like the spill-out effect, the expected influence of the exact form of 5l^(r, R) 
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for r > a is very small. We can then simplify 5V by extending its dipolar character 
also for r > a, 6V{r, R) ^ (47rnee^/3)Z ■ z, with a resulting separable form for 6H 

(6.26) 

where Z = Zi/N has been used. The term separable applies to those interactions 



satisfying V{xi^Xj) ~ f{xi) ■ f{xj), f{x) being a generic function of x, like (6.26). 
The results that would be obtained without the use of this approximation would 
include 1/N corrections to the ones we will present, like the estimate for Aec, so for 
our purposes we can work using the separable ansatz. Expressed in terms of the 



mean-field hamiltonian basis \a), the second-quantized version of (6.26) reads 



SH =^ Xd^^dlpC^c^pCsc^ , (6.27) 

a/37i5 

where d^^ = (q;|2;|7) are dipole matrix elements and A is a positive constant. 



6.4.2 Linear response theory and Random Phase Approxi- 
mation for the surface plasmon 

We will analyze the response of the nanocluster to a time-dependent weak ex- 
ternal perturbation, and obtain this response in a self-consistent way, leading to the 
RPA theory of response. The result will be written in terms of a response function, 
yet to be defined, whose poles determine the energies of the excited states. We will 



particularize to the separable potential of the dipolar field (6.27), and see what are 
its effects on the energies of the self-consistent mean-field states. We will obtain a 
secular equation for those RPA eigenenergies, and show how to solve it graphically. 

We follow [BB94] and [RS80]. Suppose we know the stationary one-particle eigen- 



states \i) and eigenenergies Ci of the mean-field, time-independent hamiltonian (6.11), 
and now introduce a time-dependent external perturbation V^(r,t) = V^(f)e*'^* of a 
given frequency u, interacting with the system through its particle density n{r), 

Hint{t) = j dW,{r,t)n{i) (6.28) 



Then, the correction to the eigenstates of the system is, to first order in Hintit), 

m= ^^ + -^Y.^3\V.ir)\e){ ^ + -7-) e--'*/^ (6.29) 
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where (pi is the single-particle wave function of state \i). The particle density becomes 
now n(r,t) = |0'j(f,t)p = no(r) + 5ra(f)e*'^*, with 

occ ^ 
i j J J 

We now define the independent particle retarded response function for the density 
change at f induced by a potential field at r', as 

occ ^ ^ 

i j J J 

In terms of the response to an arbitrary perturbation potential can be calculated 
from the integral, 

5n(?,cj) = j (S'R%i,i',uj)V^{i'). (6.32) 

This is the independent particle response, but one more ingredient has to be included 
to construct the RPA theory of response. So far the only perturbation was the 
external field, but it is clear that the induced density oscillation 6n{r,uj) will cause 
the self-consistent field to oscillate at the same frequency as well. We can define the 
change of the mean-field potential with respect to small density changes in terms of 
the functional derivative 6V{r)/6n{r').ln the case of an external long-wavelength field 
we have seen that the main modification to the mean-field potential is the presence of 
a dipolar electric field, described in (6.27). In terms of 6V{r)/6n{r') the time-varying 
mean field is then 

6Vm = f di'^^6n(i') (6.33) 
J dn{v') 

Adding this potential to the external potential, we obtain the implicit equation for 
the self-consistent density. 



6n 



RPA{i,u;) = j d?i?°(f,?2,^){K(f2) + j rf?'^^5nKPA(f)} (6-34) 

It is more useful to express the self-consistent density in terms of the external field, 
defining in this way the RPA response function R^^^ (we will always work with 
retarded response functions, so it will be assumed in the following, and not mentioned 
explicitly as "retarded"), 

5nRPA{i,uj) = j di'R^''\i,i',uj)V,{i') . (6.35) 



Prom (6.34) and (6.35) follows the implicit equation that the RPA response function 
must satisfy: 

:>RPAf~ ~i .\ _ dO/- -/ , a , f j~ jz. dO/- ~ , ,\'^^(l'2) DflPA/ 



i?«^^(f, ?',.;) = R\i,i\u;) + / rf?2rf?3i?'^(f,f2,^)^v^i?^^^(f3,f ,^) (6.36) 
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If the external perturbation were static, u = 0, equations (6.34) and (6.36) reduce sim- 
ply to the Hartree-Fock self-consistent solution to the problem, where the occupied 
eigenstates contribute to the total effective one-body potential in the hamiltonian, 
which in turn changes the occupied wavefunctions, starting a cycle that ends when 
selfconsistency between eigenfunctions and total potential is reached. The RPA ap- 
proximation is therefore nothing but the small amplitude limit of the time-dependent 
mean field approach. 

The next step is to repeat the linear response calculations, imposing again (i)a 
self-consistent mean field approach for the total hamiltonian and (ii) small amplitude 



limit for the perturbation, but working in the basis and language in which ( 6.27[ ) is 



written. The derivation of the response becomes much more tedious and lengthy (see 
e.g. [RS80, BB94]), so we will state here the main points and results, the latter being 
very similar in structure to those we have obtained: 

Take the creation and destruction operators c^, and Cq, of the mean-field states of 
the unperturbed hamiltonian. The external perturbation is chosen to be a one-body 
operator, Vx(t) = V^e"*"^* + V^e^'^^ = '^kiVki(t)c\ci. As also the rest of operators 
in the hamiltonian will be approximated to self-consistent one-body operators, one 
assumes 

• First, that the ground state wavefunction |RPA(t)) at any time is a Slater 
determinant, imposing in this way selfconsistency. 

• Second, one assumes that the external field V^(t) is weak, introducing only 
small perturbations around the stationary Hartree-Fock ground state |HF) of 
the unperturbed hamiltonian, |RPA(t)) = |HF) + 6\ip{t)). 

The combination of Slater-determinant sort of solution plus small perturbations with 
respect to |HF) limits automatically the Hilbert space out of which |RPA(t)) is con- 
structed, to the stationary Hartree-Fock ground state |HF) plus the set of all its Ip-lh 
excitations cj,Cg|HF) (cg annihilates an occupied one-particle state of |HF), while cj, 
populates a previously empty one-particle state). This shows up in the so-called 
spectral representation of the RPA response function obtained 

nRPAt ^ v-/(0l4cpK)(^l4'^9'|0) {0\cl,Cg,\i){i\clcp\0)^ 

^MPV(^) = nu-e, + eo hu + u-e, ) ' ^^'^^^ 

i>0 

where the index pairs pq and p'q' run only over particle-hole (ph) and hole-particle 
(hp) pairs, while all other matrix elements vanish in RPA order. Unsurprisingly, 



Rp^y satisfies an implicit equation analogous to (6.36) 



TjRPA _ pO I \ pO J J* tdRPA no\ 

^pqp'q' ~ ^pqp'q' ^ / j -"-pgpi qi "'P192 "p;.qi ^P9.q9.p'q' {V.OO) 



P2q2 
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In (6.38) we have introduced the spectral representation R^qpiqi of the response func- 
tion of the free system, 



E.>oK^l4c,|0)p-|(.|ctc,|0)p 
huj - ep + eg 



(6.39) 



and we have also particularized to our case, where Vx{t) is the dipolar field represented 



we can proceed further with (6.38) to get 



by the corresponding operator (6.27). Thanks to the separable form of the interaction. 



RPA/ 



R 



RPA, 



R^(u) 



UJ 



1 - Ai?o(cu) 



(6.40) 



with 



pqp'q' ph P 



+ eh hhj + t„- th 



(6.41) 



The last sum runs over Ip-lh states, \ph), whose energy with respect to the RPA 
ground state is Ae^/i = Cp — Ch- The poles of R^^^{ijj) give the excitation energies E 
we are looking for, which are thus defined by the secular equation 



S{E) 



(6.42) 



Here we have defined the RPA sum 



S{E) 



'^^£ph\dphY 



ph 



E^ 



^^ph 



(6.43) 



Nuclear physics textbooks (see, e.g., Ref. [RS80]) show how to solve the secular 



equation (6.42) graphically under the implicit assumption that the p-h spectrum 



is bounded in energy. Basically one plots S{E) and a horizontal line at the ordinate 
value y = 1/X. The intersections of S{E) with the horizontal line give the solutions 
E, which turn out to be merely small renormalizations of the p-h energies Asph ex- 
cept for the largest excitation energy Spi which corresponds to the collective mode. 



see fig. (6.3) for a schematic example. This state is nothing but the surface plasmon 



collective state. In fig.(6.3D-c) it shows how Epi is a function mostly of the initial 
low-energy Ip-lh excitations, depending only slightly on the high-energy ones, for a 
decaying dipole matrix element dph ~ l/Ae^^. This decaying behavior will be demon- 
strated in next subsection, ensuring also the applicability of this schematic picture 
even though in a metallic cluster the p-h excitations are not bounded from above on 
the scale of the plasma energy. 
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Figure 6.3: (a) Schematic example of the graphical solution of a secular equation of the 
kind (6.42). In this case the dipole matrix elements are assumed to have a dependence 
with Asph, dph ~ l/ASpf^, which corresponds to the reahstic case of metallic clusters. 
1/A is chosen to be 1/A = 0.002 (arbitrary units), while the sum includes only 4 Ip-lh 
excitations, whose excitations energies in absence of the dipolar residual interaction 
are Asph = 3,7,11,15 (arbitrary units): S{E) = {l/3)/{E^ -3^) + {l/7)/{E^ - 
rY+ (1/11)7(^2 - 112) ^ (i/i5)/(^2 _ ^52^_ jy^^ ^Yie residual interaction the 
eigenstates and eigenenergies change, with the new excitation eigenenergies E given 
by the crossing points where S{E) = 1/A. The highest excitation energy corresponds 
to the surface plasmon collective mode, Epi. (b) and (c) reflect the influence of high 
and low Ip-lh excitations in the value of Epi, respectively; in (b) Asph = 15 is absent, 
causing a slight shift of the surface plasmon energy Epi, while in (c) Asph = 3 is absent, 
resulting in a strong shift of Spi. This points to the predominant role of low-energy 
p-h excitations in the formation of the collective mode. 



6.4.3 Fast decay of dp/^. Plasmon state built from low-energy 
p-h excitations 



The secular equation (6.42) depends crucially on the form of the dipole matrix 



element dph- We will show now more rigorously that due to its fast decay dph{Ae) ~ 
l/Ae"^, cf. eq. (6.57), the plasmon excitation energy depends only on the lowest Ip-lh 



excitations. Assuming that the p-h states are confined within a hard-wall sphere, one 
can calculate dph, see Appendix D.l for a detailed derivation, dph is factorized into 
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radial and angular parts [YB92] 

The angular part is expressed in terms of Wigner-3j symbols as 



(6.44) 



L Ih 1 



- (-1)- V(2^. + + 1) X - ^ J ^ ^ J (6.45) 

and sets the selection rules for the total and azimuthal angular momenta, respectively: 



— ^/i ± 1 



nir. 



rrih 



The radial part depends on the energies of the p-h states as 

2h^ y/e^ 



(6.46) 



(6.47) 



With the help of equations ( 6.44 )-( 6.47) and using the aforementioned dipole selection 



rules (6.46), we can write the RPA sum (6.43) as 

2 



SiE) 



(6.48) 



nh,lh 



where and Up are radial quantum numbers. We have defined //^ = //i+l if lp = + 1 



and flp = Ih iilp = Ih — ^- Using the semiclassical quantizations (D.24) and (D.25) of 



Appendix D.2 we finally obtain the result shown in Figure 6.4 



In the figure we show (solid black line) the resulting i?-dependence of (6.48). The 
crossings of this curve with the horizontal dotted one with height 1/A yield the exci- 



tation spectrum within the RPA. The lowest p-h energy estimated in Appendix D.2 
is 

Ae^i, ^ . (6.49) 



The existence of a gap in the excitation spectrum of the dipole operator stems from 
the quantization of the states in a hard-wall sphere and the dipole selection rules. 
Whenever the energy E coincides with a p-h excitation energy Asph, we have a 



divergence in S{E) (see vertical grey lines in Fig. 6.4). For the lowest energies E 
of the interval considered the sum is dominated by the term associated with the 
divergence closest to E. On the other hand, for the largest energies E, the fast decay 
of dph with Asph means that the divergences are only relevant for energies extremely 
close to them. Away from the divergences the sum is dominated by the contributions 
arising from the low-energy p-h excitations. 
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Figure 6.4: RPA sum (6.48) for a Na nanoparticle with k^a = 30 (solid black line). 
The p-h excitation energies Asph (represented by vertical grey lines) have been ob- 
tained from the semiclassical spectrum (D.25). There are about 15000 degenerate 
excitations in the interval shown, and the smallest excitation Aemm — Bp/ 10 has a 
degeneracy A/" = 380. The dashed line takes only into account the contribution of 
Aemin (see Eq. (6.60)). The cutoff energy A^c separating the two RPA subspaces 
is also shown in the figure (see Eq. (6.66)). Above Asc the two sums are close to 



each other, showing that for high energies the RPA sum is essentially given by the 
contributions coming from the low-energy p-h excitations. The horizontal dotted line 
indicates the position of the coupling constant A entering the RPA secular equation 



(6.42), according to the estimate (6.61). 



To gain physical insight and proceed further in the analysis of the electron dy- 
namics in metallic nanoparticles, we introduce the typical dipole matrix element for 
states separated by a given energy difference Ae, 



d^-^iAe) 



iP-h 



7^ \dphf^{Ae-Aeph) 



1 1/2 



(6.50) 

Here, we have introduced two quantities, pP"'^(Ae) and C{e, Ae). C{e, Ae) is the local 
density of dipole matrix elements 



C{e, Ae) = J2 \dph\^^{£ - ^h)5{e + Ae - ep) 

ph 



(6.51) 
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As shown in Appendix D.4 this can be expressed as 

1 a2 



C{e,A£) 



where 



37r2A£2 \AeJ ' 
F{x) = {2x + 1) a/x(x + 1) - In (^^/x + \J x + 



For Ae <^ e, equation (6.52) simphfies to 



C(£,Ae) 



2a2 

37r2 Ae4 ' 



In expression (6.50) 



ph 



(6.52) 
(6.53) 

(6.54) 

(6.55) 



is the density of p-h excitations with energy Ae respecting the dipole selection rules. 
In order to simplify the presentation, we do not consider spin degeneracy factors. In 



Appendix D.5 we show that for Ae Sp we have 



p^-\Ae) 



{k-pa)'^ Ae 
47r2 ll' 



(6.56) 



and therefore in such a limit the typical matrix element (6.50) can be approximated 

by 



d^-\Ae)c.^-y^(^Y 



(6.57) 



To check the validity of our estimate of the typical dipole matrix element, there is 
an exact sum rule J2ph^^ph\d'^'^\'^ = {?>/ ^Ti)h^N /2m^ [LS91,YB92], satisfied also in 
the presence of Coulomb interactions between electrons. Performing the sum using 



(6.57) 



ph 



Ef 



dAeAe 



Ef 



deC{e, Ae) 



-Ae 



(6.58) 



one obtains about 70% of the exact result {3/A7[)h'^N/2mc [Bra93]. This is quite 



reasonable regarding all the approximations we made to obtain (6.57). 



In order to emphasize once more the importance of the low-energy p-h excitations. 



we present in Figure 6.4 (dashed line) the contribution to the RPA sum coming only 
from the infrared p-h excitation energy with the appropriate degeneracy factor TV 



see Eq. (D.27) in App. D.2). Indeed, we can estimate (6.43) as 



S{E) ^AT 



2A£^i„ [d^-\Ae, 



(6.59) 
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With the results (|6.49D, (\6.57\j, and (|D.27|), we obtain 



S{E) 



64A;p2 



- (epTr/ZcFa)^ 



(6.60) 



While in the lower part of the energy interval in Figure |6.4| the two curves exhibit 
considerable discrepancies, in the second half of the interval they are very close (except 
of course at the divergences). Since the collective excitation is found in this last 
interval, we see that it is mainly the low-energy p-h excitations that are relevant for 
the definition of the collective excitation. 

Since the resonance energy is known from experiments, we can obtain the value 
of the coupling constant A. Indeed, using the estimate (6.60) evaluated ai E = Itwm 
in the secular equation (6.42), we obtain 



(6.61) 



to leading order in k-pa ^ 1. Notice that this result is consistent with the estimate 
obtained from the energy- weighted sum rule given in [YB92]. For the case studied 
in Figure 6.4, our estimate (6.61) yields kpSp/X ~ 1600. As the radius a increases 
the lowest p-h energy decreases like 1/a, the degeneracy at this value increases as a^, 
and the density of p-h excitations contributing to (6.43), p^^{Ae), grows as a'^. This 
increase in the number of excitations contributing to the sum is partially cancelled by 
the 1/a behaviour of the typical dipole matrix element, resulting in a decrease of the 
coupling constant A proportional to the number of particles in the cluster, A ~ 1/a'^, 
and a value of the plasmon frequency which remains almost unaffected. Therefore, 
for a larger nanoparticle size, the divergences shown in Figure 6A would be more 
dense in energy, starting at a lower energy, and the vertical scale would increase as 
a^, thus obtaining a similar value for the plasmon excitation energy huu- 



6AA Separation of the reduced and additional particle-hole 
subspaces 

As we have shown in the last section, the high-energy part of the p-h spectrum is 
not crucial for the determination of the energy of the collective excitation. In what 
follows we make this statement more quantitative and estimate the upper-bound 
cutoff A^c of the low-energy excitations that we need in order to obtain a stable 
position of the surface plasmon. 

In order to obtain a quantitative estimate of the cutoff, we require that by changing 
it from A^c to (3/2) A^c, the position of the plasmon changes only by a fraction of its 
linewidth 7, the smallest energy scale with experimental significance. Our criterion 
leads to the condition 



(6.62) 
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with the RPA sum S that has been defined in (6.43). The additional subscript refers 
to the upper bound of the p-h energies. 



The left-hand side of (6.62) can be estimated according to 

2Ae [dP-^\Ae 



Aec 



(6.63) 



with dP~^ and p^"^ as defined in (6.50) and (6.55), respectively. Using (6.56) and 



(6.57), we obtain 



MJ 



4fcp2 / By 



M 



As mir 



As, 



(6.64) 



to leading order in k-pa. Similarly, 



37r2 



k-pa 



M 



27 



Ae 

mir 



3A£, 



(6.65) 



Using the expressions (6.19) and (6.49) for 7 and A^min, respectively, finally yields 



according to the criterion (6.62) the main result of this chapter, an estimate for the 



cutoff energy based on sound physical arguments 



Ae,, 



vr 



fVjJ 



M 



with the function qq defined in (6.20). For Na clusters, we have e-p/hjjyi 
our criterion yields a value A^c — (3/5)£:f- 



(6.66) 



0.93 and 



We have verified the robustness of our criterion (6.62) by exploring different phys- 



ical parameters, like the size of the cluster. Indeed, it can be seen in Figure 6.4 that 
above As, the RPA sum evaluated from ( 6.48[ ) (solid line) and the estimate (6.60) 
(dashed line) are close to each other, showing that for high energies the RPA sum is 
essentially given by the contributions coming from the low-energy p-h excitations. We 
have checked that this feature of A^c is independent of the size a of the nanoparticle. 

In order to have a well-defined collective excitation, the cutoff energy Ae, must 
obviously be larger than the minimal p-h excitation energy (6.49). For Na, we find 
that this condition is already verified with only = 20 conduction electrons in 
the nanoparticle, in agreement with the experiments of reference [SKM"''91] and the 
numerical calculations of reference [YB91]. 

Eq. (6.66) is the key ingredient we needed to avoid the problem of the clear 



separation between the plasmon and its environment, as discussed in section 6.3.2 
The splitting in low Asph < A^c and high-energy Asph > Aec p-h excitations gives 
physical soundness to the method based on the separation of the electronic degrees 
of freedom into centre-of-mass and relative coordinates [GGI02, WIJW06, WMWJ05, 
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WWIJ07] described in the first two sections of tlie cliapter. Now tliat tfie Fermi's 
Golden Rule calculation of the Landau damping of the plasmon (6.18) is justified, 
provided one sums over p-h pairs with Asph > Aec, we will use the estimate (6.66) 
to compare in next section the characteristic time scales of the surface plasmon and 
its high-energy p-h environment. We will show that the memory of this environment 
is very short compared to the timescale characterizing the evolution of the plasmon, 
its decay time 7(0)"^, justifying in this way a Markovian approach for the analysis 
of the plasmon dynamics in presence of such a dissipative bath, cf. [WIWJ07]. 



6.5 Dynamics of the relative-coordinate system 



The decomposition (6.15) of H suggests to treat the collective coordinate as a 
simple system of one degree of freedom which is coupled to an environment with 
many degrees of freedom. The latter are the relative coordinates described by -ffrei- 
This is the approach taken in [WIWJ07]. In this picture, the time evolution of the 
center-of-mass system (i.e., the surface plasmon) strongly depends on the dynamics 
of the relative-coordinate system. Such a dynamics is characterized by a correlation 
function which can be written at zero temperature as [WIWJ07] 

C{t) =A'Y1 Kh\'e'^''"''/^&{Ae,H - A£c) • (6.67) 

ph 

C{t) is nothing but the Fourier transform of the spectral function of this dissipative 
bath [LCD+87,Wei99] 

S(A£) = Yl \dph\'S{Ae - AepHMAsph - Ae,) (6.68) 

ph 

which has been calculated in [WIJW06] . Finite temperatures were shown to result in 
a small quadratic correction. Consistently with the results of the preceding sections. 



we employ our low energy estimates (6.56) and (6.57) to obtain 

4a \ Ae ^ 

This result is consistent with the one of reference [WIJW06] (see Eq. (34) in there) in 



E(A.)^^*(^) e(A.-A..). (6.69) 



the limit As <^ huM and for zero temperature. A schematic plot is given in fig. (6.5), 
showing the dominant role that the p-h excitations with energies close to Aec have 
in the properties of the environment, pointing to a correlation time Tcor of the bath 
of order ~ h/Asc- 

In principle we could calculate C{t) by taking the inverse Fourier transform of 



(6.69). The decay of C{t) for very long times is dominated by the discontinuity of S 
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Figure 6.5: Schematic plot of the spectral function ^(Ae), Fourier transform of 
C(t), with A^c taken to be A^c = O.le^- The sharp peak is a consequence of the fast 
decrease of dP~^{Ae) with Ae, eq. (6.57), which overwhelms the increase in the density 
of particle-hole excitations contributing to the damping of the plasmon, eq.(6.56), 
giving in this way a preponderant role in the dissipative bath to the particle-hole 
excitations with energies close to A^c- 



at A^c- This is somehow problematic since A^c can only be estimated as we did in 



Section |6.4.4[ and since the functional form of the long time decay depends on how 
sharply the cutoff is implemented. However, it is important to realize that it is not the 
very long time behavior that determines the relevant decay of the correlation function, 
but rather the typical values at which C{t) is reduced by an important factor from 
its initial value C(0). We then estimate the correlation time as the mean decay time 
of C(t), 



(''~cor) 



"^"dtUw 



1 



C(0) 



dtC{t) 



(6.70) 



If C{t) were an exponentially decreasing function, (rcor) would simply reduce to the 



inverse of the decay rate. Using the definition (6.68), we get 



(''"cor) 



/■oo 

n / dAeS(Ae)/Ae 

JAEc 

POO 

/ dAeS(Ae) 



(6.71) 



Given the fast decay of the function S ~ the above integrals are dominated 
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by their lower limit A^c and we have 



(l"cor) 



2 h 



(6.72) 



Since Ae^ is of the order of {3/5)ep, we see that the response time (or correlation 
time) of the electronic environment is of the order of its inverse Fermi energy. 

The estimate of the characteristic response time of the electronic environment is 
crucial in justifying the Markovian approximation used in [WIW J07] . In that work, 
the degrees of freedom corresponding to the relative coordinates were integrated out 
and treated as an incoherent heat bath that acts on the collective coordinate. Such 
an approach relies on the fast response of the environment as compared to the time 
evolution of the surface plasmon. The typical scale for the latter is the inverse of the 
decay rate, Tgp = 1/7. Using (6.19) and (6.66), we therefore have 



' sp 



TT 



M 



(Tcor) 9[go{eF/hwM)]^ £f 



(6.73) 



For the example of Na nanoparticles worked in Section 6.4.3 , we have Tsp/ (rcor) — kpa. 



This is a safe limit since in not too small nanoparticles, kpa 3> 1. As the size of the 
cluster increases, the applicability of the Markovian approximation is more justified. 
This is expected since the electronic bath has more and more degrees of freedom, 
approaching an "environment" in the sense of quantum dissipation. Since the physical 



parameters of alkaline nanoparticles entering (6.73) are close to that of noble-metal 



clusters, the dynamics of the surface plasmon can be expected to be Markovian in 
that case too. 



6.6 Conclusions 

We have studied the role of particle-hole excitations on the dynamics of the surface 
plasmon. A key concept in this analysis is the separation into low-energy excitations 
which lead to the collective excitation once they are mixed by the residual dipolar 
interaction, and high-energy excitations that act as an environment damping the 
resonance. Using the random phase approximation and assuming the separability 
of the residual interaction, we have established a criterion for estimating the cutoff 
energy separating the low- and high-energy subspaces. The resulting cutoff energy is 
approximately (3/5)£f for the case of Na nanoparticles. 

Since the number of electrons in the cluster is finite, the assumption that the high- 
energy particle-hole excitations act on the collective excitation as an environment, 
introducing friction in its dynamics, may be questionable. What settles this issue is 
the ratio between the typical evolution time of the collective excitation and the one 
of the high-energy particle-hole excitations. The former is given by the inverse of 
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the plasmon linewidth, while the latter is obtained from the decay of the correlation 
function of the environment. We have found that this ratio improves with increasing 
cluster size. Even for a small cluster with a = 1 nm, the ratio is approximately 10, 
justifying the use of the Markovian approximation which assumes a fast time evolution 
of the environment with respect to the one of the collective excitation. 

The relevance of memory effects in the electronic dynamics of small clusters is 
of current interest, due to the advance in time- resolved experimental techniques 
[BHMDOO, FVF+00, LKLA99, LUCSOl]. First-principle calculations have recently 
addressed this issue by comparing time-dependent density functional theories with 
and without memory effects for small gold clusters [KB06]. For very small clusters 
(N < 8) memory effects were shown to be important. It has to be noted, how- 
ever, that the memory considered in [KB06] is that of the electron g whole, 
while we are concerned in this work with the memory arising from the dynamical 
evolution of the relative-coordinate subsystem. It would be interesting to consider 
cluster sizes intermediate between the ones considered in the present work and those 
of reference [KB06] in order to study the emergence of memory effects. 
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Appendix A 
Appendix to Chapter 



1 



A.l Modes' equations of motion of quasi- ID res- 
onators. 

A. 1.1 Elasticity basics. 

We will follow ref . [LL59] . The starting point is to define the displacement vector 
of a given point Xi in the solid, representing the difference between 

the coordinates of the same point before and after the deformation. The related 
magnitude that appears everywhere in the theory of elasticity is the so-called strain 
tensor (symmetric by definition): 

^ _ 1 _^ ^ _^ 1 _^ ^\ ^^ 

2\dxk dxi dxidxkJ 2\dxk dxiJ 

Then the stress tensor is defined, considering the forces Fi exerted on a small 
volume dV due to the deformation of a body: 

FdV = f ^dV = la,,dS, - = ^ (A.2) 
J oxk J dxk 

With these two quantities one can study the thermodynamics of the deformation: 
take a deformed body and suppose there is a small variation in the deformation, so 
that Ui Ui + 6ui, then the internal work done by the internal forces is given by 

'^du.dV (A.3) 

OXk 

We will consider bodies at a fixed temperature T, so the thermodynamic potential we 
will use is the free energy F = U — TS, with a variation of F in a small deformation 
at a constant T given by dF = aikduik, implying aik = dF / duik\j,- Now, for small 
deformations, choosing as non-deformed state the one in absence of external forces, we 
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can expand F as a function of Uik, obtaining F = ^XikimUikUim, so that aik = XikimUim 
(Hooke's law). For an isotropic body the expression reduces to F = Fo+(A/2)u|j+/iu^j;. 
(A and fi are the Lame coefficients), which can be reexpressed as 

/ I \2 K 2 
F = fiyuik - -SikUuj + , K = \ + -n (A.4) 

defining the rigidity modulus K. If the deformation is homogeneous (the strain ten- 
sor is constant throughout the body), and more specifically a simple extension or 
compression along the z axis, one has 

a,^ = E X u^^ , (A.5) 

defining the Young modulus, 

^ = 

iK + /i 

The ratio between transversal contraction Uxx and longitudinal elongation Uzz is Pois- 
son's coefficient 

v = — — = . (A.7) 

A. 1.2 The case of a rod. 

We start defining a local coordinate system {^,riX), attached to the point, that 
coincides with [x, y, z) in the absence of deformation (the rod is then along the Z axis). 
Next we define Q = dcp/dl, the rate of variation of the relative rotation angle between 
the local coordinate systems of two adjacent points {dl is taken along the bar's main 
axis). For example, a pure bending has Q = Q^jV- Now we define a coordinate system 
fixed in space, independent of the rod, with respect to whom a point in the rod will 
be given by a vector r. The tangent vector at that point is given hj t = df/dl, and 
the curvature by dt/dl, with modulus 1/R, defining the curvature radius. For a slight 
bending one has 

dt d^f / d^Y d^X 



In terms of Q the free energy of a rod is deduced to be 

^ f {^hE^I + hE^l + Cnl^l (A.9) 



^ rod 2 

where Ji and I2 are the principal inertia moments of the section of the rod, and C is 
the so-called torsional rigidity. In the case of a slight bending the expression reduces 
to 

F..„..«=2/{/i(^) +/.(^) }dz (A.10) 

Now we will deduce the equations for the modes as slight deviations from the equilib- 
rium situation, and the first step towards our goal is the analysis of the equilibrium 
equations. 



A.l Modes' equations of motion of quasi- ID resonators. 
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Equilibrium conditions. 



If there are external forces F^^t applied to the rod, and external moments M^r^t 
associated with them, the rod will suffer a deformation until the internal forces and 
moments compensate the external ones, so in equilibrium one must have 



Fext ~l~ Fint 











Considering the case of slight bending, for which Vt is given by eq.(A.8), the internal 
moment Mint associated to the internal stresses is Mint bending = {—EIiY", EI2X", 0), 
where the double prime represents d'^/dz'^. Taking a small volume of length dl and 

— # — * 

supposing an external force per unit length along the rod Fext = K = {K^, Ky, K^), 
the first of the equilibrium equations results in dF/dl = —K, and the second in 
dMint + [dl A K] = ^ dMint/dl = F At. Deriving the second equation with respect 
to I, and taking into account the small value of dt/dl (slight bending), one arrives at 

d'^Mint/df = tA K, so that 

EI2 d^X/dz^ = Kx 



Eh d^Y/dz^ 

and the internal forces will be given by 

Fx = -EhX" 
Fy = ~EhY" 



(A.ll] 



Mode's equations. 

To obtain the equations of motion of an elastic medium, we start from Newton's 
equation equating the product of mass and acceleration to the internal forces the 
mass is subject to 

pu. = — (A.12) 



In the simplest case of longitudinal vibrations we have, using eq.(A.5) 



da 



pUz 



zk 



P d^Uz 



(A.13) 



dxk dz"^ E dt^ 

This is the usual wave equation with linear dispersion relation. For the bending 
(fiexural) case, the equations of m otion are obtained just by substituting —K^ and 
—Ky in the equilibrium equations (A.ll) by pSX and pSY (S = transversal section 
of the rod): 



EL 



dz^ 
d'Y 



dz^ 



-pS 



d^Y 
dt^ 



(A. 14) 
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These equations admit plane wave- kind of solutions X{z,t),Y{z,t) ~ e''-'''-'^^\ but 
with a quadratic dispersion relation, ujj{k) = ^jEIj/{pS)k'^. In the case of torsional 
vibrations, the equations follow from equating the moment C dVL(^/dz to the temporal 
derivative of the angular momentum C per unit length, C = pld(f)/dt. In the last 
expression is the rotation angle of the considered section, so that d(j)/dz = Q(^, and 
I is the inertia moment of the transversal section with respect to the center of mass, 
I = f (x^ + y'^)dS, leading finally to 

Again a wave equation with linear dispersion relation for the variable 0. 

Shape and frequency of the bending eigenmodes 

Take a rod of length L, thickness t, width w, and mass density p. There are several 
possible configurations of a rod; we choose a doubly clamped beam, with both ends 
fixed. In the experiment the fundamental mode of transversal vibration in one of the 
directions, lets say X, is excited, and its frequency uq is measured. Now we will solve 
the equation for X{z, t) following [?] with fixed ends for the rod, and obtain a relation 
between ujq and the Young modulus E, from which we will calculate its value. 
We try a solution of the form X{z,t) = Xq{z) x cos(co't + 6). Inserting it in the 



eq.(A.14) for X{z,t), we obtain an equation that Xq{z) has to verify: 



d'^XQ{z) 2 V t \ '^v ( \ ( \ ^R\ 

— j-^ — = —Xq{z) = K Xo{z) (A.16) 

The solution to this equation is of the form 

Xq{z) = A cos{k,z) + B sin(Kz) + C cosh(Kz) + D sinh(Kz) (A. 17) 

Now, the boundary conditions we have to impose in this solution are, for fixed ends, 
Xq = = in 2; = 0,L. The final form of Xq(z) in terms of k, and the equation 
for the possible values of are 

Xq{z) = y4|(sinKL — sinh kL) (cos — cosh/tz) — (cos/tL — cosh kL) (sin — sinh/t2;)| 

1 = cos cosh (A.18) 

The possible values for n can be estimated to be, for kL > 1, k ~ ^^^f ? because 
the cosh kL takes big values that have to be "compensated" by the cos/tL, for their 
product to be 1. Therefore the density of states in the k-space is approximately 
constant, ^ ^ ^. The first solution for k is — ^"^^'^ , already close to and 

from the relation between k and to {(^{k) = \/ x k^) we get 




^o = ^x-^ (A.19) 
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B.l Some details about the Standard Tunneling 
Model. 

B.1.1 Distribution function of the TLSs 

From ref. [Phi87]. The tunneling model proposes the existence of an ensemble of 
TLSs, each one characterized by a hamiltonian: 

i/o = Aga, + Aga, (B.l) 

The task is to determine the probability density distribution P(Ag, Aq) of the TLSs 
in terms of Ag and Aq. A simplified expression for the tunneling amplitude as a 
function of the potential barrier height V and width d is 

Ag = me-^ = nne-'^(2mV/?i2)i/2 2) 

where HQ is approximately equal to the average of the ground states of both isolated 
wells, and m is the particle mass. 

Typical values for A can be estimated demanding that Ag must be approximately 
equal to kT if the tunneling state is to contribute to thermal properties at a temper- 
ature T. At 1 K this requires a tunnel splitting of 10~^ eV, which with hfl ~ 10~^ 
eV gives approximately A ~ 5. This is equivalent to a bare proton tunneling across a 
barrier of 0.1 eV with d = 0.7 A. 

For Aq it is argued that the distribution function must be symmetric because both 
positive and negative values of Aq are equally likely. The scale of energy variation 
is determined by the thermal energy available at the glass-transition temperature Tg 
where the fluctuating local potentials of the liquid are frozen in the structure. Since 
Tg is between 200 and 1000 K for most glasses, this energy is about 0.05 eV, much 
larger than the thermal energy available at 1 K. The low temperature properties are 
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therefore sensitive to the center of a broad symmetric distribution, so that -P(Ag, Aq) 
can be taken as independent of Aq. 

In the case of the dependence on Aq, due to the exponential dependence of Aq on A, 
only a relatively small range of A's is sampled for a large range of Aq and over this 
limited range the distribution of A can be assumed constant. This leads immediately 
to P(Ag, Aq) = Pq/Aq. This result is slightly modified by a logarithmic factor if the 
distribution of A varies slowly with the energy, but we will ignore this dependence. 
The limits of the possible values for Aq and Aq, and the value of Pq, can be deduced 
from experiments (see next subsection). 



B.1.2 Determination of values for Pq and A* 

The value Pq ~ 10^^J"^m"^ is obtained by the experimentalists through experi- 
ments of heat release in amorphous solids, as follows (the derivations in this appendix 
follow very closely the explanations appearing in ref . [Esq98] ) : An amorphous solid is 
taken which is at a temperature Ti, and cooled down very fast to a temperature Tq. 
The slower TLSs, which were, as all the rest, in thermal equilibrium at a temperature 
Ti, do not have time to follow the cooling of the system, so they equilibrate with the 
surroundings releasing part of their energy and heating up the rest of the solid, until 
it reaches another temperature T2 > Tq. One can calculate the time dependence of 

the heat release as Q{t) = Ni{t)e = {l/2)N{t)e, where is the number of TLSs 
in the excited state, e is the energy of such a state, and N = — N-^. Then, as- 
suming that the dynamical behavior of N{t) can be described by the relaxation time 
approximation 

d[Nit) - iVo(T)] Njt) - iVo(T) 

It = ^ ^ ^ 

where t(T) is the relaxation time of a TLS with energy difference e at a temperature 
T. Taking T =const., one arrives at 



SO that the heat release is given by 



1 / e e \ e ^(^0) 

QW = -rf(ta„h — -tanh— )^ (B.5) 

Replacing the total number N of TLSs by the distribution in terms of e and u {e = 
a/A^ + and u = A^/e) and the volume V of the sample, the result is 

PqV e , , e \ , 1 



(B.6) 
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Notice that t{Tq) is a function of e and u. For low temperatures t{Tq) can be 
approximated by the FGR resuh 

so that 

Qit) = / c/eeATf tanh — — -tanh — — ] • / drpGR 77^. 

(B.8) 

For usual experimental conditions and most glasses the last integral is insensitive to 
the limits of integration, and moreover, it can be shown that the main contribution to 
the heat release comes from TLSs with e ~ IksTi and t ^ t. For the case e <^ £max, 



leads to 

TT 



Q{t) = V(T2 - T2)-e-*/-— (B.9) 

For t < 0.1 Tmax the result is independent of Tmax and we obtain 

Qit) = '^PoV{T^-T^)^- (B.IO) 

Fitting the heat release experiments to this formula determines Pq. There are other 
types of measurements that can be used for this purpose, but this is the most accu- 
rate one, as the coupling constant 7 characterizing the coupling of the TLSs to the 
acoustic waves does not appear. 

Regarding the upper cutoff of the distribution of TLSs, A*, again the experiments, 
in this case sound velocity ones, fix a minimum of about 5 K for oxide glasses, reason- 
ing as follows: The interaction of elastic waves with TLSs results in a relative change 
of the sound velocity 

v="i:^''V). (B.ii) 

where x{^) = x'{^) + '''X"{^) is the susceptibility of a single defect. At low tem- 
peratures the relaxation rates of thermal TLSs are much smaller than the exter- 
nal frequency uj, ojt ^ 1; then only the 'resonant' part of the susceptibility, x' = 
(A^/e^) tanh [e/2/cT] , is relevant. Summing according to the distribution P(Aj., A^) = 
Pq/Ax, the result is 

— = ^ \og{T/Tref) (F < < knT/h) (B.12) 

where T^e/ is a reference temperature. This result is used to fit the experiments on 
the low-temperature sound velocity of various glasses. The logarithmic law arises 
from thermal TLSs with small bias A^, i.e. with A^. close to temperature, and due to 
the fact that this behavior has been observed up to several K, it has been concluded 
that A*/^^ is at least about 5 K. 
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B.1.3 Predominant coupling of the strain to the asymmetry 

Given the energy of a TLS as e = (A^ + A^)^/^, the couphng to strain fields e 
may be represented by the deformation potential 



^ de A^dA^ A^dA^ , 
D = — = ——^ + ——^. B.13 
de e de e de 



In the standard tunneling model the second term is neglected, with a resulting cou- 
pling 

de Az dAz ,^ 
D = — ^——^. B.14 
de e de 



Experimental support is given in ref . [And86] , where the following considerations are 
made: As there is a distribution of A^. and A^, there will be a distribution n{e,D) 
in D. Due to the coupling D the TLSs infiuence the low-temperature behavior of 
the volume thermal-expansion coefficient (3, (3 oc no{D), and ultrasonic attenuation 
a, a oc r2o(-D^) (the averages are over the distribution of TLSs). From attenuation 
measurements (Z)^)^/^ !^ 1 eV, and from thermal expansion measurements (D) ^ 
10~^ eV, so n{D) is a very broad and almost symmetric distribution. Now, to justify 



eq.(B.14) the author first proves that (-D^) and uq are almost independent of pressure 
P, for values of P up to 1000 atm and more. This, together with the observed linear 
variation of a with temperature, indicates that the distribution n{e, D) is essentially 
independent of pressure. This is needed in order to get information from experiments 
of ultrasonic attenuation where a pressure is applied at low temperature (0.5 K). The 
TLSs affecting the propagation of sound are mainly those with energies close to 0.5 
K, and as increasing pressure is applied, TLSs which were in this region get out from 
it, and others, which were originally out, replace them. Given the high pressures 
exerted (up to 1300 atm), and the value of (-D^), the lack of change observed in 
the population of TLSs around 0.5 K can only be explained if, for a TLS, when a 
pressure is applied (equivalent to a variation of local strain), the energy as a function 
of pressure first decreases and then increases, id est, D has to change sign at some 
point as P increases. The most natural explanation is to assume that the main effect 
of pressure is to change the relative positions of the two wells of the double-well 



potential (B.14), so that the product —-r^ changes sign at a certain moment (due 



to the change of sign of A^), keeping the absolute value approximately constant. If 



D depended mainly on the second term of eq.(B.13), to explain the experiments the 
size of the tunneling barrier would have first to increase with P and then to decrease, 
which seems more unlikely. 
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B.2 Path integral description of the dissipative two- 
state system. 

We will here remind the reader about the fundamental concepts underlying the 



derivation of some of the results stated in section 2.3.2 For a more detailed presenta- 
tion of the basics see [FV63, Sim05, PIW], and for a detailed study of the spin-boson 
problem using a path integral approach see [LCD+87, Wei99]. 

A formulation of Quantum Mechanics different from the canonical one developed 
by Schrodinger and Heisenberg but equivalent to it was developed by Richard Feyn- 
man in the fourties [Fey48], the Path Integral formulation, based on the following 
postulates [PIW]: 

1. The probability for any fundamental event is given by the square modulus of a 
complex amplitude. 

2. The amplitude for some event is given by adding together all the histories which 
include that event. 

3. The amplitude a certain history contributes is proportional to e*"^/^, where S is 
the action of that history, or time integral of the Lagrangian. 

This approach is often used in semiclassical physics, where the transition from classi- 
cal to quantum behavior of a system is studied. In terms of path integrals the classical 
behavior is obtained imposing h ^ 0, because in this way the interferences between 
the contributions e*'^^^^*'"''^ of neighboring paths x{t) will cancel out except for those 
close to the path for which 6S/6x(t) = (imposing of course the proper boundary 
conditions at the initial and final time), which is by definition the classical trajectory. 
The first quantum corrections will be obtained by considering also the contributions 
of the paths close to the classical solution, something usually done using stationary 
phase approximation calculations [Sim05], where the otherwise unmanageable inte- 
grals reduce to Gaussian integrals. 

Feynman's ideas constitute also a very convenient framework to study the dy- 
namics of open quantum systems, where a subsystem focus of our interest is not 
isolated but coupled to many external degrees of freedom (its "bath" or "environ- 
ment") [FH65]. In the path integral approach, one starts with a generic path integral 
containing all the possible histories followed by all the variables of the composite 
"subsystem -|- environment" system. As we are only interested in the variables de- 
scribing the subsystem, the second step is to integrate out the rest of variables as best 
as we can, analogously to what is done in studies involving reduced density matrices. 
Typically, the resulting integral, now containing only the variables of the subsystem, 
contains an action S composed of two parts, the action ^isoi of the variables of the 
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b) 

























^1 1 



Figure B.l: a) Particle in a double well potential. The influence of purely quantum 
mechanical processes such as tunneling is included in the path integral semiclassical 
approach by considering the time as a general complex variable r, and searching 
for stationary solutions 6S/6x{t) = whose contributions to the path integral are 
summed; after the solution in terms of r is found, the real-time dynamics is obtained 
by performing an analytical continuation to real time. In a given history of the particle 
dynamics, each tunneling process is called " instanton" , and in b) a prototypical multi- 
instanton stationary solution is represented. Image (b) from [Sim05]. 



subsystem in absence of other degrees of freedom, plus a second term containing 
the information about the effects due to the bath, the so-called influence-functional 
'S'infl [FV63, FH65], where memory effects mediated by the bath are manifested in 
terms of "interactions" between the subsystem variables at different times. Fortu- 
nately a broad class of open systems can be considered to be coupled very weakly to 
each of the numerous degrees of freedom of their environment, allowing for a descrip- 
tion of these environmental variables in terms of harmonic oscillators. The related 
path integrals are Gaussian and can be performed exactly [FV63]. Thus the approx- 
imations will start only when dealing with S'infl, not before. 



The spin-boson model or dissipative two-state system [LCD+87, Wei99], eq.(2.6), 
is perhaps the most successful and studied example using the previous ideas, due to 
its generality and applicability to many different systems, as well as due to its star 
role in the study of environment-induced phase transitions. 

In absence of dissipation, the path integral method gives a good solution for the 



dynamics of the "tunneling particle", flg.(B.l5,), by adding up the contributions of 



all the stationary solutions SS/6x{t) = where r is taken as a general complex 
variable, and performing in the flnal result an analytical continuation to real times. 
This allows to take into account quantum paths classically forbidden, including one or 
several tunneling processes. In the jargon of the fleld each tunneling process is called 
an "instanton" or "blip", and the solution is the result of summing all possible multi- 



instanton conflgurations, like the one in flg.(B.l d). The basic approximation involved 
in the calculation is to consider the tunneling processes as completely independent, 
with the total action of a multi-tunneling path with tunneling events S'n being 
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times the action of a single-tunneling path, 5'n = x 5*1, or expressed in terms of 
bhps, there is no interaction among them. 

When a dissipative environment of harmonic oscillators is coupled to this system, 
the path integral resulting from integrating out the oscillators' degrees of freedom 
introduces, as previously told, an influence-functional S'infi[|x(r) — a;(r')|] that pro- 
vides memory to the dynamics. Now the result can be again expressed in terms of 
sums over multi-blip paths, but interactions play a role and complicate enormously 
the calculation. Leggett et al. introduced the so-called non-interacting blip approx- 
imation (NIBA), which basically limits the memory effects to the inclusion of the 
contribution to the action 5* corresponding to interaction between two consecutive 
periods separated by a tunneling process. 

The success of this approximation lies in its numerous virtues, two out of which 
we highlight here for our purposes: i) It is a controlled approximation, with a clearly 
defined range of validity (basically low temperatures and weak tunneling amplitude 
as compared to the characteristic frequency of the ground states of the wells) and 
well defined expansion parameters, ii) It turns out to be a good description in a very 
broad range of cases, from baths leading to a weak damping to baths leading to a 
suppression of tunneling and localization of the particle. 



Remember that the Hilbert space of the system described in fig.(B.l5,) is reduced 
in the spin-boson model to the two ground states of the wells, represented by a 
spin up and spin down ((o"^) = 1 for the left well ground state and (a^) = —1 for 
the right well ground state, for example). These two states become mixed by the 
tunneling amplitude Ag, introducing spin-flip processes in the dynamics. The basic 
result obtained within NIBA for the time evolution P{t) = {az{t)) starting at t = 
in a given well, {az(t = 0)) = 1, is 



Pit) 



2m 



e^'[\ + f{\)]-'d\ 



c 



where C is the standard Bromwich contour, and 



exp — 



At + Q2{t) 



/(A) = Al I cos Qi(t) 
Jo 

with Qi(t) defined by 

Qiit) ^ 

Jo ^ 
f°° J(uj) 

Q2{t) = / ^^{l-cos[ujt])coih[hwl2kBT]duj 
Jo ^ 



(B.15) 



(B.16) 



J(uj) . , , , 
— — smlutldu 

1'^ 



(B.17) 



The function J{uj) appearing in the integrals is the spectral density of the bath. 



J(^) 



dte"^'{H,,r.tit)H,ntm 



(B.18) 
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At T = 0, up to second order in perturbation theory, it can be rewritten as 

J{uj) = Y,\M''K^-^k) (B.19) 

k 

J{uj) reflects the evolution as a function of u of the amount of possible transitions 
(fluctuations) from a starting initial state (at T = the ground state) to all flnal 
states for whom an energy huj has been exchanged between our subsystem and the 
oscillators' bath. 

Many of the oscillators' baths found in nature, like phonons in solids or e-h ex- 
citations in metals, display a spectral density J{uj) ~ cu'^ up to a certain cutoff (e.g. 
the Debye frequency in the case of phonons). A very important result states that 
whenever NIBA is applicable, the power s of J{uj) is the crucial factor determining 
the effect that the environment exerts on the TLS (((^^(t))) dynamics. As explained 



with an example in Appendices C.2.1 and C.2.3, at flnite temperatures T ^ to the 
amount of accessible transitions thanks to the thermal energy present in the system 
determines J{uj,T > 0), leading to J{uj,T) ~ T*. 



B.2.1 Derivation of the spectral function for the case of the 
modes of a quasi- ID nanoresonator 



The starting point is the hamiltonian (2.10), where diUj is a component of the 



deformation gradient matrix. In the case of the bending modes of a rod of dimensions 
L, t and w, and mass density p, there are two variables X(z),Y(z) obeying eqs.( |A.14 ). 



One can thus express X{z), Y{z) in terms of bosonic operators, for example X{0) 
J2k V V {2ptwLuJk) (flfc + cik)- We can relate this variables to the strain fleld diUj 
through the free energy F: 

Frod = - J dzEIy(^—j +EI,(^—j = -J dz J dS-Xj^ul + f^J^^'^ik 

i i,k 



2 



^ IJdzj dS{^X + 9fj.)ul (B.20) 

extracting an average equivalence for one component Uij, Uij ~ 2 -y/ Ely/ (3 A + 18fi)twd'^X/ d 
The interaction term in the hamiltonian is then 

Hint = hcr^y^ X^={al + ttk) 

k * 



z\ 
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So we have approximately 9 times the same hamiltonian, once for each Uij, and the 
corresponding spectral function J{uj) will be nine times the one calculated for 



int 



^2 

ha^y^ \—={al + ak) 

k ^ 



EL 



e y (3A + l^ij)tw V 2ptwhL 
For this hamiltonian J{uj) is given by 



;4 + a,) (B.22) 



EL 



1 



(3A + 18fi)tw V 2ptwhL y/uJk 
Taking the continuum limit (-^ Ylk ^ n I '^^)" 
2L Z"''™" r AS 



5{uj - ujk) 



(B.23) 



J{u) 



i2ny 



dk 



27- 



EL 



e V (3^ + 18yu)^w V 2ptwhL 



5{u - uok) (B.24) 



Using the dispersion relation iOj{k) = ^jEIj/ {pah) x k"^ = c x k'^ we express the 
integral in terms of the frequency: 



J{UJ) 



L 



- duJk 



AX 

27^ 



EL 



(27r)2 J^^^^ y/cuj^ I e y (3A + 18p)tw \ 2ptwhL\ Uk 



5{UJ - LOk) 



(27r) 



27^ 



EL 



e y (3A + l%p)tw y 2ptwhL 



^5/2 



(B.25) 



</flex(t^) is just 9 times this, eq.(2.14). 



For the compression and twisting modes the derivation follows analogous steps. 
Longitudinal or compression modes correspond to the displacement vector Ui obeying 



the wave equation (A.13 ). Hence, there are plane waves with linear dispersion relation 



k. The strain field du appears in the hamiltonian coupled to a TLS located 



at z = 0. The quantized u(0) and du \z=o are given by 

du 1^=0 



«(0) = 



h 



2ptwLuJi. 



(4 + flfc) 



h 



2ptwL 



E 



[4 + a,) (B.26) 



The hamiltonian is: 



= ea^ + ha^Y "^^{0.1 + ttk) + Y^k{alak + h.c.) (B.27) 



k 
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With such a form for the couphng the spectral function is given by eq.(2.11). 



The twisting or torsional modes obey the wave equation (A. 15). In analogy with 
the previous case we can quantize and get the same results for 0(0) and d(j) \z=q as 
for u(0) and du \z=o, just substituting p by pi and E by C. But in our hamiltonian 
the coupling term involves du \z=o and not 90 [2=0, so the latter has to be expressed 
in terms of the former. Proceeding as with the bending modes: 



rod 



dzC 



d(j) 



dz 

dz I dS8p 



dz / dSAp 



duj. 
dz 



+ 



dUy 

dz 



dz 



Approximating further 



J dS8p 



dux 
dz 



C 



d<p 
dz 



(B.28) 



(B.29) 



the equivalence du \z=o= dux/dz = ^yC/{8ptw)d(f) \z=o is found. Therefore in this 
case one has an interaction term in the hamiltonian 



int 



^0 



C 



k 



e y 8ptw V 2pIhL Ju^k 



;4 + afc) 



(B.30) 



leading to an ohmic spectral function, eq.(2.11), as in the previous case. 



B.3 Dissipation from symmetric non-resonant TLSs 

B.3.1 Spectral function of a single TLS coupled to the sub- 
ohmic bending modes 

We follow the method of ref. [Gui85]. The form of A(uj\ the spectral function of 
a single TLS, for frequencies oj and uj ^ e can be estimated using perturbation 
theory. Without the interaction, the ground state |a) of the TLS is the antisymmetric 
combination of the ground states of the two wells, and the excited state is the sym- 
metric one, \s). We will use Fermi's Golden Rule applied to the subohmic spin-boson 

hamiltonian, H = ea^ + hXaz J2k ('^I + '^fc) + Xlfc fy^{k)a\ak, where is the 
annihilation operator of a bending mode k. Considering only the low energy modes 
uj{k) <^ e, to first order the ground state and a state with energy uj{k) are given by 

1^) - k) - ^4ls) + ••• > 

\k) ^ al\a) - ^\s) + ... . (B.31) 
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We estimate the behavior of A{uj) by taking the matrix element of (t^ between these 
two states, obtaining (remember that u!{k) oc k"^) 

A{uj) + ... , u{k) < e. B.32 

The expression in the numerator is proportional to the spectral function of the cou- 
pling, J{iS) — a^-^ujcosf^- Now we turn our attention to the case ^ £, where 
the ground state \g) and an excited state |A;) can be written as 



\k^l^k ti \ , 



l*> - + ^^^W + ■.. ■ (B.33) 



— Is 

The matrix element (Olcr^lA;) is ~ (fi^I)^ ' leading to 

A{u^) - ""'^f + ... , (B.34) 

B.3.2 Value of A'^^_^JuJo) 

Now we will add the contributions of all the non-resonant TLSs using the probabil- 
ity distribution P(A^, A§) = Po/Ag [AHV72, Phi72]. For the case of weak coupling, 
Qift < 1/2 and ujq > {2abf'uJco-i which is the one found in experiments, one has 

J hluJo+r{a;o)] J-A^ ^0 \^o) 

n[u>o-r{oJo)] /-Ag p (y^.,/7J-(^f:)2 



obtaining the result A*g_j.gg(a;o) ~ 2Pai,-\/uJco/^o- 



B.3. 3 The off-resonance contribution for T > 

Using the same scheme, the modifications due to the temperature will appear in 
the density of states of absorption and emission of energy corresponding to a "dressed" 
TLS, A{e,uJ^T). Now there will be a probability for the TLS to be initially in the 
excited symmetric state, proportional to exp[— e/ZcT], and to emit energy hw 
giving it to our externally excited mode, thus compensating the absorption of 

energy corresponding to the opposite case (transition from \a)\k,n) to n — 1)), 
but contributing in an additive manner to the total amount of fluctuations, which 
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are the ones defining the linewidth of the vibrational mode observed in experiments, 
fixing the value of Q~^{uj,T). The expression for A{e,uj,T) is given by 

A{u;,T) = ^5^$^|(zk.|/)re-^5[no;- (B.35) 

i f 

We consider a generic state l^a) = \a)\kini, kjUj, ...) or \is) = \s)\kini, ...,kjnj, ...), 
and states that differ from it in huj, 

\fa±) = \a)\kini,...,kjnj±l,...) , \fs±)\s)\kini,...,kjnj±l,...) (B.36) 

As for T = 0, we will correct them to first order in the interaction hamiltonian 
Hint = hXaz v^^l^^I + cik), and then calculate the square of the matrix element 
ofa„ \(i\az\f)\\ 

Elements \{ia,s\o'z\fa+,s+)\'^ correspond to absorption by the "dressed" TLS of an 
energy hjjj from the mode kj, while elements \ {ia^s\o'z\fa-,s-)\'^ correspond to emission 
and "feeding" of the mode with a phonon fihjj. The expressions for the initial states 
are 

\ia) = \a)\kini,...,kjnj,...) \a)\kini, ....kjUj, ...) + ^^^^^^^^^\s)\...kini - I...) 

\is) = \s)\kini, ...,kjnj, ...) \s)\kini, kjUj, ...) + ^ ^^j^-^^\a)\...kini - 1...) 



- + (B.37) 

Vfc, 

and, for example, the state \fa+) is given by 
|/a+) = \a)\kini,...,kjnj + l,...)^ 

sr^ fi'X\/'niUJi 

|a)|A;ini, fcjnj + 1, ...) + - 2e I'^^l""^^^ ~ + I---/ 



/itUj + 26 



|s)|...fcjnj +2...) (B.38) 



with similar expression for the rest of states. The value of \{ia,sWz\fa+,s+)\ (absorp- 

tion) is \{^a,Mz\fa^,s^)\ - VK±i)^ 



(?!a)j)2-4e2 



with rij = 0, 1, while for emission the 
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result is \{ia,s\crz\fa-,s-)\ = {h!^Z^ ' ^^^^ '^j ^ 1' ■■■■ Taking the limits we con- 
sidered at T = {hcuj <S s and Huj ^ e) the results are the same as for T = 
for absorption, except for a factor {nj + 1), and wc also have now the possibility of 
emission, with the same matrix element but with the factor n^: 

¥1772 , nwj > e 

Emission ^ ^ (B.39) 

Now we have to sum over all initial states, and noting that the first order correction 
to the energy of any eigenstate is 0, the partition function, Z, is easy to calculate, 
everything factorizes, and the result is, for example in the case hcuj <^ e: 



Z ^ ^ £2 

iy^j^ni=0 nj=0 



kT 



3 



1 — e fcT 



^ OO OO f- I I 



i^j,ni=0 nj=l 
£2 

1 _ e fcT 



fcT 



(B.40) 



In y4abs(£, cUj, T) we have added the contributions from the matrix elements 

|(i„|<7,|/„+)|2eW + |(i,|a,|/,+)|2e^ = |(i„|<7,|/„+)ne^ +e^] (B.41) 

The sum of exponentials cancels with the partition function of the TLS (appearing as 

a factor in the total Z), leading in this way to the expression above (the same applies 
for Aem,(^, — 7")). The total fiuctuations will be proportional to their sum, which 
thus turns out to be at this level of approximation oc cotanh[/ia;j/kT]: 

Ad.iss{e, (jj, T) = Aahs{e, ujj, T) + A^rn{e, -ujj, T) = Adiss{e, ujj, T = O)cotanh ^ 

(B.42) 

In fact, this result applies for any other type of modes, independently of its dispersion 
relation, provided the coupling hamiltonian is linear in and (al^ + at) and every- 
thing is treated at this level of perturbation theory. Moreover, it can be proven that if 
one has an externally excited mode with an average population {rij) , and fluctuations 
around that value are thermal-like, with a probability oc exp[— jn^ — {nj)\/kT)\, one 
recovers again the same temperature dependence, oc cotanh[/ki;j/kT]. 
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B.4 due to the delayed response of biased TLSs 

(relaxation mechanism) 

The derivation presented here is based mainly on pages 396-406 of ref. [Esq98]. 
We want to deduce the formula for the relaxation contribution to the attenuation of 
acoustic waves: 

Consider a set of TLSs with a difference of energy between their two eigenstates 
e = 82 — El and with instantaneous occupation fractions rii and n2, that obey 

^ = ^^ + ^^ (B.44) 
rii = -z/2ini + z/i2n2 

where i>i2 and 1^21 are the transition rates due to the interaction of the TLS with the 
phonons, which reflect the finite response time of the TLS to a perturbation. This 
delayed response is responsible, together with the finite bias Aq > 0, of the dissipation 



(B.43), as we shall demonstrate. We will simplify the elasticity tensors and treat them 
as scalars. If we apply a periodic stress a(t) = aoc"*'^* we perturb the energy levels, 
in the linear response approximation, according to 

Eaicr) = EaiO) - VoXfjaij , = £^(0) - VoPijUij (B.45) 

where 

• f is the unit cell volume, 

• Xfj = —{l/vo)dea/daij is the elastic dipole of state a (also defined as the strain 
Uij created by a homogeneous distribution of c defects per mole, A°- = duij/ dc"), 

• P?j = —{^/vo)dea/duij is the so-called double force tensor (defined too as the 
stress (jjj created by a homogeneous distribution of c defects per mole, p^- = 

dai,/dc") 

Note that from the first definition of one can see that it is approximately propor- 
tional to the deformation potential 7^" , 7^" ~ {l/2)dA'^/duij = —vopfj . This last 
quantity is related to our coupling constant 'jij ~ 7 of the interaction hamiltonian 
i^int through jij = 'iipffj — ■yjj). So returning again to the applied periodic stress, 
it perturbs the energy levels, causing a variation in the TLS populations and the 
corresponding strain, according to ujj^^ = CaXfj ■ We want to find an expression 
relating Aujj^^ and cr{t), from which we can deduce the associated dynamic com- 
pliance A5'(ti;) = Au/a. S{u!) is the inverse of the relative change of the elastic 
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stiffness, SC(ij) = a/Au, the magnitude giving us the frequency shift and quahty 

factor (eq.( B.43"| )) due to relaxational processes of TLSs. 

We start dividing the instantaneous occupation fraction ni{t) as 



ni{t) =n'l + Ani(t) + Sni{t), 



(B.46) 



where is the equihbrium value in the absence of stress, Ani(t) the perturbation to 
the equihbrium value due to the applied stress. 



Ani(t) 



de da 



Anie" 



-iuit 



(B.47) 



ni(t) = n° + Ani(t) is the instantaneous equilibrium value, 6ni(t) is the deviation 
from the instantaneous equilibrium value due to the delayed response of the defects, 
and Ani(t) = Ani(t) + 6ni{t) is the deviation from the static equilibrium in the 
absence of applied stress, 

Ani(t) = Ariie"'^*. (B.48) 

The definitions of Ani(t) and 6ni{t) imply that their amplitudes are complex, because 
they include a phase lag between the excitation and the response, due to the finite 
time needed for reaching the instantaneous equilibrium according to (B.44). The 
strain Aufj^^ we need to compute AS^u) is that due to Ani, which is 



Am = AiA^i + X2An2 



c de 
f da 



Ani 



(B.49) 



Now we want to extract the value of the complex phase appearing in eq.(B.48), and 
for that sake we have to manipulate a bit the rate equation (B.44): we begin using the 
condition for the equilibrium values of the populations (detailed balance condition), 

dni{t) 



dt 







ni{t)u2i = n2{t)ui2 



With this relation we transform ( B.44[ ) into 

duiit) 6ni(t) 



dt 



T 



(B.50) 



(B.51) 



where r ^ = 1^12 + 1^21- (B.51) indicates that the rate of change of the population is 
proportional to the deviation from the instantaneous equilibrium. Noting that 

dni{t) dAni{t) 



dt 



dt 



-iujAriiit) 



and 



5ni{t) = Ani(t) - Ani{t) 



(B.52) 



(B.53) 
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and using (|B.51|), one obtains 

UJT 



^^1 

Aril = = Aril 

1 — iujT 



+ 1, 



Ll + (cur)2 l + (tur)2j 
So the dynamic compliance can be written as 



(B.54) 



ASH = ^ = ^5|l(|i)'-^ (B.55) 

When one compares this result to the static compliance, it appears that the dynamic 
response function is obtained multiplying the static one by (1 — iojt)'^ . With the 
dynamic modulus 5C{uj) an equivalent expression can be derived, with the roles of u 
and a exchanged. 



SCijhki^) = ^ = — V j^j. — z : — B.56 

a/3 ^ ■' 



Now, this can be rewritten observing that 

drtn 



leading to 



(3na{ni3 — dajs) and ^^nQ, = l (B.57) 



SQM^) = -f3n^nJ^ - - j^)-^ (B.58) 



To get to (B.43) we need a couple more transformations of this expression. First, if 



the energies of the eigenstates under strain u are given by 



then eq.(2.21) is obtained, so that, together with the definition e = 62 — £1, and the 
equivalence nin2 = (l/4)sech^(£/2kBT), results in 

SQjhkiuj) = -(—)\ij^hkjsech\e/2kBT:)-^— (B.60) 

This is the result for a concentration of c defects per mole of TLSs with fixed e and 
A^. In the case of an amorphous solid one has instead a probability distribution in 
terms of the parameters A^. and A^, g{Arc, Az)dAxdAz = Po/A^dAxdAz, which can 
be expressed with the variables e = a/ {AxY + (A^)^ and u = A^/e as g{e,u)dedu = 
Po/{u^/l — v?)dedu , with e running from to emax and u from Umin to 1. Writing 
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ir(A) 


CO 





A/co 



Figure B.2: Evolution with Aq of the different quantities determining the approxi- 
mations to be taken in the integrand loqt j\V + (cjo^)^] of eq.(B.43). 



SCijhki^^) in terms of e and u, the inverse quahty factor associated to the variations 
of the Young Modulus E, the component of Cijhk entering the equations of motion of 
the vibration, is, as anticipated. 



Q~\uj,T) = Im((5E)/E 



pv^ 



de 



du 



u^ 



1 



UJT 



u kTcosh2[e/2kT] 1 + (cur) 



(B.61) 



B.5 Derivation of Q-f(a;o,r), eq.(2.23) 



As discussed after eq.( B.43[ ), we have to sum over underdamped TLSs, e > 
[SOab^cUcoT]^/^, using the approximation for F, T{e,T) ~ 30ab^/uJ^T/ ^/e. Moreover, 
if (jJq > F(e = [SOaby/uJ^T]"^/^ , T) then in the whole integration range ujq ^ F(e, T) ^ 
ujQT{e, T) > 1, so that uqt/[1 + (o^ot)^] ^ l/(a;or), see fig.(|R2|). Q^ci (^o, T) follows: 



gJ(wo,r) 



de 

[iOab^/lJ^T]^/^ Ju. 



du 



^/l^^T{e,T) 



u 



(B.62) 



For temperatures T ^ [SOab^cUco^]^^^, which holds for reasonable T and sizes, the 
integral, which renders a result of the kind Q~^{uoq,T) ^ {1/ujo){A^/T — BT^^^), 
can be approximated by just the first term, obtaining eq.(2.23). In any case for 
completeness we give the expression for B: 



B 



500Po7''/^ (1 + i^)^/^(l - 2z/)^/3 / P \ 1/3 



Ell/3(3_5;y) 



E 



(B.63) 



Also for completeness we give the result for higher temperatures, although for current 
sizes the condition F(£: = [30ab^/uJ^T]'^/^ , T) > luq implies values of T above the range 
of applicability of the Standard Tunneling Model. Now for some range of energies 
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r(e,T) > cjq; and the range of integration is divided into two regions, one where 
ujqT ^ 1 and one where the opposite holds: 



Qr^M.T) = %f du^^^^x (B.64) 



u- 



de LJoT{e,T) + [ de 

Jllfiaty/UJ^2T/^]2 i^OT[e, I ) J 



The final result is Q~^\{ijJq,T) ^ —IPq'^'^loq/T + A\/T/loq — CT/uJq, with A defined 



by eq.(2.23) and C by: 



1500Po7M1 + z/)^(1-2z/)V pV/2 

All the results for Q~^{ijJq^T) have to be multiplied by the fraction of volume of the 
resonator presenting amorphous features, Vamorph/ (twL). 



B.6 Derivation of eq.(|2^ 

To calculate classically the energy stored of a torsional mode t) = A sin[(2j — 
1)ttz/ (2L)] sm{ujjt) we just calculate the kinetic energy in a moment where the elastic 
energy is zero, for example at time t = 0. If an element of mass is originally at position 
{x,y,z) {x,y transversal coordinates), with a torsion (t){z^t) it moves to 



r(t) = ( a/o;^ + ■ cos 



X 



arccos 



+< 



y'^ ■ sin 



arccos 



B.66) 



The kinetic energy at time t = is 



dz 



t/2 



dx 



-t/2 



w/2 ^ 

dy ■ 

-w/2 ^ 



dr 



(B.67) 



Substituting the expression for f{t) in the integrand, one arrives at 



dz 



t/2 



dx 



t/2 



w/2 



w/2 



dy--pA sin 



(2j-l)7r 
2L 



{x' + y') 



= —A^u']pL{t^w + wH) 
48 ^ 

In terms of the creation and annihilation operators 



(t)j{z,t) 



h 



i{kjZ—u)t) 



(B.68) 



(B.69) 



2LpIujj 

so the mean square of its amplitude is (0^) = A^ /2 = h(2n+l) /[2LpIu!j]. Substituting 
this in eq.(B.68) the eq.(2.27) for Eq is obtained. 
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C.l Elect ron-phonon coupling: friction in terms 
of the susceptibility % 

C.1.1 Damping of a phonon mode due to Coulomb interac- 
tions between charges in a device 



We have seen in subsection 6.4.2 how the response of system to a given external 
perturbing field can be expressed, for weak perturbations, in terms of a susceptibility 
or response function. To be more concrete, we started from an interaction hamiltonian 
Hint acting as a perturbation on the initial one, Hq, and its eigenfunctions (pi, 

Hint = J drV{r,t)p{r) (C.l) 
where p(r) = X]kii'[c^^j^,Cke*^''"+h.c.] is the charge operator, and the potential is an 



oscillating function of frequency uj, V{r,t) = F(f)e*^*. We saw in 6.4.2 that to linear 
order in Hint the variation in the particle density n(f , t) = J2T'^ Wiijy with respect 
to the case in absence of perturbation can be expressed as 



5n{r) = j dYx{i,i\oj)V{Y). (C.2) 



This can be taken as the definition of the retarded density- density response function 
X for the charge density change at r induced by a potential field at r'. In terms of a 
given approximation for the single-particle eigenstates \i) and eigenenergies Cj of the 
total hamiltonian Htot = Hq + Hint, it is expressed as 



occ ^ ^ 
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The response frequencies of the system are given by the poles of the response function, 
which depend in turn on the approximation of the basis. For convenience we will 



generalize eq.(C.3) in terms of the A^-body states |a) of our system composed of N 



fermions. |a) can be given, for example, by a Slater determinant of one-body states 
\i). We assume that the system starts in its ground state |0) (T=0 case). We define 
the fermionic creation operators c] creating one-body states \i), so that 



5^(a|4c,|0)(0|cjQ|a) 



^0 



+ 



hhJ 



(C.4) 



Now we are interested in the probability per unit time that the field Vij^t) trans- 
fers energy hw to the system. According to Fermi's Golden Rule of second-order 
perturbation theory the probability we are looking for is given by 



27r 



(C.5) 



Heref = E„(j|^(f,i) 



i)(^iCj is the operator corresponding to V{Y,t) in the second- 



quantized basis c\. Eq.(C.5) can be expressed in terms of x by providing the latter 



with an imaginary part through the addition of a small term it] to and then using 
the relation Im[l/(x — irf)] = 7t6{x). In this way we obtain the expression for F used 
in the thesis, 



T{V,u;) 



drdr'V {r)V {r')lmx{r , uj) 



(C.6) 



Now the goal is to use a good approximation for the susceptibility function of our 
system. 



C.1.2 Clean metal susceptibility 

Suppose that our system is a metal. When the presence of impurities, electron- 
phonon collisions and electron-electron interactions is such that the collision frequency 
z/ = 1/r of an electron is much less than the frequencies u we are interested in, ut <C 1, 
and/or the mean free path is much longer than the length scales L we study, vpT ^ L, 
we are in the clean metal limit [NP99], where the susceptibility can be expressed in 
terms of elementary excitations. If one neglects electron-electron interactions (free 
electron gas), the retarded susceptibility is given by the bare bubble diagram: 

X^{q,uj) = -/S ^. , C.7 

K,<T 

where Fis the volume of the system, np^e) is the Fermi function, — ej^_^- are the 
excitation energies of the quasiparticles and t] —>■ 0. The response of the system is 
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determined by the excitation of these one-body quasiparticle states, corresponding to 
zeroes of the real part of the denominator. The main effect of the inclusion of Coulomb 
interactions among electrons is the appearance of collective excitations, formed by a 
superposition of one-body states, called plasmons in the case of the electron gas. 
Their presence is manifested for example in an RPA calculation. RPA is the result 
of summing, for the interacting electron gas, the most important diagrams of each 
order [BF04], which turn out to be the ones of the bubbles, and these diagrams 
increase their importance as the density increases, so RPA is the high density limit 
of the interacting electron gas. The RPA result is the following: 

i - vcoui{<i)xo iq, 

where Vcoui represents the Coulomb potential. The response will show peaks again 
in the zeroes of the real part of the denominator of x^(q, cu), but now also when 
1 — VcomKq)Xo^(Q) = 0, which corresponds to collective excitations. 



C.1.3 Dirty metal susceptibility 

When due to the presence of impurities, electron-phonon collisions and electron- 
electron interactions the collision frequency u = 1/r of an electron is higher than the 
frequencies u we are interested in, ujt ^ 1, and the mean free path is much shorter 
than the length scales L we study, vpT -C L, one enters the so-called hydrodynamic 
regime [NP99]. A given quasiparticle is subject to many collisions during one period 
of the exciting external field. It is thus completely " thermalized" well before it can 
sample the periodicity of the field. Therefore the best way to proceed is to assume 
that the collisions act to bring about everywhere a state of local thermodynamic 
equilibrium, and use just macroscopic "local" quantities (density, current, pressure, 
etc) to describe the state of the system; the response to the external field may be 
obtained by using the usual laws of thermodynamics and hydrodynamics. For high 
enough temperatures one arrives at an expression for the response function x(q, cj) 
which contains a peak corresponding to the plasmon and a continuous background 



arising from thermal diffusion. This background is the one described by eq.(4.15). 



Taking into account the plasmon contribution, eq.(4.15) becomes 



xq,w =77^ 2 C-9 

Uq'^ — uj^T — luj 

The plasmon contribution to the dissipation of the vibrational modes we study in 
this thesis can be safely neglected. Basically one has to see if Dq^ ^ uP't in each 
of the layers (graphene and inversion layer in the Si-Si02 interface). In both cases 
one can take Dv ~ 10'^, r = Ijvp ~ 10^^/10^ ~ 10"^'^ s (assuming I ~ 100 nm), 
g ^ 1/L ~ 10^ m"\ D = vpl ~ lO^^mVs. The density of states (DOS) of a stack 
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of N graphene layers is i/*^ = N'y/h'^Vp, being 7 ~ 0.3 eV the interlayer hopping 
element. The DOS of the approximately two-dimensional electron gas lying in the 
inversion layer Si-Si02 is i'2Deg ^ ttl * /nh?, with m* ~ O.Sm^- the effective mass 
of the electrons in Si. Using these values, N ~ 10, and for the frequencies of the 
oscillators analyzed, u ~ 100 MHz, one immediately verifies that indeed 

Dq'' > u, uj\ ^ Imx(q,^) ~ ^ (C.IO) 

The inequalities are quite large, so there is place for substantial variations oil 01 N 
and they will still hold, justifying our use of eq.(4.15) and the limit Dq^ ^ u, uP't. 



C.1.4 Microscopic derivation of eq.(4.15) 



We sketch the analysis presented in [AS06]. This form of the susceptibility can 
be derived assuming a distribution of short-range scattering centers represented by 
potentials, creating an V^(f) = — fi)- To calculate propagators and phys- 

ical magnitudes one has to perform in a final step a disorder average {...)dis = 
f DV P\V]{...), where the probability measure P\V] describes the statistical prop- 
erties of the potential V . In most applications it is sufficient to implement a Gaus- 
sian distribution, P\V] = exp J d'^rd'^r'V {y)K^^ {r — v')V{v') , where 7 mea- 
sures the strength of the potential and K describes its spatial correlation profile, 
(y {y)V {r')) = '~f'^K{r — r'). Very often one finds that the finite spatial correlation 
of V does not matter, in which case one may set K{r) = 6{r), so that 



(...) 



dis 



DV exp 



272 



d\V^ 



iC.ll] 



For length scales longer than the mean free path, so that the electrons experience 
multiple scattering, details of the form chosen for the potential V"(r) are quickly 
erased. 

In practice, this average is very difficult to perform, and a series of techniques 
have been developed for mapping it into a different form. For example, in the replica 
field theory, the expectation value of an operator O is written in terms of a source J 
introduced conveniently in the action of the path integral representing the replicated 
partition function, Z^[J] = J D{^, ^) exp [- '£a=i "^[^"^ ^^ J] 



{0) = 



R 



1 B7^\ T] I 

-(d/dJ)lnZ[J] = lim = lim -y(0(^^^'^))^ (C.12) 

^ ' R^o R dJ R^oR"^ 



a=l 



When performing now the disorder average of the functional Z^[J] over the distribu- 
tion -P[l^], the result is the presence of a new effective interaction term in the action 
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of the system, with respect to the case without disorder: 

D(*,fr)exp[-^5|v=o[*",^"]+E]YE / d'r^li^mi^mrXi^) 

a=l a=l mn 

(C.13) 

Then, one can apply perturbation theory to this interaction term to know what 
are the disorder averaged values of observables, or in our case of the susceptibility, 
using diagrammatic theory and an ansatz for the imaginary part of the self-energy 
correction to the fermion propagator, ImS(a;„) = — sgn(a;„)/2r, where r is the so- 
called elastic scattering time, related to the elastic mean free path I by r = l/vp- 
One can see that this parameter is related to the potential strength 7^ of the disorder 
potential by = 27rz/7^, where v is the density of states. Usually 7^ is expressed 
through r from the outset. The main results are 



(G(x,y;T)>ii, = G(x,yiT)|K=oe-l"-''l''"' 
{X(r.r))., = (C.14) 



where x(r, r) is the density-density correlation function (the susceptibility of a dirty 
metal) . 



C.2 Imx, Q temperature dependence of fric- 

tion 

In order to make a connection between theory and experimental data a first point 
that has to be kept in mind is the following: the linewidth observed in the 

frequency domain of the response to an external driving force, does not correspond 
exclusively to an irreversible loss of energy of the mode (damping), but to all the 
processes which alter the phonon population of an excited mode of a certain frequency, 
including both increases and decreases of population. That is, Imx, correspond 
to the fluctuations of the given mode due to its coupling to other degrees of freedom 
of the system, giving a finite lifetime of the mode (ko,nko) subject of study due to 
any kind of process, including those increasing n\f^. Friction, instead, corresponds 
just to processes involving a decrease of njcQ. 



C.2.1 Temperature dependence of Im^, Q ^ due to excitation 
and relaxation of e-h pairs 

Friction experienced by, for example, traveling charged particles through a host 
metal due to the excitation of e-h pairs has been subject of intensive research since a 
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long time [Gui84] . In the graphene resonator's case a similar effect produces absorp- 
tion and ffuctuations of energy, but there is an important difference with the previous 
example: we are studying the linewidth of a single bending mode of a well-defined 
wavevector q and frequency cjq, so we will not sum over frequencies to calculate the 
total absorption and fluctuations of the resonator's vibrational modes due to this 
mechanism. More specifically, we want to calculate 

T{u^, T) = j d^i j c/3r'ReVscr(?, cUq)ReVscr(f', cUq)Imx[? - u:^, T] (C.15) 

And for this we need Imx[r — r',ci;q,T]. To obtain the dependence with T we can 
recall the fact that, embedded in Imx, what we are looking for is the sum of e- 
h transitions that absorb an energy cUq from the vibrational mode (q, Wq) through 
Hint plus desexcitations that transfer an energy u;q to that mode. We will begin by 
analyzing the T = case, where only the former can take place, so that dissipation and 
fluctuations are indistinguishable. We will see that for a given frequency uj -C Ep/h, 
the e-h transitions available are proportional to u, independently of the dimension of 
the fermionic bath. It is a direct consequence of the exclusion principle through the 
restrictions imposed by the Fermi distribution functions np{E) = 1/[1 + exp[{E — 
Ep)/kBT]] appearing. Taking for example a 3D Fermi gas, with dispersion relation 

E = 22- the probability of exciting an e-h pair of frequency u will be proportional 
to the number of such pairs available, P{uj), 

P{u) = j d^^^j d^k'np{E{k)) [l - np{E{k'))j6(u - :^M_^^Mj (c.l6) 
Expressing the integrals in terms of energies and making the approximation hw -C Ep, 



Piuj) oc 



/ dEyjE{E + nw)np{E)[l - np{E + hw)] 
Jo 

i-Ep 

= / dE^E{E + hu) Ep X hw Lu (C.17) 

J Ep — hui 

In the case we took a fermionic bath of a different dimension, the only difference in 
the calculation would arise in the factor which in the 3D case is E{E + fko), which 
will be in general a function f{E,huj). But as long as this is a smooth function of 
E and hw, for hu <^ Ep we can approximate f{E,hw) ~ f{Ep,0) and arrive again 
at P{uj) ^ UJ. At T = only excitations of e-h pairs from the ground state can 
occur, so the net absorption rate involving excitations of frequency u will be indeed 
proportional to u. 

At finite temperatures there will be absorption and emission processes. First we will 
calculate the number of e-h excitations of a certain frequency u which will cause 
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energy absorption, as in the T — case: 

E{k) - E{k') 



P{u;,T)^ J d"k j d''k'nF{E{k)/kBT) 1 - nF{E{k')/kBT) 5(u; ^ 

(C.18) 

For frequencies and temperatures such that hcu <^ ksT <C Ep one can approximate 

/•oo 

P{uj,T) oc / dEf{E,huj)nF{E/kBT)[l-nF{[E + hw]/kBT)] 
Jo 

poo 

f« f{EF,0) / dEnF{E/kBT)[l-nF{E/kBT)] 
Jo 

poo 

= f{EF,0) X kBT X / dxnF{x)[l-nF{x)]r.T (C.19) 

Jo 

So for excitation frequencies much smaller than the temperature the number of avail- 
able excitations absorbing energy is proportional to T, and roughly independent of cu. 
The number of desexcitations compensating this absorption processes has a similar 
expression: 

POO 

P'{cu,T)oc dEf{E,hw)nFi[E + hw]/kBT)[l-nFiE/kBT)], (C.20) 
Jo 

and is also proportional to T, leading to Imx(T), Q''^{T) ~ T. 

C.2. 2 Temperature dependence of energy loss due to excita- 
tion and relaxation of e-h pairs 

The net amount of energy absorbed by the fermionic environment of e-h pairs of 
frequency a; at a finite temperature, A{u!,T), for hcu <^ /c^T <^ Ef, is proportional 
to 

A{u;,T) oc P{u,T)-P'{uj,T) (X j dEf{E,noj)^nF{E/kBT)-nF{\E + hoj]/kBT)^, 

(C.21) 

and it can be approximated as 

A{uj,T) oc f{EF,Q) j dE^nF{E/kBT)-nF{[E + njj\/kBT)^ 

= f{EF,0) X kBT X / dxriFix) 

J-EwlknT 



The net absorption of energy due to pairs of frequency uj is therefore to a first ap- 
proximation proportional to uj and independent of T. 
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C.2.3 Extension to a generic system + bath 

Suppose that the subsystem subject of study, characterized by the index k, is 
coupled linearly to a bath of excitations, 

H^nt oc (4 + Cii) X ^ \^{b{ + b^), (C.23) 
q 

where ct and fet represent subsystem and bath creation operators, respectively. We 
can calculate at T = 

/oo 
dte''''{H,ntit)Hi^tm ~ , (C.24) 
-oo 



which as we saw in Appendix B.2 reflects the evolution as a function of u, up to second 
order in perturbation theory, of the amount of possible transitions (fluctuations) from 
a starting initial state (at T = the ground state) to all final states for whom an 
energy hw has been exchanged between our subsystem and the bath. As temperature 
grows, for T ^ u, the amount of accessible transitions thanks to the thermal energy 
present in the system will consequently grow as T'^, so that Imx(T), Q^^{T) ~ T'^, at 
least when the perturbation Hint is not too strong. 

C.3 Charge impurity density in Si02 and the Si02- 
Si interface 

From refs. [SzeSl, NB82]. In ref. [SzeSl] they speak of four types of interface 
trapped charges: 

1. Just in the interface, due to the interruption of the periodic lattice struc- 
ture, there is the so-called interface-trapped charge Qu. Densities as low as 
Qit ~ 10^° cm~^ are achieved by low-temperature (450°C) hydrogen annealing. 
These charges are associated to traps, called interface traps, which change oc- 
cupancy with gate bias changes and have energy levels distributed throughout 
the bandgap. Interface traps characteristic of thermal oxidation are donor type 
in the upper half of the bandgap. 

2. Close to the Si-Si02 interface there are also charge centers called oxide fixed 
charge, Qfc, predominantly positive, immobile and with fixed charge when a 
gate voltage is applied. Experiments locate it within 30 A of the interface. 

3. There are positive mobile ionic charges, Qm, related to trace contamination by 
alkali metal ions, causing reliability problems in semiconductor devices operated 
at high temperatures and voltages. Films impervious to mobile ions, such as 
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Si3N4 can prevent mobile ionic charge contamination of the oxide during device 
hfe, but there is no such film in our case. Effective densities Qm ~ 10^° cm~^ 
can be achieved without films. 

4. There are oxide traps associated with defects in Si02, which are usually electri- 
cally neutral, and are charged by introducing electrons and holes into the oxide, 
originating oxide trapped charge. 

5. During the diffusion of oxigen creating the oxide layer, if the process is fast 
enough, a considerable amount of dopant impurities (B, P, As or Sb) may be 
left in the oxide, although their diffusion coefficients in silicon dioxide are very 
small compared to the ones in silicon. In [YOTT79] experiments were made 
where up to 20% of the phosphorous was trapped inside the oxide during the 
oxidation in wet oxygen process at temperatures between 750 and 900°C. 

These charges give rise to a charge density in the oxide layer, p^x- Unfortunately, 
the published literature lacks the information on this quantity for oxides used in 
MEMS applications. What has been extracted from high-frequency capacitance- 
voltage curves of MOS capacitors is the effective interface charge density Qf, which 
corresponds to 



In [FKDB02] they mention that the breakdown field value for Si02 is ~ 1 V/nm. 
This means that in typical graphene experiments, where the potentials are about 100 
V applied through an oxide layer 300 nm thick, the field is 1/3 of the breakdown field, 
so most probably all the potentially chargeable impurities within the oxide will be 
charged. Indeed, in [KHDS02] a field of about 0.03 V/nm was sufficient to polarize 
the Si02 dielectric, and in [RFTJ02] they mention that the bulk oxide traps become 
charged with injection of electrons from the nanotube at gate fields above 0.03 V/nm. 

C.3.1 Estimate of the charge concentration, and comparison 



Prom the experimental Qf values in [AP87], where a MOS setup with a 200 nm 
oxide layer, and p-typc silicon (the most common choice for mass-produced chips) 
with doping Na = 4 ■ 10^^ cm^'^ was the subject of the measurements, the authors 
in [WS04] extract the value of the quantity 




(C.25) 



with numbers given in other references to fit experi- 
ments 




(C.26) 
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where tox is the thickness of the oxide layer, assuming several distributions of charge 
in the oxide. If one chooses the charges to be concentrated in a layer of 20 nm starting 
from the interface and Qu = 0, the result is 

Pox ~ 1.5 ■ 10^^ charges/cm^ (C.27) 

If instead of that one assumes a uniform distribution of charge throughout the oxide, 
the result is 

Pox ~ 3 ■ 10^^ charges/cm^ (C.28) 



C.3.2 Thickness of the charged layer in the doped Si gate 

From [SzeSl]. When one applies a voltage between the graphene and the p-doped 
Si electrode, the first effect at low voltages (~ 1 V) is to populate the intergap 
states associated with the impurities, creating a region of a certain thickness where 
all impurities are charged, so there the charge density is just the doping Na, which for 
a typical low-doped Si wafer is A^^ ~ 10^^ — lO^^cm"^. For that scale of voltages the 
thickness soon reaches a maximum value Wm before the complete depletion regime is 
reached. This value can be calculated as 



W„.^J'J*SI^I^, (C.29) 

where ~ 11.9eo is the dielectric constant of Si, ~ 10^° cm~^ its intrinsic doping 
(everything at room temperature), and q the electron charge. For a low-doped Si 
wafer one has Wm ~ 450 nm. The complete depletion regime is reached when the 
voltages applied are enough to curve the bands close to the Si-Si02 interface until 
the Fermi level reaches the bottom of the conduction band. Due to the high density 
of states of the band the charge density will increase dramatically in a narrow region 
~ 10 nm thick as the voltage is further raised, overwhelming the amount of charge 
corresponding to the intergap states, so for high voltages (for example typical values 
V ~ 100 V) the approximation of a 2D charged Si layer is reasonable for typical oxide 
thicknesses tsi02 ~ 300 nm. 

The voltage that really matters is the one at the Si-Si02 interface, which can be 
safely taken as of the order of the potential applied between the electrodes. 



C.3.3 Information about the Si02 surface and structure. Con- 
clusions from experiments with carbon nanotubes 

The strong adherence of the graphene layer to the Si02 substrate is probably due 
to the presence of Hydroxyl groups in the surface of the latter, which form hydrogen 
bonds with the carbon layer. When setting the graphene sheet into periodic motion. 
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a repeated breaking and healing of some these bonds (hke a Velcro) can be a source of 
short-wavelength phonons which carry away part of the energy initially stored in the 
resonator, thus damping its motion. How many OH's are there? In ref. [SG95] they 
analyze Hydroxyl groups on a well-defined silica surface, and mention that isolated 
Si-OH will survive anneals at over 1000 "C, presumably because there is nowhere for 
the hydrogen to go: it needs to find another hydrogen to be removed as H2 or even 
H2O. As the concentration of H at the surface is reduced this becomes increasingly 
improbable. Adding hydrogen atoms from e.g. a hydrogen plasma helps dcsorb the 
remaining isolated Si-OH, but in any case the surface density of a previously fully hy- 
droxylated surface after high-temperature anneahng is ~ 1 OH/nm^ = lO^^OH/cm^. 
This is in agreement with another experiment [DPX98], which estimates a surface 
density of isolated silanol groups on a fiat silica surface of about 9 ■ lO^^silanol/cm^, 
thus supporting the idea of either a simultaneous breaking of a great number hydro- 
gen bonds, for a resonator ~ l/im wide, or no breaking of bonds at all. 

The presence of oxide trapped charges within the Si02 and of other surface impu- 
rities has attracted the attention in several studies where its infiuence and possible use 
in carbon nanotubc field-effect transistors has been contemplated [FKDB02,RFTJ02, 
KJV+03, FTK+06]. The most interesting conclusion for us is the one of [KJV+03], 
where they show after systematic studies that there is a monolayer or submono- 
layer of water molecules hydrogen bonded to the Si-OH silanol groups all over the 
Si02 surface, which cannot be removed by pumping in vacuum over extended peri- 
ods of time at room temperature. These molecules can store charge when a volt- 
age is applied between the CNT and the underlying p-doped Si backgate, and are 
the main cause of hysteresis in I-V cycles, much more than the underlying oxide 
trapped charges within the Si02 substrate. These bulk charges were invoked incor- 
rectly in [FKDB02, RFTJ02] as the main cause of hysteresis. 

Returning to the surface water molecules, they manage in [KJV"'"03] to estimate 
the charging per H2O to be roughly 0.1 |e| per molecule. 

They also suggest ways to get rid of the molecules by combining i) thermal anneal- 
ing in a dry atmosphere with simultaneous passivation of the surface with PMMA, 
which links strongly to the silanol groups and is also hydrophobic, ii) the permanent 
immersion of the passivatcd surfaces in dry oxygen atmosphere. In [FTK+OG] they 
also point out that a PMMA coating of the oxide surface changes locally the dielectric 
constant (it is 3.3-3.9 at 60 KHz, and 2.2-3.2 at 1 MHz, a bit smaller than the one of 
the oxide, 3.9). 

In [FTK+06] they say that logically the regions in the oxide which will be first 
charged will be those subject to the highest fields, which in the case of nanotubes will 
be close to the nanotubc. They also estimate that trapped charges close to the CNT 
are expected to produce potential fiuctuations on the length scale of the screening 



length Xscr ~ x h^^tcNTtox ~ 40 nm iecNT ~ 20, tcNT ~ 3 nm is the diameter of 
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the CNT and tox = 100 nm is the oxide thickness in their case). This formula for 
estimating the screening length is derived in [YOL92]. 



C.4 Screening of the potentials at the graphene 
sheet and Si gate 

The equations for the selfconsistent potentials Vscr{z, r — r', u) as a function of the 
bare potentials vl{z, f — f', u) are given by 



Vscr{d,r - r',uj) 
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(C.30) 




+ 


/ dri 

Jc 


/ c/f 2f coui(c^, r - ?! 
Ic 






- r2 


iu)Vscr{d,r2 


- y\uj) 




+ 


/ dr-i 
Jg 


/ (ir4Wcoui(c?, r - f 3 
Ig 


^] 
?' 




- r4 


a;)fscr(0,f4 


-~r\uj) 




+ 


/ dri 

Jg 


'G 


^] 




- ^2, 


U;)Wscr(0,f2 






+ 


/ drs 

Jc 


/ (ir4f coui(0, f - rg, 
Ic 






- f4, 


u;)t;5cr(c?,f4 





where for example in the first equation (rf, r — v' ,uj) represents the bare potential 
experienced by a point charge e in the graphene layer due to the presence of charges 
in that same layer, while r — r', cj) is the bare potential experienced by a point 
charge e in the graphene layer due to the presence of charges in the Si plane, f coui is 
the two-dimensional bare Coulomb potential. These equations simplify considerably 
in the q space, assuming 2D translational invariance: 



Vscrid, q, u) = v^{d, q, a;)e'''^ + v^{d, q, uj) + 

+ VqX'^'i^, Uj)Vscr{d, q, u) + V gC' '^'^ x'^ Uj)Vscr{0, q, u) 

< 

Vscr{0, q, = ^'^(O, q, uj) + v^{d, q, uj)e'"^+ 

+ fge"'?'^x'^(q, u)vscr{d, q, u) + Vqx'^iq, uj)vscr{0, q, u) 

(C.31) 

, where Vg = 27re^/|q| is the Fourier transform of the Coulomb potential in two 
dimensions, and where f^(0,q, cj) and q, cj) have been expressed in terms of 

VQ{d,q,ij) and t>^(0,q, u;). Now we can calculate Vscr{d,q,u) and fscr(0, q, tu) in 
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terms of the rest of the variables, 

The dependence on q and uj has been omitted for the sake of clarity. Now, if we 
are interested only in the long wavelenght limit VqX^ I'^qX^ ^ 1) the last equation 
simplifies to 



Vscr{d) 



X (C.33) 




(x^e-^-^ - x^e«<^) Vq (-x° + X^) 
[-X^ + X^) Vq (x^e-'i'' - x^e'i") 

Values of ^;^(0,q, a;) and q,a;) 

Now we will calculate the parts of these terms which will give rise to a coupling 
to the vibration. When the graphene layer is set into motion with a bending mode 
of wavevector q and amplitude Aq, the potential (0, f , t) of a point charge e in the 
Si plane due to the charge in the graphene layer is 

v^{Q,v,t) = ^ / dv'voU^-v\z')p{v\z',t) 

J c 

27re^Po 



dr 



'^Jc V(f - r'y + {d + ^qe'(^^^'-'^^*))2 
1 f 2ne'po ^1 r ^^,2ne'poA^e'(^''--^')d 



2 Jo Vir-fT + d^ 2Ja ^ff - F)^ + d^V^'' 



^ f{i) + ne^PoA^e-'^ie'^^^-''^^^ {CM) 

where in the third line an expansion for small has been performed. The Fourier 
transform for a; 7^ is 

v^(0, k, a;') = ne^poA^e-'^^ik - q)5(a;' - a;q) , |q| = 1/L (C.35) 

Similarly, the potential of a point charge in the oscillating graphene sheet due to the 
charge in the Si plane VQ{d), is 

v^{d,r,t) = l [ dr' , ^ = ^ /(r) + ne' poA^e'^"^'-^'^ 

° ^ ^ 2 7g V(r-f')2 + (d + Aqe^(^f-^*))2 ' ^ ^ 

(C.36) 
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leading to the same expression as eq.(C.35) but without the factor e '^'^ 



Substituting daSSfaSTl ) in eq.(|C.33l), one obtains eq.(|4.17l). 



(C.37) 



C.5 Dissipation due to two- level systems in graphene 
resonators 



C.5.1 Model and vibrating modes of a 2D sheet 

The hamiltonian describing the couphng of the effective TLS's and the oscillating 
graphene sheet is given by [SGC07a] 

H = ea, + 7^a, ^k(&k + &!k) + Yl ^Abu (C.38) 

k k 

where e = ^ (Ag)^ + (Ag)^, fej. represent the phonon creation operators associated to 
the different vibrational modes of a sheet, and Ak(&k + ^!k) represents the coupling 
to the strain tensor Uik- There are two types, compression modes (longitudinal waves) 
and bending modes. The equations governing the former are [LL59] 
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the equations simplify to wave equations: 



E d'^u. 
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(C.39) 



(C.40) 



(9t2 p(l - 1/2) d^2 ' 2p(l + U) 

so that the speed of the longitudinal {vi = \/ Ej p{\ — z/^)) and transversal (ft = 
y/E/2p{\ + v)) waves are different. Therefore we will have two types of compression 
modes, longitudinal and transversal, each parametrized with a vector k in the plane. 
The expression for the local deformations associated to the longitudinal component 
in second quantization is 




h 



ib[ 



bli) 




{bi + b^i) (C.41) 



and similarly for the transversal one. They give rise to an interaction hamiltonian 



int = 

j=x,y i=l,t k 
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(C.42) 
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which has an associated superohmic spectral function 



(C.43) 
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In the case of the bending modes the governing equation is 



ptw 



Et^ 



(C.44) 



(C.45) 



so there is a non-hnear dispersion relation, u; — k"^ y/ Et"^ / 12 p{l — i/2). To obtain the 
approximate relation between the coordinate z which will have the bosonic opera- 
tors associated, and the strain tensor Uik appearing in the interaction part of the 
hamiltonian, we equate two expressions for the elastic free energy. 
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(C.46) 
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(C.47) 



2(1 + 

and obtain an average relation 

(Uij) t^ 

With this information Hint reads 

^ ''%£,J?.3?^'^y3(l - .)(9+ ^) V2nML-p 

(C.48) 

Proceeding as with the previous modes, the spectral function associated turns out to 
be ohmic: 



t 



lm\2 



Im 



e 



a = (7 

TT \ 



(C.49) 



We will restrict ourselves to the study of dissipation caused by the ohmic bath of 
bending modes, which will prevail at low temperatures over the superohmic bath of 
compression modes. 
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1 layer 


t = 10 nm 


a 


7-10^ 
8 ■ 10^^ 
~ 0.03 


2.5 ■ 10^ 

3 ■ 10^3 

~ 3 ■ 10-5 



Table C.l: Estimate of the parameters of the bath, for 7 ~ 1 eV. 



C.5.2 Losses due to the ohmic bath. Temperature depen- 
dence 

An estimate is given of the properties of the bath, namely a, Umf and Uco, for 
two cases, the one of a single layer graphene sheet (whose thickness we take to be 
t ~ 3 A), and a stack of layers 10 nm thick, using the value 7 ~ 1 eV, in ta- 
ble ( C.l[ ). This value of 7 is a big overestimate, since it corresponds to the value 



of the coupling of TLS's embedded in a vibrating amorphous structure to the os- 
cillations, and as stated previously this is not our case, and the coupling will be 
reduced by some power of the factor (a/d) ~ 10-^. We thus conclude that a is ex- 
tremely small, and there is no appreciable renormalization of the tunneling amplitude, 
= Ag • (Ag/cuco)"/^^-") ~ Ag [LCD+87]. In ref. [Wei99] an extensive analysis of 
the correlation function C{t) = {az{t)az{0)) of a TLS in presence of an ohmic bath is 
carried out, and for the cases we consider (a -C 1, biased TLS) the following results 
are obtained for the quasielastic (relaxational) peak around u = 0: 



1. For e < TTakT: 



2. For TiakT < e < kT: 



with 



3. For e > kT: 



kT A^ 



aeK-) = 2 ,1 , - 4e-^/^^(Ag)Ve^^^ , (C.53) 

COsh^f^)^ +7r ^^+1? 
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with 

7r = TTO; 



coth (^^) Al/he ^ TtaAl/he (C.54) 
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These expressions for Creii^^) substitute the one found usually in the context of sound 
attenuation in amorphous sohds due to 3D acoustic modes, C^I^{uj) = ^^^^2°^^^2fcT) t^^+r^ ? 
and they will enter into the expression for the dissipation as follows: 

Q-\uj,T) = ^^l del du^==ujC{u) (C.55) 



EkT Jn , uy/l 



u 



where u = \-/e, e^ax ~ 5 K, and {u\/l — u'^)~^ comes from the probability density 
of TLS's in an amorphous solid, like Si02. 

When kT > emax there will be TLS's only in the regimes of eqs. (C.50) and (C.51 ), and 
moreover, thanks to the very low value of a and the range of u typical of experiments 
we will have ■jr < u for all values of u, so that 7r/(t^^ + 7^) ~ Ir/^"^, and 



+ de du — 2(1 -M^) \ 

So for high temperatures it will give a constant contribution to dissipation, plus a 
weak (due to the prefactor a^) linear T dependence. 

In the case of lower temperatures, kT < emax-, there are some TLS's in the regime 
corresponding to eq. (C.53), and again one can approximate 7r/(ti^^ + ll) ~ Ir/^'^ 
thanks to the low a, leading to 



EkT i Jo J^^^^ uy/1 - 2TThujakT 

+ de du 2{l-u^) 

+ nde r du ^ 4(l-.^)e-A^^| 

P7' r47r Stt ^ 
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Thus, at low temperatures the contribution from TLS's to attenuation is linear in 
T. These tendencies come basically from the ohmic character of the vibrational bath 
and the probability density function of the TLS's. In case the TLS's resided in a 
crystalline material, or for some other reason their distribution function differed from 
the one of amorphous solids, changes in the behavior of Q~^{T, u) are to be expected. 
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D.l Dipole matrix element from single-particle mean- 
field states 

The derivation given in this section is borrowed from [Wei06] . We want to calculate 
the dipole matrix element dap defined in (6.14) 



d. 



(a 



za 



z9(a — r) H —Q{r — a 



(D.i: 



For that sake we need a description of the eigenstates |a) of the mean-field Hamilto- 



nian (6.11). The TDLDA numerical calculations briefiy presented in Sec. 6.2 nform 



us on the shape of the self-consistent mean field potential V{r), as shown in Fig. 6.1 



It suggests that for analytical calculations, V{r) can be approximated by a spherical 
well of radius a and finite height Vq, V{r) = Vo0(r — a), where Vq = ep + W, ep 
and W being the Fermi energy and the work function of the considered nanoparticle, 
respectively. 

Because of the spherical symmetry of the problem, the one-particle wave functions 



Ud(r 



r 



(D.2) 



decompose into radial and angular parts given by the spherical harmonics y^"^(^,0), 
where / and m are the angular momentum quantum numbers. The radial wave 
functions Ud{r) satisfy the reduced Schrodinger equation 



2mp dr"^ 



+ V{r) Ud{r) = eud{r) 



(D.3) 



with the conditions Ud{0) = and \imr-,oo[ud{r) / r] = 0. This yields the single- 
particle eigenenergies e in the mean-field potential V{r). Thus, the matrix elements 



147 



148 Appendix D: Appendix to Chapter^ 



dap can be separated into an angular part A^J^'^ and a radial part 7li^i^{€a,£i3), 

d,H = ATj;'ni^i,{ea,ep). (D.4) 



With (D.l), we have 

'.27r 



PTT PZ7V 

ATj;'= dOsine t/0y-'^*(^,0)cos^F,7(^,0) (D.5) 
Jo Jo 

which can be expressed in terms of Wigner-3j symbols [Edm60] as 

In Ir \\ I In la 1 



The Wigner-3j symbols contain the dipole selection rules = =t 1 and = m^^. 

Concerning 7li^i^{ea,£f3), in the limit of strong electronic confinement Vq ^ e^, 
the single-particle wave function vanishes outside the nanoparticle, enabling us to 
take the following approximation 

nup{ea,ef3) - [ dru*^ijr)ru,^i^^{r) (D.7) 

^0 

The commutator of r with the self-consistent Hamiltonian Hsc = p"^ /2me + V{r) reads 

[r,Hs,] = —Pr, (D.8) 
me 

where pr is the conjugated momentum to the variable r. Calculating the matrix 
element of this commutator between two single-particle eigenstates \a) and \(3) of Hsc 
then allows to write 

{a\r\f3) = —j^ Aa\pr\P) . (D.9) 

Now we use \pr,Hsc\ = —ihdV/dr to obtain 

{a\Pr\P) = —{a&P), (D.IO) 



and substituting in (D.9) one gets 



so that the radial part of the dipole matrix element can be reexpressed as a function 
of dV/dr as 

h"^ r dV 

nu.iea^s,) ^ I drul,^{r)-u,^,^{r) , (D.12) 
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which for an approximate step potential such as ours, V{r) = Vo0(r — a), leads to 



(D.13) 



For Vo — > oo and r < a, the regular solutions of (D.3) satisfying Uei{a) = are given 
by 

V2 rji{kr) 

where ju{z) are spherical Bessel functions of the first kind [A1S70], and k = \j2met/h. 
The condition ji{ka) = yields the quantization of the single-particle states in the 

he radial wave func 



infinitely deep spherical well. Inserting now the radial wave function (D.14) into 



(|Dl3|, 



mea{e/3 - ji^+i (k^a) j/^+i (k^a) ' 



(D.15) 



Now, for r —>■ a"*", Ud{r) — > and V{r) = Vq oo, so that the radial Schrodinger 
equation (D.3) takes the form 



2mp dr"^ 



Uei{r) = , 



With ko = v^2m7K)/ h, the physical solution of this equation is 

-kor 



Then 



dUf 



dr 



-koUdir — > a'^ 



(D.16) 



(D.17) 



(D.18) 



and imposing the continuity of the radial wave function and of its derivative with 
respect to r at r = a, we obtain 



du. 



dr 



-koUdir ^ a ) . 



(D.19) 



Therefore, using (D.14) 



lim 

Vb— ►oo 



V^meVo 

h 



ji{ka) 



dji{kr) 



dr 



djiiv) 



r=a dr) 



ri=ka 



(D.20) 



Using the recurrence relation [A1S70] 

I 



V 



-jiiv) 



drj 



ji+iiv) 



(D.2i: 
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and the fact that ji{ri) = ji{ka) 
element [YB92] 



Note that the dipole selection rule la 



0, we reach the final result for the radial matrix 



- Ij3 ± 1 appearing in the angluar matrix 



element (D.6) implies that Sa 7^ Sp, and thus (D.22) does not diverge. Furthermore, 



it decreases with increasing energy — ^/jI of the dipole transition |a) 



D.2 Lowest energy of the particle- hole spectrum 



Using the large ka expansion of (D.14) (semiclassical high-energy limit), the quan- 
tisation condition reads 

ka = ^ + (D.23) 

with / and n non-negative integers. The energy of a single-particle (hole) state is re- 
lated to its wavevector kp(k.), its total angular momentum Ipi^h)-, and its radial quantum 
number Up^h) as 



2mp 



TT 

k-pa 



p{h) 



+ '"'p(/i) 



Thus, the energy of a p-h excitation entering the RPA sum (6.48) is 

\^ f Ip ~ h , \ f lp + Ih , 

k^a (— + -.-Mx(— + -. + -'^ 



^£ph — 



(D.24) 



(D.25) 



Notice that using the exact quantum mechanical spectrum in (6.48) would not change 



significantly the result depicted in Figure 6.4, since the approximation (D.23) is very 



reliable for states close to the Fermi energy. We have also checked that generating the 



p-h excitation energies randomly in the RPA sum (6.43) does not affect the physical 



picture of Figure |6.4 Indeed, the main ingredient to understand such a picture is the 
fast decay of the dipole matrix element with the p-h energy. 



The expressions for the dipole matrix element (6.47) as well as for the typical 



dipole matrix element (|6.57|) diverge in the limit of a small p-h energy. It is therefore 

-mm that 



crucial for our analysis to determine the appropriate minimal p-h energy Asr 
renders this divergence unphysical. 



This can be achieved by imposing the dipole selection rules in (D.25) and that 



the energy difference is minimal. The first condition dictates that 



/p ± 1 and 



rrih = m. 



V 



Therefore there are two ways of obtaining the minimal energy difference: 



Up = rih with Ip 



//i + 1 and Up = Uh + 1 with Ip = Ih — 1- In both cases we have 



Ep kh 
kpa/7T kp 



(D.26) 
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Since we are interested in states close to the Fermi level, we can simplify (D.26) to 
expression (6.49). 

If we consider sodium clusters with k-pa = 30 (a = 3.3 nm and ^ 4000 con- 
duction electrons per spin direction), we have /^.e^in ~ £f/10. This is a much larger 
energy than the lowest one we can observe in the numerically generated excitation 
spectrum (see Fig. I in Ref. [MWJ02]). However, the two results are reconciled once 
we take into account the large degeneracy yielded by our approximate quantisation 



condition (D.23). 



The degeneracy of p-h excitations with minimal energy is given by twice the 
number of pairs (Ih, Uh) compatible with kh = kp and Asph = A^min- Indeed, we 
have seen that there are two possible particle states p starting from h and verifying 
the above-mentioned conditions. For each n between I and kpa/ir, there is a value of 
/ = 2{kpa/7r — n) and therefore the number of degenerate p-h excitations with energy 

A^min is 



A/- = 2 5^(2/ + l) = 2 ^ 



n=0 



n=0 



4(^-nl+l 



TT 



kpa 



TT 



(D.27) 



This degeneracy factor has to be included in Figure |6.4[ and it is crucial for the 
determination of the collective excitation. 



D.3 Density of states at a fixed angular momen- 
tum /, pi{e] 



For our derivations we need to have an estimate for the /-fixed density of states, 
which we will call pi{e). It corresponds to the number of eigenstates per energy unit 
whose radial part Ukiir) satisfy (D.3) with the conditions Md(0) = and Yrnir^ooiueiiT) I f] 
0. It is important to notice that the variable r is limited to positive values and that 
the centrifugal potential possesses a singularity at r = 0. This prevents us from a 
naive application of the WKB approximation [Mer70] to treat this radial problem, 
as would be done for a standard ID Schrodinger equation. First we have to rewrite 



(D.3) in terms of a variable where such a singularity is absent, and then we will apply 



WKB. The change of variables is 



X 



Inr, , Xki{x) = e^''^Uki{r) 



and (D.3) becomes non-singular 



with 



dx"^ 
2m, 



+ n2(x) Xki{x)=Q. 



l\2 

2, 



(D.28) 

(D.29) 
(D.30) 




Figure D.l: Left: Phase-space trajectory of a periodic motion of energy e in a con- 
fining potential. Right: Radial potential V^'^^-^ for fixed /, of the spherical cavity as a 
function of the radial coordinate r. 



Using now the WKB approximation for Xki{x) amounts to change the centrifugal 

(D.31) 



potential in (D.3) according to the Langer modification [Lan37, BM72] 

1\2 



+ (/ + -) . 

Therefore the modified effective radial potential becomes 



h\l + 1/2)- 
2mpr2 



V{r) 



(D.32) 



where in our case V{r) is the spherical billiard potential 



V{r) 



0, r < a , 
cxD, r > a , 



(D.33) 



A plot of V^^^{r) is provided in fig. (D.l). A semiclassical theory can be developed 



to calculate pi{e), obtaining an expression whose first corrections lead to an extra 



oscillatory component p"' 
Wei06], 



added to the Thomas- Fermi result for pi{e) [WMWJ05, 



dpdr6{e — Hi{r,p)) 



(D.34) 



where Hi{r,p) is the hamiltonian with potential (D.32). But for our purposes in this 
thesis the Thomas-Fermi result suffices, because the quantities we want to calculate 
will be integrals whose integrands include products of the kind Pi{£)pi±i{£), where 
the effects of the oscillating corrections are strongly suppressed. 



To calculate (D.34) we start showing that the phase space integral on the right 



hand side equals the period of the classical trajectory of energy e of the particle 
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confined in the V^^^{r) shown in fig.(D.l) 



j dpdr6{e - Hi{r,p)) = Ti{e) 



(D.35) 



The delta function can be expressed as 



(r' ,p') such that 

H{r',p')=e 



S(r — r',p — p') 
\VH{r\p')\ 



(D.36) 



where 5{r — r',p — p') fulfills J dpdr6{r — r',p — p')f{r,p) = f{r',p'), for a generic 
function f{r,p). Therefore the surface integral (D.35) is converted by the sum over 



delta functions into a line integral over the trajectory in phase space whose energy 
is fixed at the value e. The denominator can be reexpressed with the aid of the 
Hamilton equations satisfied by the classical trajectory: 



dp 

'dt 



dH 
dr 



dr dH 
dt dp 



(D.37) 



so that \VH {r',p') \ = ^[dH/drf + [dH/dpf becomes ^dr"^ + dp'^/dt = dl/dt. Sub- 
stituting in (D.35): 



/ dpdr 



{r ,p ) such that 

H{r\p')=e 



5{r — r',p — p') 
\VH{r',p')\ 



dl 



dl/dt 



dt = Ti{e) . 



(D.38) 



The next step is to calculate ti (e) . For that sake we use the classical action functional 

'•r{t=T,) 



2m 



dt + eTi 



dr ■ p . 



(D.39) 



r(t=0) 



which satisfies dS/de = ti. We will calculate S using its last form in terms of r and 
p, recalling the fact that for a trajectory of fixed energy 



p = sj2m{e~Vi^f^{r)) 



The value of S can then be calculated as 

'•r(t=T) 



dr-p = 2 / drJ2mie-Vi^^^{r)). 

r(t=0) Jr.{e) 



(D.40) 



(D.41) 



The value of the turning point r^{e) on the left of the classically allowed reg ion i s 
given by the condition V^^^ {r^{e)) = e, i.e., r_{e) = h{l + 1 / 2) / y^2rn^ , see fig.(D.l). 
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Substituting the form of the potential Vj^-^-^i^r) in the integral (D.41), after a straigh- 
forward calculation we get 



(D.42) 



and the period ti is obtained performing the partial derivative with respect to e, 
obtaining in this way the desired Thomas-Fermi estimate for the partial density of 
states: 



ri(e 



{2mea^/h?) - {l + \ 



(D.43) 



D.4 Local density of the dipole matrix element 



Equation ( 6.51[ ) defines the local density of dipole matrix elements connecting 
states at energies e and e + Ae. We used particle and hole states in our definition, 
since this is the main interest of our work. But note that the calculation presented in 
this appendix is not restricted to that case and can be easily extended to any states. 
The result would be of course unchanged. 

Local densities of matrix elements of arbitrary operators have been thoroughly 
studied as they can be easily connected with physical properties, ranging from far- 
infrared absorption in small particles [MR97] to electronic lifetimes of quantum dots 
[GJS04]. A semiclassical theory for the local density of matrix elements has been 



developed [FP86, Wil87, EFMW92], where (6.51) can be expressed as a smooth part 
given by correlations along classical trajectories plus a periodic orbit expansion. We 
will not follow here this general procedure, but continue to use the simple form of 



the dipole matrix elements (6.44) for states confined in a hard-wall sphere and the 



semiclassical approximation applied to the radial (fixed /) problem [MWJ02,MWJ03, 
WMWJ05]. 



(D.4, D.6, D.22) 



Using the /-fixed density of states (D.43) we can write with the help of equations 



poo pco \ 2 / \ 2 

^ lama ^ lama 



J] p,, (6)p,, (6 + Ae) (a"") ' [Ru, (6., e,)) ' 



(D.44) 



lamajlpmp 
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A^^™'' 7^ only if rria = mp and /q, = //j ± 1, therefore 



C(e,Ae) 



2n2\2e(e + Ae) 



E PiSe)pi,{e + ^e){ATX-y (D.45) 



We first perform the sum over Ip: 



Za=Za±l 



-Pi„+i(e + Ae) + — — — -— — — — pi^_i(e + Ae(p.46) 



(2/„ + l)(2/„ + 3) 
and then over mo,: 

E (/a + l-m«)(/« + l + 



77T.Q. /q: 



la 



(2/„-l)(2/„ + l) 

(/„ + l)(2/„ + l)(2/, + 3), 
L(2/,-l)(2Z, + l), (D.47) 



yielding 

C(£,Ae) 



2^2 \^ ee' 



E f^z.(^) + l)f^z.+i(e') + lhQi,-i{e')] (D.48) 



where /^ax is the maximum allowed Ih for an energy £, while e' = e + Ae. In the 
semiclassical limit we can take Z/i ~ Z/t + 1 3> 1 and convert the sum into an integral. 
Thus, 



~ 3\rna) Ae^ 27re 27r(e + Ae) 

Jo 



1 /2h\2 1 n 



Qn'^ \ maJ Ae^ 



(D.49) 



where y"^ = 2mea'^ /It? and y'^ = + 2mAea^//i^. Changing variables now to u = Z^, 
du = 2ldl, we get 



C(e,Ae) = (- 



1271^ \ma J Ae^ 
Using the results of ref. [GROO] we obtain 



du\J y"^ — u\J y'^ — u 



du\J y'^ — u^/y'"^ — u 



(^'2 _ y> _y 



In 



y' + y 



{y'^ + y^)yy' 



(D.50) 



(D.51) 



And substituting the values of y and y' one arrives at the desired expression (6.52). 
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D.5 Density of particle- hole excitations 



The density of p-h excitations with energy Ae is defined in (6.55) and can be 
written as 



deh X](2^/^ + ^)^iHi^h) X [gi^+iieh + Ae) + gi^^i{eh + Ae)] . 

(D.52) 



Using the semiclassical density of states (D.43) and performing the sum in the hmit 
//j 1, we obtain 



Ae"^ ( 2m^a'^ 



2 \ 2 



Fjeh/Ae) 



de 



(D.53) 



where the function F has been defined in (6.53). Performing the remaining integral 



over the hole energy in the limit Ae <^ ep is straightforward and leads to the result 
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